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The  polarisation  null  and  maximum  characteristics  of  simple  compound  targa 
are  studied;  with  the  objective  of  assessing  the  applicability  of  these 
concepts  to  target  discrimination.  The  polarisation  characteristics  of 
single  scattering  centers  ere  derived  in  theoretical  detail  and  shown  to  bi 
highly  distinctive  according  to  scattering  center  type.  However,  the 
dynamic  behavior  of  the  characteristics  with  changing  aspect  for  compound 
targets  is  shown  to  be  very  complex.  It  is  concluded  that  these  concepts 
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are  net  generally  useful  for  target  discriminations  if  an  extended  target 
is  viewed  at  a  single  frequency.  Extensions  are  suggested  to  make  the 
approach  useful./k 
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SECTION  1 


INTRODUCTION 

This  is  the  final  report  on  work  performed  under  Contract  No. 
F3060Z-81-C-0254,  August  1981  to  August  1982,  with  the  Rome  Air  Development 
Center,  Rome,  NY. 

The  objective  of  this  effort  was  to  study  the  polarization  null 
patterns  and  maximum  patterns  of  simple  aircraft  models,  at  practical  radar 
frequencies,  and  to  assess  the  feasibility  of  using  the  nulls  or  maxima  for 
target  classif icatio.i  in  both  benign  and  hostile  electronic  environments. 

The  approach  to  this  problem  included  combined  theoretical  and 
experimental  studies.  The  null  and  maximum  characteristics  were  computed  for 
theoretical  target  models  of  increasing  complexity:  single  scattering  centers 
of  various  types,  combinations  of  scattering  centers,  and  complete  modeled 
aircraft.  This  study  provided  insight  into  the  behavior  of  polarization 
characteristics  as  functions  of  aspect,  frequency,  and  target  geometry. 
Considerable  effort  was  expended  to  assure  realistic  modeling  of  the  target 
returns.  Measurements  were  performed  using  a  dual-polarization  35  GHz  radar. 
Three  aircraft  models  of  slightly  differing  shapes  and  each  at  four  different 
size  scales  were  measured.  Measurements  were  made  over  0*-90*  aspect 
variation  and  at  0*  and  10*  elevation  angles.  The  data  is  presented  in  the 
form  of  scattered  Stokes  vectors,  from  which  the  nulls  and  maxima  are 
derived.  Agreement  between  computed  and  measured  polarization  characteristics 
was  good,  so  that  conclusions  could  be  drawn  with  confidence  from  the  more 
readily  manipulated  computed  characteristics. 

This  study  was  limited  to  the  monostatic  case,  which  results  in  a 
significant  simplification  due  to  the  symmetry  of  the  scattering  matrix.  The 
bi static  case  is  theoretically  more  complicated  and  will  not  exhibit  some  of 
the  symmetry  relations  between  null  and  max  polarizations  derived  here  for  the 
monostatic  case.  However,  the  dynamic  behavior  of  nulls  and  maxima  can  be 
expected  to  be  similar. 
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The  study  was  also  originally  limited  to  single  frequency 
illumination.  It  was  found  that  this  resulted  in  a  limitation  on  the 
usefulness  of  polarization  characteristics  for  target  identification,  for 
targets  distributed  over  many  wavelengths.  Hence  some  tentative  extensions 
were  mads  to  multifrequency  illumination.  More  work  in  this  area  is  expected 
to  yield  a  highly  effective  target  discrimination  technique. 

It  is  known,  and  is  detailed  in  this  report,  that  scattering  centers 

of  different  types  (e.g.  plates,  dihedrals,  edges)  have  highly  distinctive 

polarization  characteristics.  It  is  therefore  possible  to  infer  sea tterer type 

from  a  set  of  radar  returns  at  judiciously  chosen  transmit  and  receive 

polarizations.  Combinations  of  uncoupled  scattering  centers  (such  as  an 

aircraft  at  practical  radar  frequencies)  also  have  distinctive 

characteristics,  dependent  on  the  characteristics  of  the  parts  and  on  the 

spatial  relations  between  them.  It  might  be  expected  that  two  aircraft,  which 

are  similar  in  overall  appearance  but  differ  in  detail,  can  exhibit 

sufficiently  distinctive  polarization  behavior  to  permit  discrimination.  This 

study  shows  that  this  is  indeed  the  case:  at  any  particular  aspect  and 

frequency,  two  such  aircraft  do  exhibit  different  polarization 

characteristics.  However,  the  dynamic  behavior  with  aspect  or  frequency  of 

the  characteristics  is  so  complicated  that  their  usefulness  for  target 

discrimination  is  problematic  without  further  processing.  This  is  the  major 

result  of  this  study:  the  instantaneous  polarization  null  and  maximum 

behavior  of  compound  targets  is  too  complex  to  permit  effective  target 

discrimination  based  o.i  single  frequency/ aspect  data,  except  possibly  over 

certain  (narrow)  aspect  bins.  In  particular  it  is  shown  that  none  of  the 

polarization  characteristics  are  invariant  with  aspect,  as  has  at  times  been 

speculated  in  the  literature.  The  difficulty  is  very  similar  to  (and,  in 

fact,  has  the  same  cause  as)  the  rapid  fluctuations  of  RCS  versus  aspect  - 

while  the  patterns  may  be  distinctive,  they  are  too  complex  for  direct  use  in 

a  classifier.  There  are  two  ways  to  overcome  the  difficulties  of  rapidly 

fluctuating  null  patterns:  either  by  averaging  over  aspect/frequency  diverse 

data;  or  by  utilizing  the  measured  fluctuations  to  resolve  parts  of  the 

target.  The  simpler  averaging  method  is  discussed  to  some  extent  In  this 

report.  It  Is  seen  that  discrimination  may  be  possible  at  those  aspects  where 

a  distinctive  scattering  center  is  dominant.  The  second  method  is  also 
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considered  to  show  great  promise.  It  *s  similar  to  range  or  cross-range 
profile  methods,  but  with  the  single  amplitude  in  a  resolution  cell  replaced 
by  a  multi-dimensional  polarization  characteristic.  Further  discussion  of 
this  method  is  beyond  the  scope  of  this  report;  it  is  mentioned  here  to.  soften 
the  blow  of  the  rather  negative  conclusions  drawn  here  regarding  target 
discrimination  under  restrictive  operating  conditions. 

The  polarization  characteristics  upon  which  most  emphasis  Is  placed 
in  this  report  are  the  co-  and  cross-polarized  null  and  maximum  polarizations 
and  their  respective  powers.  In  particular,  there  Is  a  set  of  five  real 
parameters,  m,  Zyt  2v,  2t  ,  and  2*  ,  an  amplitude  factor  and  four  angles, 
in  terms  of  which  these  polarizations  and  powers  are  defined.  The  usefulness 
of  these  five  parameters  lies  in  that  they  form  a  set  of  independent  and 
complete  target  descriptors,  which  are  invariant  under  polarization 
transformations.  That  is,  the  five  parameters  are  characteristic  features  of 
the  target  at  a  particular  aspect  and  frequency  and  are  independent  of  the 
transmit  and  receive  polarizations  used  to  measure  them.  They  are  also 
invariant  under  target  roll  about  the  RLOS.  However,  they  are  not  invariant 
under  more  general  target  orientation  changes  (aspect),  nor  under  frequency 
changes.  The  completeness  of  the  five  parameters  holds  only  for  static, 
deterministic  targets,  for  which  the  scattered  radiation  is  fully  polarized. 
For  time-variant  and  random  targets  (meaning,  really,  that  the  relative 
orientations  of  the  scattering  centers  of  a  distributed  target,  such  as  chaff, 
are  changing  during  the  measurement  interval),  the  scattered  radiation  becomes 
depolarized,  requiring  four  additional  parameters.  This  is  related  to  the 
fact  that  five  independent  real  numbers  characterize  the  relative  phase 
scattering  matrix,  which  fully  describes  coherent  scattering,  while  nine 
independent  real  numbers  characterize  the  Nueller  matrix  which  fully  describes 
incoherent  scattering. 

Only  three  parameters,  namely  m,  2r»  2v,  are  required  to 
characterize  a  plane-symmetric  coherent  target  (with  the  RLOS  In  the  plane  of 
symmetry).  The  aircraft  models  studied  here  at  0*  elevation  are  of  this 
type.  This  case  is  studied  In  detail,  because  of  the  more  readily  assimilable 
insight  it  provides. 
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The  polarization  nulls  and  maxima  were  not  measured  directly. 

Rather,  two  (fixed)  orthogonal  transmit  polarizations  were  used.  At  each 
transmit  polarization,  the  received  power  was  measured  for  each  of  six 
different  receiver  polarizations.  (Only  one  transmit  polarization  is  required 
for  the  plane  symmetric  cases.)  The  procedure  and  data  recording  was  done 
under  computer  control.  From  the  power  measurements,  the  Stokes  vectors  and 
polarization  characteristics  are  derived.  Great  care  was  taken  to  develop  a 
good  polarization  calibration  procedure. 

In  Section  2,  the  theoretical  background  is  presented,  including  a 
discussion  of  the  characteristics  of  simple  compound  targets.  In  Section  3, 
the  computed  polarization  characteristics  are  presented  for  each  of  the 
aircraft  models.  The  measurement  system  and  measurement  results  are  given  in 
Section  4.  In  Section  5,  the  effects  of  noise,  chaff,  and  jamming  on 
polarization  null  and  maxima  concepts  are  discussed.  The  use  of  polarization 
agility  as  an  ECCM  to  reduce  the  effect  of  jamming  is  analyzed.  A  summary  is 
given  in  Section  6;  followed  by  References  and  an  Appendix  containing  some 
mathematical  details. 

Wolfgang-M.  Boerner  gives  an  excellent  history  and  extensive 
bibliography  of  polarimetry  (See  Appendix,  page  8.12). 


SECTION  2 


CHARACTERIZATION  OF  POLARIZATION  PROPERTIES  OF  SIMPLE  TARGETS 


In  this  section  definitions  and  standards  of  polarimetry  are 
presented*  and  the  polarization  properties  of  elementary  and  simple  compound 
targets  are  investigated.  The  literature  abounds  with  conflicting  and 
sometimes  self-inconsistent  definitions  of  the  descriptive  parameters  of 
polarization  phenomena.  An  effort  has  been  made  here  to  present  a  consistent 
set  of  definitions  which  is  compatible  with  other  engineering  standards.  The 
results  are  almost  identical  to  the  work  of  Chan  and  Boerner  [1],  although 
there  are  some  notatlonal  differences,  due  to  the  independent  history  of 
subject  development  at  Sperry  Research  Center  (SRC). 

Much  of  the  material  in  this  section  is  not  new,  being  a  compendium 
of  earlier  development.  However,  emphasis  will  be  put  here  on  the  concepts 
which  are  useful  to  this  report.  Section  2  concludes  with  a  detailed 
investigation  of  the  polarization  properties  of  dumbells  and  cylinders.  Those 
results  are  extremely  useful;  first,  for  developing  a  feeling  and 
understanding  for  polarization  concepts;  and  second,  for  the  direct 
applicability  to  more  complicated  targets. 

The  polarization  properties  of  radar  targets,  utilizing  null  and 
maxima  concepts  ("optimal  polarizations"),  were  investigated  by  Kennaugh  [2] 

In  the  late  1940's  and  early  1950's,  and  received  comprehensive  treatment  by 
Huynen  In  the  1950's  [7]  and  in  his  1970  Thesis  [3].  More  recently,  Poelman 
has  built  an  experimental  polarlmetrlc  radar  and  made  progress  in  using 
polarization  properties  to  Increase  the  detection  capability  of  radar  systems 
[17].  (See  References  [16,4]  and  Volume  II  of  this  report  [18]  for  an 
extensive  collection  of  references  and  historical  notes.)  Some  of  the 
standards  used  in  this  report  differ  from  the  above  works  for  the  sake  of  ease 
and  compatibility  with  engineering  practice,  as  will  be  noted. 


2.1  REPRESENTATION  OF  POLARIZATION  STATES 


For  a  plane  TEH  wave  propagating  In  the  +2  direction,  we  write 

F(z,t)  -  Re  |  Eq  E  ej<“t-kz+*>  J  ,  (1) 

where  the  (complex)  unit  polarization  vector  Is 


2  2 

with  |  ex  |  +  |ey  |  ■  1.  We  will  associate  x  and  y  with  horizontal 

(H)  and  vertical  (V)  polarizations.  For  linear  polarization,  we  have  «  -  0  , 
and  we  can  define  the  rotation  Y  In  the  x-y  plane:  tan  ¥  -  |e  |/|e  |  . 

Jr  A 

For  circular  polarization,  we  have  |  ex  |  -  |ey  |  «  1/  and  «  »  ±  »/ 2. 

We  choose  «  -  -  w/2  for  right  circular  (RC)  and  a  *  »/2  for  left 
circular  (LC).  That  Is: 


This  convention  agrees  with  engineering  usage  and  Is  contrary  to  that  of  Hqynen 
and  of  physics  texts.  This  convention  corresponds  to  the  sense  of  rotation  of 
the  E  vector  in  time  at  a  fixed  point  In  space.  The  general  elliptical 
polarization  of  Equation  (2)  can  be  related  to  the  motion  of  a  point  on  the 
perimeter  of  the  polarization  ellipse  illustrated  In  Figure  1.  The 
Illustration  Is  for  right-sense  polarization,  for  which  the  elllptlclty 
angle  t  Is  defined  to  be  negative.  For  linear  polarization  t  -  0  ;  and  for 
left-sense  polarization  t  >  0  .  For  consistency,  we  require  the  restrictions 
-t/4  <  t  <  w/4  and  -  */2  <  't  <  w/2.  Using  t  ,  f  we  can  represent 


the  phase  a  Is  adjusted  so  that  the  phase  of  ex  Is  zero  as  In  Equation  (2). 
The  reason  for  this  form  Mill  become  clear  In  the  next  section  Mhen 
transformations  are  discussed. 


The  normalized  polarization  vector  can  be  represented  In  terms  of  power 
measurements  as  the  normalized  Stokes  vector  s  : 

lexl2+ley|2  s0 

|  ex! 2  -Kl2  Si 

s  -  x  y  -  1  (4) 

2  Re  {exey*}  s? 

-  2  im  {exey*}  S3 


(The  above  definition  of  s  applies  to  a  right-handed  system  y,  k.  Use 

of  the  ruonostatlc  convention  (see  Section  2.3  and  Appendix  A--3)  also  Introduces 

a  left-hand  system.  In  a  left-hand  system,  s^  has  the  opposite  sign.)  For 

2  J  2  2? 

the  fully  polarized  wave  of  this  section  Sq  «  1  ■  s![  ♦  s£  +  S3 

which  does  not  hold  for  partial  polarization.  Evaluating  the  Stokes  vector 

for  the  forms  (2)  and  (3),  we  obtain 


COS  ?T  COS  2f 
cos  ;!t  sin  2 f  , 

sin  2t 

L  J(5) 

as  is  shown  by  Huynen,  but  with  a  different  order  of  the  s^  . 

The  right-handed  triplet  (s^SgtSj)  Is  conveniently  plotted  as  a 
point  on  the  surface  of  the  normalized  Poincare  sphere,  as  In  Figure  2.  The 
Cartesian  coordinates  on  the  sphere  are  (s^^.Sj)  and  the  i*olar 
coordinates  (r,©,4)  are  (1,  »/2  -  2r  ,  2T  ).  Orthogonal  polarizations  are 
antipodal  on  the  sphere,  where  and  are  orthogonal  If  fj* -  0. 

Note  that  two  polarization  vectors  which  differ  only  by  a  relative  phase  (e^* 
In  Equation  (3)),  are  mapped  Into  the  same  Stokes  vector  *  .  Table  1 
summarizes  the  coordinates  of  some  basic  polarizations. 
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LINEAR  48° 


FIG.  1  Normalized  polarization  ellipse  for  right-sense  polarization. 


FIG.  2  Normalized  Poincari  sphere. 


POLARIZATION  COORDINATES 


2.2  POLARIZATION  TRANSFORMATIONS 


The  transformation  from  one  polarization  basis  to  another,  «  QE^, 
can  be  (and  has  been)  written  In  many  forms.  Most  convenient  Is  to  write  Q 
In  the  form  _ 


Q  - 


<2 

* 

ql 


» 


(6) 


a  unitary  transformation  with  the  property  Q*QT  -  I  and 
Det(Q)  -  (qjj2  +  Iq^2  •  1  (See  Appendix  8.1  for  some  details  and 
generalizations.)  If  two  polarizations  are  considered  equivalent  when  they 
differ  only  by  a  relative  phase,  then  there  exists  a  unique  Q  of  the  above 
form,  which  takes  one  orthonormal  basis  pair  of  vectors  Into  another.  In 
particular  we  have 


(7) 


for  transforming  from  a  linear  to  a  circular  basis. 


where 


with 


(8) 


ke  see  that  right-circular  polarization  Is  transformed  as 


and  left  circular  as 
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C  .responding  to  a  transformation  Q  Is  a  4  x  4  transformation*  V, 
of  the  Stokes  vector.  This  Is  worked  out  In  Appendix  8.2.  Reference  to 
Figure  2  of  the  Poincare  sphere  suggests  the  decomposition  of  Q  Into  a 
product  of  three  transformations 

Q  -  Qf  QT  Qv  »  i») 

corresponding  to  the  Euler  rotations  2v,  w/2  -  2t  ,  and  2Y  of  the  Poincare 
sphere,  with  corresponding 

V  -  vt  vT  Vy  .  (10) 

These  matrices  are  derived  and  tabulated  In  Appendix  8.2.  The  angles  t  and  T 
have  already  been  introduced  in  Section  2.1  as  descriptors  of  the  polarization 
vector;  namely  elllptlclty  and  spatial  rotation.  The  angle  v  ,  called  the 
"skip  angle"  by  Huynen.  will  be  seen  to  be  more  Intrinsically  associated  with 
the  target  scattering  matrix.  Q  changes  the  relative  phase  between  the  two 

A  V 

components  of  E  ,  without  changing  their  magnitude.  The  subset  of  unitary 
transformations  containing  only  HI  and  x  will  be  called  "polarization 
transformations". 


2.3  THE  POLARIZATION  SCATTERING  MATRIX.  MUELLER  MATRIX,  AND  CANONICAL  FORMS 


The  general  polarization  scattering  matrix  In  a  linear  basis  Is  written 


(ID 


for  the  monostatic  case  a^  -  «yH  (see  Appendix  8.3  and  [5]).  Further,  we 
will  always  deal  with  the  relative  phase  scattering  matrix  so  that  we  can  write 


f  I  *H  I  6  ”  l*HV  I  \ 
A  -  jPv 

\\\rt  IM#  J 


(12) 


which  consists  of  five  distinct  real  numbers.  (This  matrix  Is  usually  called  S.) 
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The  scattered  field  is 


F  s  -  A  E  *  .  (13) 

In  the  so-called  "monostatic  convention"  (Appendix  8.3),  the  H  and  V 

A  f 

components  of  E  and  E  have  the  same  polarity,  although  the 
propagation  vectors  are  opposite  directed  (Note  that  this  convention  makes  the 

AAA 

(H,V,k)  triplets  have  opposite  handedness!) .  With  this  convention  the 
received  voltage  is 


where  Er  is  the  receiver  polarization,  and  the  power  is  P  =  VV*  .  A 
convenient  way  to  write  this  is 


V  = 


rTAt 


(14) 


where  t  and  r  are  the  transmit  and  receive  unit  polarization  vectors, 
respectively.  By  reciprocity,  a  single  antenna  receives  the  same  polarization 
it  transmits  (r  *  t). 


Under  a  change  of  basis,  say  t  -  Qt *  and  r  »  Qr'  ,  we  have 


rTAt 


r  ,TQTAQt '  .  r,TA't\ 


with 


A' 


qtaq 


(15) 


a  m 

In  particular,  suppose  that  r  -  t  -  ER  -l,  I  ,  right  circular  in  a 
linear  basis,  and  that  A  •  I  ,  the  scattering  matrix  for  a  simple  specular 
return  (e.g.  a  sphere).  Using  linear  coordinates  we  obtain  V  *  [j  l][i]a 
as  expected,  since  a  single  bounce  reflector  changes  the  sense  of  circular 
polarization,  which  cannot  be  received  by  the  same  antenna.  In  a  circular 


0, 


basis,  we  have  r 
resulting  in  V  - 


Cl 


*'  “  «ie»  -(i)  ■  <>Ic«ic  '  ll‘(]  l)(j 

01  (°  c)Q-  “ o]  0)  -  0,  as  desired. 
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We  now  want  to  find  a  canonical  form  of  A  under  a  unitary 
transformation.  In  Appendix  8.4,  it  is  shown  that  A  can  be  diagonalized  by 
solving  the  eigenvalue  problem  Ae  ■  xe*,  obtaining 


S 


0 


ej8QTAQ 


(16) 


where  b  is  an  arbitrary  phase  angle  which  cancels  that  of  det(A),  and  where 
m  and  s  are  real.  The  significance  of  this  is  that  we  can  investigate  the 
behavior  of  the  simpler  form  Sq  and  obtain  the  behavior  of  A  by 
application  of  a  unitary  transformation.  As  shown  in  Appendix  8.4,  a  by-product 
of  the  derivation  of  (16)  is  the  observation  that  invariants  of  any 
transformation  Q  (when  A  is  symmetric)  are  the  norm  and  determinant  of  A, 


2  2  2  -i 

all  I  +  la22 1  +  2la12l  J 

det(A)  -  an  a22  -  a122 

(For  the  more  general  Qg  ,  which  may  have  det(QG)  -  e^2. 


where 


norm(A) 


4 


invariant 


we  can  absorb 


e^#  in  Equation  (17)  into  Qfl,  with  the  result  that  |det(A)|  is  invariant.) 

The  columns  of  Q  are  the  orthonormal  eigenvectors  e^.e,,.  (The 
orthogonality  is  obtained  only  if  A  is  symmetric.)  The  entries  m  and  s 
are  related  to  the  eigenvalues  by  m  -  |Xj  !,  s  -  |*2|,  and  where  we  choose 
m  to  correspond  to  the  largest  eigenvalue. 


Then  following  Huynen,  we  introduce  the  angle  y  ,  writing  SQ  in  the 


form 


S°  '  "(  2  ) 
\0  tan  Y  / 


,  0  <  U I  < 


i  T" 


(17) 
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It  is  also  possible  to  form  a  4  x  4  matrix,  M  ,  representing  scattering 
which  operates  on  the  Stokes  vector.  This  has  been  called  the  Mueller  matrix. 
Similar  to  Equation  (13)  we  write 


where  st  and  s$  are  the  Stokes  vectors  of  the  transmitted  and  scattered 
radiation  respectively.  (Since  these  have  opposite  handedness  in  accordance 
with  our  monostatic  convention,  remember  that  s3  is  defined  with  opposite 
sign  for  the  transmitted  and  scattered  vectors;  see  Appendix  8.3.)  The  received 
power  can  then  be  shown  to  be 

P  =  — srT  M  s^.  (18) 

The  receiver  Stokes  vector  is  defined  to  have  the  same  handedness  as  though 
the  receiver  were  operated  as  transmitter. 

For  the  corresponding  transformations  Q  ,  V  .  we  can  form  (as  '.n 
Equation  15), 

P  -  “5“  srT  M  st  *  *T"  srT  vTf1  v  st  -  -f-  srT  M'st 
with  M'  .  VTHV  .  (19) 

In  particular,  corresponding  to  the  canonical  form  S  ,  there  is  a  canonical 
form  MQ  .  This  can  readily  be  found  by  writing  out  the  Stokes  vectors  of 
t  and  s  in  s  ■  S  t  for  four  independent  cases  of  t  and  solving  for  M 
in  ss  «  MqSj.  .  The  result  is 
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This  can  also  be  written  in  the  form 


4  cos  y 


2 

1  +  cos  2y 

2  cos  2y 


2  cos  2y  0  0 

1  +  cos22y  0  0 

0  sin22Y  0 


In  particular,  the  unit  scattering  matrix  S, 


■  i  o 


0  -s1nz2 


corresponds  to 


(20  * ) 


(Note,  Huynen  includes  a  factor  of  1/2  in  his  definition  of  M  .)  The 
Mueller  matrix  is  symmetric  If  and  only  If  the  scattering  matrix  is 
symmetric.  A  more  elegant  way  to  relate  A  ,  M  ,  Q  ,  and  V  is  given  at 
the  end  of  Appendix  8.2 


2.4  POLARIZATION  NULLS  AND  MAXIMA 

The  received  voltage  can  now  be  written  as  (see  Equation  15) 

V  -  r‘TSo  t‘  , 

v  where  r',t*  are  the  polarization  vectors  In  the  transformed  basis  (we  have 
dropped  the  phase  term  e'**,  since  it  is  arbitrary).  The  co-polarized 


voltage  (r*  «  t 


Vc  -  t,TS  t*  .  t2  m  +  t|s. 


The  cross-polarized  voltage  (r,Tt'  «  0)  Is 


*  * 


Vx  “  tjt^m  - 


since  r*  is  orthogonal  to  t'  if  *  tj,  r^  ®  -t^.  We  can  now 
immediately  obtain  the  maxima  and  nulls.  Noting  that  |tjf  +  | t^ | ^ 
write, 

|vc|<  I ti! 2m  +  I l2l ?s  “  Iti^C^-s)  +  s  , 

where  the  equality  holds  if  t^  and  t2  have  the  same  phase.  Since  m  >  s, 
the  right  side  of  the  inequality  is  maximized  when  t^  =  1  ,  say 


t,  =  1, 


hence  | t 2 j 


t2  -o. 


Hence  the  co-pol  max  polarization  and  power  are 


t' 


ith  PCM  =  |Vf;M  |  -  m 


CM 


(23) 


Writing. 


Vq  j  <  m  -  |t2  f(m  -  s)  , 


we  see  that  the  right  side  is  minimized  for  | t?  |  =1,  say  t2  =  1,  but 
that  the  left  side  is  maximized  (i.e.  equality  holds)  if  tj  and  t2 
have  the  same  phase.  Therefore,  this  is  a  saddle  point  (maximized  subject  to 
constant  relative  magnitude,  but  simultaneously  minimized  subject  to  constant 
relative  phase).  Hence  the  co-pol  saddle  point  polarization  and  power  are 


V 


(:) 


with  P 


CS 


c2  2.4 
s  =  m  tan  y 


(24) 


Note  that  tM  and  t^  are  orthogonal,  and  that 

r»  o 

norm  (S)  «  norm  (A)  . 


PCM  +  PCS 


m  +  s 


(25) 


We  see  directly  from  Equation  (22)  that  for  each  of  these  polarizations  *  0 
The  co-pol  max  and  saddle  points  also  the  cross-pol  null  points. 


L 
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His. .» 5 


2  2 

For  the  co-pol  null,  we  solve  Vq  *  0  -  m  tgS, 
obtaining  the  two  solutions 

t  '  - - — 

"  /m  +  s 

Using  the  definition  of  y  in  Equation  (19),  notice  that,  significantly,  the 
co-pol  null  locations  contain  information  concerning  the  relative  magnitudes 
of  the  eigenvalues  of  A  ,  whereas  the  other  null  and  max  locations  are  points 
on  the  canonical  axes,  and  hence  only  serve  to  orient  the  diagonalizing 
transformation,  Q  . 

The  cross-pol  max  location  is  slightly  more  difficult  to  derive, 
although  the  result  is  simple. 


This  is  maximized  when  the  right-most  term  is  as  negative  as  possible.  This 
happens  when  t^/t^  -  *  j  and  |  t^ j  «  (  t2 |  «  1/  /2  .  Hence,  for  the 

x-pol  max  polarizations  and  power,  we  have 


(Note  that  Pxx  »  m*  at  |  y  } •  -J-  )  . 

Also  note  that  for  arbitrary  transmit  polarization,  the  total  power  received  by 
the  co-  and  cross-polarized  channels  of  the  receiver  is 

2  2  2  2 
P  "  lVd  +  *  0,2  <  I  +  tan4*  1 12 1  )  1  m2 
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These  "optimal  polarizations"  can  be  plotted  on  the  Poincare  sphere  In 
the  canonical  basis  using  definition  (4).  The  results  are  that  tM'  goes  into 
a  point  at  s^',  t$'  antipodal 1y  at  -s^',  the  two  t^  go  into  ±s3', 
while  th«?  co-pol  nulls  go  Into 


All  these  points  are  In  the  same  plane,  s^'  -  0;  that  is,  they  all  lie  on  a 
great  circle  of  the  Poincare  sphere.  This  is  illustrated  In  Figure  3  and 
summarized  in  Table  2.  The  points  are  unique,  except  for  two  degenerate  cases, 
which  are  discussed  below.  The  points  M,  Np  define  the  polarization 
(or  "Huynen")  fork.  The  diagram  for  the  scattering  matrix  A  is  obtained  by  an 
Euler  rotation  of  the  Poincare  sphere,  as  defined  by  the  transformation,  Q. 

The  maximum  power  received  is  Pmax  «  PM  »  m  ,  which  occurs  for 
cc-polarized  reception  at  t^.  An  effective  radar  cross  section  is  defined  by 
Kennaugh  [6]  (there  is  a  factor  of  4  error  in  [6])  as 


(1  +  tan2Y)2 


(28) 


A  subsit  of  the  above  quantities  completely  characterize  the  polarization 
properties  of  the  target.  To  obtain  the  scattering  matrix  using  power 
measurements,  it  is  necessary  to  establish  the  orientation  of  the  Huynen  fork, 
its  opening  angle  4-r  ,  and  the  scale  factor  m  .  Any  combination  can  be  used, 
for  example  t^ ' ,  tN',  and  P^  . 


(a)  NON -DEGENERATE  CASE 


DEGENERATE  CASES 


(e)  7  -  4T  «s*T 


PIG.  3  Great  circle  cut  (Sj  ■  0)  on  Poincart  *>here,  showing  nulls  and  maxima  for 
canonical  (positive  real  diagonal)  scattering  matrix. 


i 

t 

I 


I 

I 


TABLE  2 
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|v  I A 


There  are  two  degenerate  cases  to  this  analysis.  These  are  also  shown 
In  Figure  3^  The  first  Is  the  case  *  -  0  .  The  canonical  scattering  matrix 
Is  Sq  •  (  o  o)  *  *****  *s»  one  *he  e*9enva^ues  *ero.  The  co-pol 

nulls  collapse  onto  the  saddle  point,  and  the  x-pol  max  points  "flash"  Into  the 
great  circle  (s1  -  0)  orthogonal  to  the  Sj  axis,  with  ?xx  ■  m  /4.  If 
q  «  I  so  that  A  -  S0,  this  Is  the  case  of  a  horizontal  dipole  with  a  co-pol 
max  at  H,  co-pol  null  at  V,  and  x-pol  max  at  any  45  elliptical  polarization. 


The  second  degenerate  case  is 

matrix  is  S  -  m 
w/4 


generate  < 

(l !  ‘ 


*  t/4  .  The  canonical  scattering 

The  co-pol 


that  Is,  the  eigenvalues  are  equal. 

2 


nulls  coincide  with  the  x-pol  maxima  with  PyN  -  P 


maxima  "flash"  Into  the  great  circle  s. 


XN 

0 


XX 


m 


The  co-pol 


orthogonal  to  s^.  If  Q  «  I 


so  that  A  *  Sq  ,  this  is  the  case  of  (circularly  symmetric)  specular 
reflection;  with  a  co-pol  max  at  any  linear  polarization  and  x-pol  maxima  at 
right-  and  left-circular  polarizations. 


The  important  case  of  a  simple  double  bounce  scatterer:  A 


is  not  one  of  our  canonical  forms.  It  is  obtained  from  S 


74 


by 


*  *T 

A  -  Q  S  ,.Q 
v  w/4Mv 

with  .  n 

v  -  *  -jp 

•nls  corresponds  to  a  90*  rotation  of  Figure  3c  about  axis  s^'.  The  result 
s  that  the  perimeter  s2'  -  0  Is  now  the  great  circle  of  co-pol  maxima, 
corresponding  to  elliptical  polarization  with  ¥  »  0  ,  including  pure  H,  V,  RC, 
and  LC.  The  co-pol  nulls,  coincident  with  the  x-pol  maxima,  occur  at  the  *  45* 
linear  points. 
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2.5  NULLS  AND  MAXIMA  FOR  GENERAL  SYMMETRIC  TARGET 


In  order  to  gain  an  understanding  of  the  behavior  of  nulls  and  maxima, 
we  will  first  investigate  plane  symmetric  targets.  If  the  plane  of  symmetry 
Is  aligned  with  H-  (or  V-)  polarization,  then  the  scattering  matrix  Is 
diagonal.  (Note  that  If  the  target  Is  rotated  about  the  RLOS,  we  can  simply 
rotate  the  polarization  reference  plane  to  maintain  alignment  with  the  target 
symmetry  plane.  Hence,  If  by  H-  (or  V-)  polarization  we  always  mean  that 
Th  (orTy)  Is  aligned  with  the  symmetry  plane,  then  our  result  will  apply 
to  all  plane  symmetric  targets,  such  that  the  symmetry  plane  contains  the 
RLOS.)  The  general  form  of  the  scattering  matrix  Is  then  (see  Equation  (17))  : 


where 


A 


*  -  *T 
Q  SnQ  *  m 
v  0  v 


f 

eJ2v 

v° 


t  -  Q  t'  , 


(29) 


t  is  the  polarization  vector  In  the  H-V  basis,  and  t'  Is  the  polarization 
vector  in  the  canonical  basis. 


The  matrix  A  is  completely  characterized  by  the  magnitude,  m,  and 
the  two  angles  y  and  v  .  4y  was  Identified  In  Section  2.4  as  the  opening 
anglfc  of  the  Huynen  fork.  The  angle  v  was  named  the  “target  skip  angle"  by 
Huynen,  because  It  Is  related  to  how  closely  the  target  resembles  an  odd-  or 
even-bounce  scatterer:  If  v  ■  0  ,  the  target  is  single  (odd)  bounce;  if 
v«w/4  the  target  Is  double  (even)  bounce.  Note  that  v  Is  given  by  the 
relative  phase  of  a^  and  avv  :  4v  -  phase  (aHH)  -  phase  (avv); 
and  that  y  is  given  by  the  relative  magnitude:  tan2Y  -  |aw^aHH  I  *  In 
Equation  (29),  A  is  written  with  the  assumption  that  |  a^  |  >  |avv  |. 

(In  case  |aVy|>  1^1  .  we  should  reverse  the  meaning  of  H  and  V  and 
maintain  the  same  analysis). 
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From  Table  2  we  can  Immediately  obtain  the  optimal  polarizations, 

using  t  -  Q^t 1  .  In  particular,  (after  removing  the  arbitrary 

phase  factor),  with  the  Stokes  representation  (1,  1,  0,  0).  That  Is,  the 
co-pol  max  always  occurs  at  H  (or  V  If  |ayV  Ja^  |)  for  plane 
symmetric  targets.  (For  the  degenerate  case  ja^j-  |*yy|«  the  co-pol 
max  occurs  at  any  linear  polarization.) 

(e-JY  sin  y\ 

_ I*  The  Stokes  vector 


cos  y. 


The  Stokes  vector 


representation  can  be  obtained  either  by  applying  the  definition  (Equation 
(4))  or  by  transforming  from  Table  2  using  (Equation  (A-6)): 


-  cos  2y 
+  sin  2v  sin  2t 
*  cos  2v  sin  2v 


cos  2  tN  cos  2 
cos  2  tn  sin  2 
sin  2  t„ 


The  right  side  of  the  equality  comes  from  Equation  (5),  giving  the  Stokes 
vector  In  terms  of  the  elliptical  polarization  parameters  f ,  t.  Solving 
for  ,  t  ^  ,  we  can  plot  the  co-pol  nulls  on  the  polarization  chart, 
which  Is  a  projection  of  the  Poincare  sphere  with  s^  at  the  center.  (An 
equal  area  projection  Is  used,  see  Appendix  8.5.)  The  nulls  and  maxima  for 
the  gamut  of  possible  symmetric  targets  are  plotted  on  polarization  charts  in 
Figure  4.  In  these  charts,  H  Is  at  the  top,  V  at  the  bottom,  tC  In  the 
center,  and  the  perimeter  Is  linear  polarization.  The  central  meridian  Is 
2v  *  0*,  while  the  perimeter  Is  2v  »  *  90*.  Only  the  left-sense  half  of  the 
sphere  is  represented,  so  that  usually  only  one  of  the  two  co-pol  nulls  or  the 
two  x-pol  maxima  Is  plotted;  the  others  of  the  pair  would  appear  symmetrically 
on  the  right-sense  chart.  Going  across  the  page,  the  charts  are  for  targets 
with  the  same  ?  (or  1  *Vy / j)  but  with  different  values  of  v  (or 
relative  phase  (aHH/ ayy);  while  going  down  the  page,  the  charts  are  for 
targets  with  different  y  and  tha  same  v  .  Note  that  Is  always  at  the 
top  (it  could  also  be  at  the  bottom)  and  that  all  the  optimal  polarizations 
lie  on  the  projection  of  a  great  circle  (a  meridian),  except  for  the 


...2V*  ANY  VALUE 
0 


fcf-N 


FIG.  4  Polarization  nulls  and  maxima  for  plana  symmetric  targets  for  various 
values  of  parameters  v,  y .  One  hemisphere  shown.  Note:  shown  for 


Uhh  I  >  l>vv  I  For  l*vv  '  •>  I*hh  *  turn  disc  upside  down. 


degenerate  cases.  In  addition  to  its  polarization  chart  description,  the 
value  m  is  required  to  characterize  a  target.  The  parameters  m,  y,  v,  are 
aspect  and  frequency  dependent,  so  that  a  complete  characterization  would 
consist  of  null  or  max  loci  as  functions  of  aspect  and  frequency  together 
with  m  as  a  function  of  aspect  and  frequency. 

The  position  of  the  co-pol  null  on  the  polarization  chart  can  also  be 
interpreted  directly  as  representing  the  parameters  v,  y  of  the 
scattering  matrix.  The  longitude  corresponds  to  2v,  the  latitude  corresponds 
to  (2y  -  90’). 

From  Figure  4  we  can  observe  several  important  properties  of  the 
optimal  polarization  locations  pertinent  to  target  classification:  (1)  As 
noted  above,  the  co-pol  max  and  saddle  (and  the  x-pol  nulls)  are  always  at  H 
or  V;  (2)  the  x-pol  maxima  always  occur  on  the  equator,  which  is  45* 
elliptical  polarization  or  a  linear  real  combination  of  RC  and  LC;  (3)  the 
co-pol  nulls  may  appear  anywhere.  We  may  think  of  the  locations  ©  >© . 
and  ©  as  having  respectively  0,  1,  and  2  degrees  of  freedom.  Hence  we  can 
already  anticipate  that  for  two  synroetric  targets,  the  co-pol  null  location 
will  be  the  best  discriminant  (after,  of  course,  the  polarization  reference 
has  been  established). 

There  are  two  ways  in  which  target  asymmetry  will  modify  Figure  4. 
First  is  the  transformation  Qy  ,  which  is  the  same  as  or  is  equivalent  to 
rotating  the  target  by  f  about  the  RLOS,  and  results  in  a  rotation  of  Figure 
(4)  by  2*  about  the  center  point.  This  asymmetry  would  be  removed  in 
practice  by  establishing  the  polarization  reference  plane  to  be  parallel  to 
the  symmetry  plane,  hence  restoring  Figure  4.  That  is,  polarimetry  can  be  and 
should  be  made  Invariant  under  target  rotation  about  the  RLOS.  The  second  way 
in  which  a  target  may  be  asymmetric  corresponds  to  a  rotation  of  Figure  4 
about  $2  tbe  angle  2  t  ) ;  that  is,  the  polarization  chart  rotates  out 
of  the  plane  of  the  paper.  Depending  on  the  degree  of  this  asymmetry,  it  mqy 
turn  out  to  be  possible  to  discriminate  between  a  symmetric  and  a 
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non-symmetric  target  (or  between  two  non-symmetric  targets)  on  the  basis  of 
this  rotation.  Whether  or  not  this  Is  true  will  depend  on  the  mechanism  which 
underlies  the  structure  of  the  scattering  matrix,  namely,  the  spatial  and 
amplitude  relations  between  the  scattering  centers  that  make  up  the  scatterer 
and  the  dynamic  (with  frequency  and  aspect)  behavior  of  these  relations.  This 
will  be  Investigated  In  a  later  section.  It  may  turn  out  that  the  asymmetry 
character  Is  not  a  good  discriminant,  In  which  case  It  can  be  removed  by 
choice  of  (elliptical)  polarisation  basis,  thus  again  obtaining  Figure  4.  In 
any  case,  the  general  behavior  of  Figure  4  can  be  considered  fundamental,  as 
It  can  be  obtained  for  any  target  by  a  suitable  choice  of  the  transmit  and 
receive  polarizations. 


2.6  SIMPLE  APPROXIMATE  DESCRIPTIONS  OF  SCATTERING  CHARACTERISTICS: 

SINGLE-BOUNCE,  DOUBLE-BOUNCE,  DIPOLE. 

We  can  Identify  some  of  the  Charts  in  Figure  4  as  belonging  to  certain 
types  of  scattering  centers.  Chart  (a)  represents  the  single-bounce 
circularly  symmetric  specular  scatterer  with  A  «  I  (2v  »  0,  Zy  •  90*),  while 
(c)  Is  the  double-bounce  scatterer  (e.g.  a  dihedral).  At  low  frequency,  the 
specular  return  from  a  doubly  curved  surface  (e.g.  prolate  spheroid)  might  be 
represented  by  (d).  A  horizontal  edge  or  dipole  is  represented  by  (j).  A 
vertical  edge  Is  represented  by  (v1)  turned  upside-down. 

The  symmetric  target  rotated  about  the  RLOS  Is  represented  by  the  same 
figures  rotated  by  2H*  about  the  center  point  (axis  s^)* 

When  discussing  polarization  characteristics  of  a  target,  it  Is 
convenient  to  make  use  of  approximate  descriptions.  We  collect  these 
descriptions  in  this  separate  subsection  for  the  sake  of  reference.  We 
restrict  the  concepts  to  symmetric  targets  aligned  with  the  H-V  polarization 
basis. 
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By  single-bounce  we  will  mean  a  scatterer  with  a  scattering  matrix 
that  is  approximately  the  unit  matrix,  with  no  phase  difference  between  the  HH 
and  VV  terms.  The  specular  return  from  a  sphere  or  plate  is  of  this  type. 

That  is,  2v  ■  0  and  Zy  -  90*  for  a  “single-bounce"  scatterer.  The  co-pol 
nulls  occur  at  LC  and  RC;  that  is,  2v  -  0,  180*  and  Zy  -90*;  or  in  terms  of 
the  polarization  ellipse,  at  t  ■  *45*,  with  ¥  undefined.  The  nulls  are  thus 
located  at  the  center  of  the  polarization  chart,  as  In  Figure  4a.  When 
discussing  target  type,  we  may  Identify  a  "single-bounce"  scatterer  as  one  for 
which  these  conditions  hold  only  approximately. 

By  double-bounce  we  will  mean  a  scatterer  having  a  matrix  with 
approximately  unit  diagonal  entries  of  opposite  sign.  An  inside  corner 
(dihedral)  is  of  this  type.  That  is,  2v  -  *90*,  2?  -  90*  for  a 
"double-bounce"  scatterer.  The  co-pol  nulls  occur  at  45*-lineer  polarization, 
that  is  2v  m  *90*  and  2-r  «  30*  ;  or,  In  terms  of  the  polarization  ellipse 
at  t  »  0,  H1 »  *45*.  The  nulls  are  located  at  the  equator  end  points  of 
the  polarization  chart,  as  in  Figure  4c. 

By  dipole  we  will  mean  a  scatterer  with  only  one  non-zero  (HH  or  VV) 
entry  in  the  scattering  matrix.  A  horizontal  dipole  has  an  HH  entry.  Its 
co-pol  nulls  coalesce  at  VV,  that  is  2-r  «  0  with  2v  undefined;  or  in  terms 
of  the  polarization  ellipse,  at  t  •  0,  f  ■  0  .  The  null  Is  located  at 
the  bottom  of  the  polarization  chart,  as  in  Figure  4g. 

A  vertical  dipole  has  a  scattering  matrix  with  only  a  non-zero  VV 
entry.  Its  co-pol  null  occurs  at  HH,  and  thus  has  the  polarization  chart 
representation  of  Figure  4g  turned  upside-down,  with  the  co-pol  null  at 
x  -  0,  ¥  -  90*.  Although  the  canonical  form  description  is  defined  for 

0  <  2Y  <  90*  ,  It  is  convenient  to  think  of  the  vertical  dipole  co-pol  null  as 
occurring  at  Zy  -  180*  .  That  Is,  for  the  sake  of  the  polarization  chart 
description,  2-r  will  always  be  measured  with  respect  to  the  south  pole. 
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2.7  DYNAMIC  BEHAVIOR  OF  NULLS  AND  MAXIMA  FOR  OUMBELL  TARGETS 


We  now  begin  our  investigation  of  the  dynamic  properties  of  the 
optimal  polarizations;  that  is,  the  null  and  maxima  locations  as  functions  of 
aspect  and/or  frequency.  First  note  that  the  scattering  matrix  of  a  single 
scattering  center  is  only  a  slow  function  of  aspect  or  frequency.  For 
compound  targets  the  dominant  effect  on  the  scattering  matrix  is  due  to 
changes  in  the  relative  phase  between  scattering  centers.  Accordingly  we  will 
begin  our  study  with  a  simple  dumbell  target  composed  of  two  independent 
scattering  centers,  each  of  which  is  of  the  general  type  discussed  in  Section 
2.5  and  each  of  which  can  be  considered  to  be  approximately  constant  over  a 
certain  region  of  aspect  or  frequency  variation.  The  modeling  of  compound 
targets  as  the  phase-dependent  sum  of  independent  scattering  centers  is 
discussed,  for  example,  in  Reference  [8].  It  is  valid  when  scattering  centers 
are  separated  by  more  than  a  few  wavelengths  (as  is  the  case  for  the  targets 
considered  in  this  report).  The  variable  is  the  relative  phase  between  the 
two  scattering  centers  which  can  be  considereu  to  be  due  to  aspect  variation 
as  the  dumbell  rotates  in  the  plane  of  symmetry  or  (equivalently)  due  to 
frequency  variation.  For  aspect,  o  ,  variation  the  relative  phase  is 

=  (4irL/ x)cos(«q  +  Aa) ,  where  L  is  the  separation  between  scatterers  . 

For  frequency  variation,  we  can  write  *(4it/Xq)(1  +  Af/fgHL  cos  a^) . 

We  will  observe  null  and  max  motion  as  functions  of  varying  0^  .  The 
horizontal  and  vertical  scattering  coefficients  of  scatterer  1  are  given  by 

j®ui  j®Vi  ^eH2  ^°V2 

h^e  ni  and  v^e  *  and  for  scatterer  l  by  h^e  and  v^e 


Combining  these  and  manipulating  to  remove  the  common  constant  phase,  the 
scattering  matrix  can  be  written 


a  Pj012/2  +  A  e’^12 

A,e  +  A9e 


/  2 


HH 


VV' 


with 


HH 


Ml  2  -Ml  2 

h^e  +  h2e 


j(0  +  aovh)/2  -j(0  +  AoyH)/2 

aVV  =  vle  +  v28 


(31) 


where 


aovh  =  evi  '  9V2  "  °H1  +  °H2 
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This  Is  the  most  general  dumbell  target  composed  of  constant  Individual 
scatters,  each  with  Its  symmetry  plane  coincident  with  the  symmetry  plane  of 
the  rotating  dumbell.  (A  completely  general  dumbell  would  be  asymmetric,  with 
the  Individual  scatterers  not  aligned  with  the  plane  of  rotation.) 

The  scattering  matrix  and  null  and  maxima  locations  were  computed  for 
various  combinations  of  the  parameters  h^,  h^,  v^,  v2,  dyH  as  d 
varied  from  -*  to  +*  ,  i.e.  as  the  RCS  goes  through  one  scintillation 
period.  One  case  Is  Illustrated  In  detail  In  Figure  5.  This  is  the  equal 
area  polarization  chart  described  in  Appendix  8.5.  The  center  Is  LC 
polarization,  H  is  at  the  north  pole  and  V  at  the  south  pole.  The 
egg-shaped  curved  Is  the  (left-sense)  co-pol  null  locus  as  the  relative  phase 
d  between  scatterers  goes  through  360*.  The  symmetric  dashed  curve  is  the 
locus  of  the  other  (right-sense)  co-pol  null,  which  appears  on  the  back  of  the 
Poincare  sphere.  The  Huynen  fork  is  indicated  for  a  particular  value  of  0  , 
Identified  by  the  points  ^),  and  (J)  .  These  points  lie  on  a 

meridian.  As  d  varies,  the  cross-pol  max,  ©  ,  moves  along  the  equator. 

The  co-pol  max,  ©  ,  Is  always  at  H  (except  when  the  co-pol  nulls  move  to 
the  upper  hemisphere,  in  which  case  is  at  V). 

As  discussed  in  Section  2.5,  the  longitude  of  the  co-pol  null  position 
corresponds  to  the  parameter  2v  of  the  scattering  matrix,  while  the  latitude 
corresponds  to  (2-r  -  90*).  (It  Is  convenient  to  think  of  (2y)  as  taking 
values  from  0  to  180*  ,  even  though  it  is  defined  for  0  to  90*  -  values 
of  2y  >  90*  mean  that  the  co-pol  max  i  s  at  V  . ) 

A  number  of  typical  dumbell  cases  are  illustrated  in  Figure  6.  Only 
the  left-sense  co-pol  null  loci  are  shown.  (The  point  ©  and  the  ^^locus 
can  be  found  using  the  construction  of  Figure  5.)  Some  general  observations 
can  be  made:  When  one  of  the  scattering  centers  is  dominant  (Figures  6a  and 
6b),  the  co-pol  null  makes  only  small  excursions.  The  relative  phase,  oVH, 
between  the  horizontal  and  vertical  scattering  constants  affects  the  average 
value  of  2v  (Figure  6b).  The  case  illustrated  in  Figure  6a  is  for 
hj  »  vi  ,  which  is  the  cnse  if  scatterer  1  is  a  specular  point.  The 


dumbell  of  Figures  6c  end  6d  has  the  same  scattering  strengths,  but  now  1  is 
predominantly  horizontal  and  2  predominantly  vertical.  The  null  locus 
excursion  is  small  in  y  but  large  in  v  .  The  corresponding  variation  of  m 
(maximum  co-pol  voltage)  is  illustrated  in  Figures  7a,  b.  c,  d.  Because  one 
scatterer  is  dominant,  the  fluctuation  in  m  is  small. 

Figures  6e  and  6f  show  what  happens  if  both  scatterers  are  about  equal 
vertically  but  scatterer  1  is  dominant  horizontally.  This  case  could 
correspond  to  a  specular  scatterer  1,  and  some  kind  of  vertical  edge-like 
structure  for  scatterer  2.  The  null  trajectory  in  this  case  is  much  larger, 
showing  variation  in  both  the  v  and  y  parameters.  The  corresponding  plots 
of  m  in  Figure  7e  and  7f  show  larger  fluctuations  in  m.  Finally,  the  last 
two  cases  (Figures  6g,  6h  and  7g,  7h)  show  the  case  of  two  nearly  equal 
interfering  scatterers.  Constructive  and  destructive  interference  effects  are 
very  large:  For  Figure  6g,  the  horizontal  and  vertical  oscillations  are  in 
phase j  resulting  in  no  movement  of  the  null  locus;  but  for  Figure  6h,  the 
oscillations  are  out  of  phase,  resulting  in  very  large  excursions  of  the 
co-pol  null. 

It  is  important  to  note  that  a  symmetric  target  is  completely 
described  by  a  co-pol  null  locus  and  the  corresponding  m-variaticn. 

A  significant  result  of  this  simulation  is  that  even  for  simple 
compound  targets  the  null  trajectory  can  appear  just  about  anywhere  on  the 
polarization  chart  ind  may  exhibit  large  variations  in  both  parameters  2v 
and  Zy  .  Large  xcursions  in  the  null  locus  occur  in  just  the  same  way  as 
large  fluctuations  in  the  RCS  when  observed  as  a  function  of  aspect  or 
frequency  -  namely,  as  a  result  of  Interference  effects.  The  relationship 
between  RCS  fluctuations  and  null  position  fluctuations  is  complicated  because 
it  involves  the  relative  phase  of  the  HH  and  VV  components  of  the 
scattering  matrix  as  well  as  the  "relative  phase"  between  the  H  and  V 
scintillations.  However,  it  can  be  said  that  if  RCS  fluctuations  are  rapid, 
then  null  location  fluctuations  are  equally  rapid. 
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It  Is  possible  to  anticipate  from  Figure  6  the  null  locus  behavior  for 
more  general  2-scatterer  targets:  Suppose  we  observe  the  target  through 
several  scintillation  periods,  and  that  the  magnitudes  of  the  scattering 
centers  change  slowly  as  a  function  of  aspect  (or  frequency),  as  would  be  the 
case  for  a  typical  radar  target.  We  then  expect  the  shape  and  centroid  of  the 
loci  of  Figure  5  to  change  as  the  null  moves  around  on  the  locus,  resulting  In 
a  spiral-like  or  smeared  out  locus.  This  will  be  seen  to  be  the  case  for  the 
cylinder  (Section  2.9)  which  over  some  range  of  aspect  behaves  like  two  or 
three  scattering  centers  with  slowly  varying  strengths. 


2.8  DEPOLARIZATION  DUE  TO  AVERAGING  OVER  A  SCINTILLATION  PERIOD 

When  the  separation  between  two  scattering  centers  is  many 
wavelengths,  the  fluctuations  In  RCS  and  In  null  locations  Is  rapid  for  even 
small  changes  In  aspect  (or  frequency).  In  that  case  It  may  be  impractical  or 
undesirable  to  follow  these  fluctuations  in  a  measurement;  Instead,  an  average 
measurement  is  taken.  The  result  of  averaging  Is  that  the  received  power  Is 
depolarized.  This  Is  equivalent  to  the  statement  that  for  the  received  Stokes 
vector  defined  In  Equation  (4)  we  have 

s02  *  sx2  +  s22  ♦  s32  .  (32) 

A  result  Is  that  It  is  no  longer  possible  to  observe  a  co-pol  null,  but  rather 
a  co-pol  minimum.  A  complete  treatment  of  the  depolarized  case  is  beyond  the 
scope  of  this  report.  It  suffices  for  our  present  purposes  to  note  that  the 
location  of  the  co-pol  minimum  (when  the  fluctuations  are  not  too  large)  is 
approximately  given  by  the  average  of  the  instantaneous  null  locations.  The 
co-pol  minimum  location  is  Indicated  by  the  point  A  in  each  of  Figures  6  for 
the  average  over  one  cycle  of  a  fluctuation  for  the  dumbell.  The  points  A 
were  computed  from  the  average  depolarized  Stokes  vector  In  just  the  same 
manner  as  the  instantaneous  null  locations  were  computed  from  the  fully 
polarized  instantaneous  Stokes  vectors. 
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Data  averaging  so  as  to  remove  the  large  fluctuations  In  the 
polarization  characteristics  Is  pursued  to  some  extent  In  Section  3  with  the 
aircraft  data. 

For  purposes  of  target  classification,  the  fluctuations  may  serve  as  a 
discriminant,  since  they  are  related  to  target  size  and  structure.  However, 
the  fluctuations  are  exhibited  first  and  foremost  by  the  received  power.  It 
is  not  clear  whether  null  location  fluctuation  adds  any  additional  information 
about  the  target.  It  is  therefore  likely  that  the  average  null  location  is 
the  more  useful  quantity  to  be  derived  from  polar imetry.  These  comments  apply 
to  the  present  case  of  unresolved  scattering  centers.  When  scattering  centers 
are  resolved,  so  that  they  can  be  characterized  individually,  then  their 
polarization  properties  (discussed  in  Section  2.5)  would,  of  course,  be 
readily  apparent. 


2.9  CO-POL  NULL  LOCUS  FOR  RIGHT-CIRCULAR  CYLINDER 

The  simple  aircraft  models  which  are  the  ultimate  subject  of 
investigation  in  this  report,  all  consist  of  right-circular  cylinders  (RCC) 
with  fins  attached.  Accordingly,  we  now  investigate  the  polarization 
chai'act"  sties  of  the  RCC  by  itself. 

The  backscatter  modeling  of  the  RCC  used  in  this  report  is  given  for 
reference  in  Appendix  8.6.  The  cylinder  is  a  symmetric  target,  and  we  assume 
in  this  discussion  that  the  H-V  polarization  basis  has  been  aligned  with  the 
plane  of  symmetry  (E^  in  plane  of  symmetry,  £y  perpendicular  to  plane 
of  symmetry).  The  aspect  change  Is  a  rotation  In  the  plane  of  symmetry.  The 
HH  and  VV  r*rf'r  cross  sections  for  a  cylinder  with  length  to  diameter 
r&'i .  i,0  *  -  ka  -  2.4  (where  a  is  the  radius  and  k  -  2* lx)  are 

given  in  Figure  8,  as  functions  of  aspect.  Not  shown  is  the  relative  phase 
between  the  HH  and  VV  responses,  although  of  course,  this  is  an  Important 
contributor  to  t’.~  nolarlzatlon  characteristics. 
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Netr  nose-on  (a  >0),  the  specular  response  Is  that  of  a  circular 
disc,  which  Is  polarization  Independent.  Near  broadside  (a  -  90*)  the  return 
Is  again  specular  and  nearly  polarization  Independent.  Since  this  case  Is  In 
the  resonance  region,  there  Is  a  slight  difference  In  the  HH  and  VV 
responses  at  broadside,  due  both  to  the  polarization  dependent  specular  and 
creeping  wave  contributions  to  the  response.  At  Intermediate  angles, 
scattering  may  be  considered  to  occur  at  the  three  visible  edges  of  the 
cylinder.  The  edge  returns  are  polarization  dependent;  Individually  they  vary 
slowly  with  aspect,  but  due  to  the  varying  relative  phase,  the  total  response 
fluctuates.  The  HH  and  VV  response  fluctuations  are  different  In  both 

phase  and  amplitude.  As  a  consequence,  the  co-pol  null  location  will  change 

with  aspect. 

The  scattering  matrix  behavior  as  a  function  of  aspect  (0  to  90*)  Is 
illustrated  in  Figure  9  In  terms  of  the  canonical  parameters  m,  2v,  and  2?. 
These  three  curves  completely  describe  the  scattering  properties  of  the 

cylinder.  The  polarization  maximum  m  (max  received  voltage).  Figure  9a, 

% 

traces  the  maximum  of  the  HH  and  VV  response  curves  of  Figure  8.  Because 
the  target  Is  symmetric,  m  occurs  at  either  H  or  V  on  the  polarization 
chart  (i.e.  at  tne  poles).  Figure  9b  Is  2v  ,  which  Is  the  relative  phase 
between  the  HH  and  VV  returns.  Note  that  2v  -  0  from  <*  ■  0  to  a  -  35*, 

where  the  response  Is  polarization  Independent.  At  a  -  90*,  2v  Is  again 

nearly  zero,  though  not  quite  because  of  the  slight  polarization  dependence  of 
the  broadside  specular.  At  Intermediate  aspects,  2v  is  a  smoothly  varying 
function,  going  through  several  360*  cycles.  Also  note  that  there  are  two 
solutions  for  2v  ,  so  that  another  curve  could  be  drawn,  displaced  by  180*. 
Figure  9c  shows  that  2*  vs.  aspect  makes  only  small  deviations  from  2y  -  90*. 
Phenomenologically,  the  cylinder  behaves  as  a  single-bounce  scatterer  (2v  -  0, 
2y  •  90*)  near  nose-on  and  broadside;  while  at  Intermediate  angles  It 
oscillates  between  approximate  single-bounce  and  double-bounce  (2v  «  *90*. 

Zy  •  90*)  behavior. 
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The  cross-pol  maximum  polarization,  as  discussed  In  Section  2.5,  Is 
defined  by  the  parameter  2v  (the  longitude  on  the  polarization  chart)  and 
occurs  on  the  equator  of  the  polarization  chart.  The  co-pol  null  polarization 
Is  determined  by  both  2v  and  2y  •  (Since  the  co-pol  null  locus  has  one 
more  degree  of  freedom  than  the  cross-pol  max  locus,  the  cross-pol  max  locus 
will  usually  not  be  discussed  -  It  Is  simply  the  projection  onto  the  equator 
along  longitude  lines  of  the  co-pol  null  locus.) 

A  partial  co-pol  null  locus  Is  plotted  on  the  polarization  chart  for 
aspects  between  35*  and  52.5*  (Figure  10),  The  longitude  of  the  co-pol  null 
is  2v  and  the  co-latitude  Is  2y.  Below  a  -  35*  ,  the  co-pol  null  Is  at 
the  center  (LC)  and  above  50*;  it  makes  severzl  traversals  around  the  Poincare 
sphere,  near  the  equator  (as  indicated  by  Figures  9b  and  9c)  before  coming  to 
rest  near  LC  (or  RC)  again  for  a  -  90*  . 

Several  comments  are  In  order  about  these  null  locus  plots  to  avoid 
confusion.  First,  recall  that  there  are  actually  two  loci,  displaced  by 
2v  •  180*.  We  always  plot  only  one  locus,  the  other  would  appear 
symmetrically  about  the  H-V  axis  on  the  polarization  chart.  Second,  the 
choice  of  one  null  or  the  other  In  tracing  the  locus  Is  completely  arbitrary. 
The  computer  algorithm  used  to  generate  these  curves  chooses  the  null  closest 
to  the  previously  computed  null,  thus  tracing  out  a  smooth  curve;  If  a 
different  choice  were  made,  or  if  the  aspect  steps  between  samples  were  too 
coarse,  the  locus  would  exhibit  discontinuous  skips  from  one  side  of  the 
polarization  chart  to  the  other.  There  Is  no  significance  to  this  skipping; 
it  is  not  incorrect,  but  makes  a  confusing  picture.  Third,  the  polarization 
chart  shows  only  one  hemisphere  of  the  Poincare  sphere;  here,  nominally  the 
left-sense  hemisphere.  In  presenting  tha  null-locus,  we  have  here  superposed 
the  left-  and  right-sense  hemispheres,  in  order  for  the  locus  to  show  up  as  a 
smooth  curve.  This  means  that  In  Figure  10,  as  the  locus  moves  from  the  LC 
center  to  the  perimeter,  it  continues  on  the  back-side  of  the  Poincare 
sphere.  In  later  figures,  where  the  full  locus  Is  shown,  the  apparent 
figure-eight  behavior  Is  actually  an  undulating  curve  continuously  turning 
around  the  sphere. 
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For  a  longer  cylinder,  the  response  fluctuations  are  more  rapid.  The 
HH  and  VV  radar  cross  sections  for  a  L/D  »  10  cylinder  at  ka  »  2.4  are 
shown  in  Figure  lib;  the  responses  of  the  L/D  *  5  cylinder  are  duplicated  in 
Figure  11  a  for  easy  comparison.  The  complete  characterizations  of  the 
scattering  matrices  for  these  two  cases  in  terms  of  m,  2v,  and  2y  are  shown 
in  Figure  12.  The  left  column  (a,b,c)  duplicates  previous  results  for  the 
L/D  ■  5  cylinder,  the  right  column  (d,e,f)  are  the  corresponding  parameters 
for  the  L/D  *  10  cylinder.  The  behavior  of  the  two  cylinders  is  very 
similar,  the  main  difference  being  the  more  rapid  variation  of  the  larger 
cylinder.  The  (2v,  2y)  loci,  that  is  the  co-pol  null  loci,  are  plotted  on  the 
polarization  charts  in  Figure  13.  The  plots  have  been  broken  up  into  four 
intervals  of  aspect  a  ,  to  display  the  loci  more  clearly.  As  before,  the 
left  column  is  the  L/D  *  5  cylinder,  the  right  column  is  the  L/D  »  10 
cylinder.  The  tick  marks  on  the  curves  correspond  to  2-1/2*  aspect  change. 

The  locus  moves  around  the  Poincare  sphere  in  an  undulating  fashion,  except 
that  for  a  •  0  to  40*  the  null  is  fixed  at  the  center  (circular 
polarization).  We  see  that  the  cylinder  may  act  at  some  aspects  like  a 
single-bounce  scatterer  (null  near  CP),  and  at  other  aspects  like  a 
double-bounce  scatterer  (null  near  45*  linear  polarization,  i.e.  the  extremes 
of  the  equator).  One  traversal  around  the  Poincare  sphere  looks  similar  to 
some  of  the  dumbell  results  of  Figure  6,  as  was  expected,  since  the  cylinder 
consists  of  three  scatters;  one  strong,  one  weaker,  and  the  third  much 

weaker.  The  null  loci  are  generally  in  the  upper  (H)  hemisphere,  reflecting 

the  facts  that  the  VV  response  of  the  cylinder  is  somewhat  larger  than  the 
HH  response  most  of  the  time  and  that  the  HH  response  minima  are  deeper. 

At  broadwide,  the  VV  response  Is  smaller  than  the  HH  response,  so  that  the 
final  null  position  is  in  the  lower  hemisphere  (not  clear  in  the  figures). 

For  a  very  thin  cylinder  (ka  small,  L/D  large),  we  expect  behavior  like  a 
dipole  (null  at  V)  at  broadside. 

The  aircraft  models  to  be  studied  in  this  report  scale  in  the  UHF 
region  to  higher  values  of  ka  ,  namely  ka  ■  5.24  and  ka  -  8.15 

with  L/D  -  10  .  The  behavior  of  cylinders  at  those  scalings  are  shown  in 

Figures  14  and  15.  The  organization  of  these  plots  is  similar  to  that  of  the 
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FIG.  14  Parameters  m,  2v  and  2y  vs.  aspect  for  RCC,  L/D  IQ  at  ka  -  5.24  (left  column) 
end  ka  -  8.15  (right  column). 


previous  two  figures;  here  the  left  column  is  for  ka  «  5.24  and  the  right 
column  for  ka  -  8.15  ,  with  L/D  ■  10  in  both.  The  results  are  similar  to 
the  cases  with  lower  ka  ,  except  that  the  variation  with  aspect  is  more 
rapid.  An  interesting  detail  of  the  locus,  which  (by  coincidence!)  didn't 
show  up  in  the  lower  ka  cases  Is  a  spiral  motion.  .Part  of  the  Figure  15b 
locus  is  reproduced  in  Figure  16  to  show  this  behavior  more  clearly.  This 
effect  had  been  predicted  in  the  dumbell  study  earlier  in  this  section.  It  is 
caused  by  one  dominant  scatterer  interfering  with  a  weaker  scatterer,  through 
several  interference  cycles,  while  the  relative  strengths  of  the  scatterers 
change  slowly. 

Note  that  for  scaling  of  the  aircraft  models  to  S-band  or  X-band 
(about  a  factor  of  10  higher  in  ka  )  the  fluctuations  of  the  null 
characteristics  would  be  a  factor  of  10  or  so  more  rapid  yet.  It  is 
clear  that  any  target  discrimination  scheme  based  upon  null  characteristics  of 
an  unresolved  target  could,  at  these  frequencies,  not  be  used  when  using  a 
point  by  point  comparison.  Rather,  average  or  statistical  features  would  have 
to  be  used.  For  example:  For  a  cylinder,  Zy  is  on  the  average  slightly 
greater  than  90*,  and  doesn't  go  much  below  90*.  The  presence  of  other  large 
scatterers  on  the  structure  are  likely  to  alter  this  picture.  It  is  therefore 
pertinent  to  examine  more  complex  targets,  as  is  done  in  Section  3. 
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2.10  EXTENSION  OF  PREVIOUS  RESULTS  TO  A  NON-DIAGONAL  SCATTERING  MATRIX 


So  fa*'  Me  have  discussed  symmetric  targets  which  have  a  diagonal 
scattering  matrix,  and  which  can  be  parameterized  by  the  three  real 
quantities  m,  2v,  and  Zy.  The  most  general  monostatic.,  relative  phase 
scattering  matrix  Is  symmetrical  and  requires  two  more  real  parameters 
(corresponding  to  the  complex  off-diagonal  term).  From  the  discussion  In 
Appendices  2  and  4  and  Section  2.4,  the  general  scattering  matrix  A,  Is 
related  to  a  diagonal  matrix  D,  by  the  similarity  transformation 


where 


(Q*Qt  )  A(QyQT  ) 


(33) 

(34) 


(35) 


The  effect  on  the  Stokes  vector  representing  a  point  on  the  Poincare  sphere  is 

a  rotation  of  |2t.  about  the  S2  (45°  linear  polarization)  axis, 

followed  by  a  rotation  of  2  f  about  the  (circular  polarization)  axis. 

In  particular,  the  polarization  maximum,  m,  will  be  rotated  by  2x  and  2  ¥ 
from  its  position  at  H  or  V  for  the  diagonal  matrix.  Now,  a  rotation  of 
2  V  of  the  Poincare  sphere  corresponds  physically  to  a  rotation  by  V  of  the 
target  or  of  the  polarization  reference  plane)  about  the  RLOS.  A  rotation  of 
2t  corresponds  to  a  change  of  x  in  the  ellipticity  of  the  polarization 
ellipse  (Figure  1).  Considering  our  polarization  chart  representation  of  the 
null  locus  (e.g.  Figure  5),  the  chart  would  be  rotated  out  of  the  plane  of  the 
page  and  then  rotated  in  the  page  about  the  center.  This  operation  destroys 
some  of  the  convenience  of  the  polarization  chart.  For  instance,  the  two 
co-pol  nulls  will  now  no  longer  be  located  symmetrically  about  the  HV  axis. 
Therefore,  in  our  work  with  non- symmetric  targets,  we  will  usually  remove 
the  t  and  ¥  transformations  from  the  matrix  and  continue  to  plot  the 
polarization  characteristics  in  the  standardized  way  with^M)  at  the  top  (or 
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bottom).  The  complete  description  of  the  scatterer  will  then  require  the  x 
and  4'  history  in  addition  to  that  of  m  and  the  co-pol  null  locus. 
Alternatively  we  can  present  five  separate  plots  of  m,  2\>  ,  2*  *  2x  , 
and  2  4'. 


The  solution  for  4*  and  x  may  be  found  by  writing  out  Equation  (33) 
In  full .  The  result  is 

2Re((aHH  +  ®w)  a8v) 

tan  24*  »  2  2  (36) 

laHH  |  -  laW  | 


2Im((aHH  -  avv)anv) 

tan  2x  =  cos  24'  2  2  (37) 

laHH  I  -  laVV  I 


The  diagonalized  matrix  is  then  given  by 

On  =  c*2  ann  +  s*2  ayy  +  2c*s*  a^v  , 

022  *  c2  ayv  +  s2  aHH  -  2cs  aHV »  (38) 


where 

c  ■  cosx  cos4'  +  j  sinx  sinf  , 
s  -  cosx  sin'y  -  j  sinx  cos'?  . 

The  co-pol  max  polarization  before  the  transformation  is  given  by  (-4*  ,  -x). 
After  the  transformation,  the  co-pol  max  is  either  at  H  or  V  ,  as  for  the 
diagonal  scattering  matrices  discussed  before.  24' and  2t  are  both 
restricted  to  between  ±90*. 
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2.11  EFFECT  OF  POLARIZATION  CALIBRATION  ERRORS  ON  NULL  MEASUREMENTS 


As  described  in  Section  4,  great  care  was  taken  in  the  measurement 
program  to  assure  proper  calibration  of  the  two  channels  of  the  dual -polarized 
radar.  The  calibration  procedure  Is  rather  complicated  due  to  the  fact  that  a 
complex  2x2  calibration  matrix  has  tc  be  determined,  as  opposed  to  the 
single  real  constant  required  for  the  single-polarization  case.  It  is  of 
interest  to  know  just  how  well  such  a  calibration  needs  to  be  carried  cut. 

First,  consider  a  calibration  error  of  the  form 


that  is,  the  vertical  channel  introduces  additional  attenuation  and  a  phase 
change.  The  measured  scattering  matrix  will  he  of  the  form  =  EA, 
where  A  is  the  actual  scattering  matrix.  The  constant  k  is  just  the 
absolute  calibration  factor,  the  same  as  -equlred  for  single-polarization 
measurements,  and  will  not  be  discussed  further. 

The  effect  on  t.he  apparent  polarization  characteristics  is  seen 
readily:  c  changes  the  relative  magnitude  of  the  a^  and  avv 
components  of  the  matrix  and  hence  changes  the  value  of  2y  .  That  Is,  the 
measured  value  y1  is  related  to  the  true  value  y  by 

2  2 
tan  y'  m  t  tan  y  . 

This  causes  a  distortion  in  the  latitude  of  the  polarization  chart. 

The  effect  of  Is  to  change  the  apparent  phase  between  the 
aHH  and  avv  components  and  hence  change  the  value  of  2v  .  That  Is, 
the  polarization  chart  will  be  rotated  about  the  HV  axis.  In  particular,  a 
point  at  LC  (chart  center)  will  move  off-center,  hence  a  single-bounce 
scatterer  will  appear  as  not  quite  single-bounce.  To  get  some  feeling  for 
these  errors,  the  dumbell  cases  of  Figure  6  have  been  recalculated  for  an 


error  of  10  percent  In  relative 


amplitude  and  30°  in 


phase  (c 


0.9,  *e 


30* ). 


Th-i  results  are  shown  in  Figure  17.  We  see  that  the  effect  of  calibration 
errors  are  not  serious  unless  a  precise  knowledge  of  2v  or  2y  is  desired. 


If  the  error  matrix  also  has  off-diagonal  terms,  the  effect  is  more 
complicated.  If  the  off-diagonal  errors  are  symmetric  (which  is  unlikely), 
then  the  error  contains  additional  rotations  in  t  and  ¥  ,  as  discussed  In 
Section  2.10.  If  the  errors  are  non-symmetric,  the  effect  depends  on  just  how 
the  data  is  manipulated  to  obtain  the  polarization  characteristics.  Clearly, 
it  is  best  to  obtain  both  a^v  and  ayH  by  independent  measurement,  and 
then  average  them  to  remove  the  non-symmetric  part  of  the  error  matrix. 


FIG.  17  Copol  loci,  ah  owing  effect  of  error  in  VV  channel  of  10%  magnitude,  3C°  phase. 
Compere  this  with  figure  6. 


i 


\ 

■j 


J 


I 


-63“ 


SECTION  3 


COMPUTED  POLARIZATION  PROPERTIES  OF  SIMPLE  AIRCRAFT  MODELS 


This  investigation  into  the  polarization  characteristics  of  al-craft 
models  followed  two  paths  simultaneously  and  interactively:  the  direct 
measurement  of  polarization  responses  and  the  computation  of  the  responses, 
1,e.,  target  modeling.  The  reasons  for  pursuing  the  target  modeling  approach 
are  compelling;  through  it,  the  dependence  of  polarization  characteristics 
upon  target  parameters  (such  as  the  shapes,  sizes,  and  interrelationships  of 
its  parts)  and  observation  parameters  can  be  studied  readily  without 
restriction  and  without  the  nuisance  of  noise  and  other  measurement  hazards. 
Phenomena  yield  more  readily  to  understanding  so  that  predictions  can  be  made 
for  other  targets  and  observation  conditions.  The  problem  with  target 
modeling  (as  a  solitary  method)  is  that  one  is  never  certain  that  the  modeling 
is  adequate,  since  the  scattering  process  is  usually  highly  complex.  Hence 
the  suspicion  could  arise  that  the  derived  polarization  properties  are  mere 
artifacts  of  the  target  models.  Of  course,  similarly,  if  only  measured 
responses  are  analyzea,  the  results  could  be  coincidental  properties  of  the 
particular  targets  or  measurement  conditions.  Our  dual  model ing/measurement 
approach  attempts  to  deal  with  these  pitfalls.  First:  an  effort  was  made  to 
model  the  target  resoonse  accurately,  though,  with  computationally  efficient 
algorithms.  Second:  the  computed  aircraft  model  responses  and  their  derived 
polarization  characteristics  were  compared  with  the  measurements  taken  of 
these  models.  The  agreement  was,  for  the  most,  part  very  good.  The 
Inevitable  discrepancies  were  explained  and  shown  to  have  no  effect  on  the 
general  conclusions.  Conclusions  could  then  be  drawn  with  some  confidence 
from  the  computed  responses. 


One  difficulty  remains:  the  targets  studied  here  are  rather  simple, 
consisting  of  fewer  scattering  centers  than  an  actual  aircraft.  Since  the 
object  of  this  study  is  to  examine  the  applicability  of  polarization  concepts 
to  discrimination  between  actual  aircraft,  one  must  avoid  generalizations 
based  on  results  which  are  artifacts  of  target  (relative)  simplicity.  An 
example  of  this  is  the  occurrence  of  a  dominant  distinctive  response,  which 
for  these  target  models  is  present  at  only  one  aspect  (due  to  the  return  from 
the  leading  edge  of  a  wing).  An  actual  target  may  have  more  of  these  or  may 
have  a  distinctive  response  at  a  strategically  important  aspect,  such  as 
nose-on. 


The  details  of  the  target  modeling  effort  is  collected  in  the 
appendices.  At  the  size/wavelength  scales  of  interest  here,  a  complex  target 
can  be  considered  to  consist  of  non-interacting  scattering  centers.  The 
scattering  centers  are  the  edges  and  specular  surfaces  of  the  target  model, 
and  also  the  join  regions  (inside  corners),  which  give  the  important 
double-bounce  (or  n-bounce)  returns.  Simple  expressions  obtained  fron  PO  and 
GTD  are  used  to  model  the  returns,  when  applicable.  Where  these  expressions 
fail  (such  as  at  grazing  incidence)  they  are  modified  to  obtain  agreement  with 
known  responses  for  plates  or  cylinders.  The  known  responses  include 
published  experimental  results  and  exact  diffraction  integral  computations, 
such  as  the  Sperry  Research  Center  STIE  results.  J.t  was  found  possible  to 
describe  the  polarization  and  aspect  dependent  ref.pcnse  of  scattering  centers 
adequately  with  simple  expression*,  although  some  of  these  are  ad-hoc. 

The  targets  and  their  computed  polarization  characteristics  are  given 
in  the  remainder  of  Section  3.  The  measured  responses  and  the  comparison  of 
measurements  with  computations  are  given  in  Section  4. 
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3.1  DESCRIPTION  OF  THE  AIRCRAFT  MODELS 


The  general  features  of  the  simple  aircraft  models  were  specified  by 
RADC.  Three  targets  were  to  be  studied,  each  at  four  size-scales.  Two  of  the 
targets  were  to  be  similar  in  size  and  shape  with  a  difference  in  detail  such 
as  wing  shape  or  position,  the  third  target  was  to  be  dissimilar.  The  target 
rodels  constructed  (and  computed)  all  consisted  of  a  cylindrical  fuselage  and 
flat  quadrilateral  wings.  They  are  plane  symmetric. 

Models  1  and  2  are  suggestive  of  the  F-5  and  MiG-21,  respectively. 
These  are  roughly  the  same  size  and  differ  in  the  sweep  of  the  wings  (30*  and 
60°,  respectively).  Model  3  is  suggestive  of  the  MiG-25;  it  is  larger,  has  a 
relatively  larger  wing  area,  but  has  a  wing  angle  of  30°  like  Model  1.  The 
size-scales  were  chosen  to  correspond  to  frequencies  ranging  from  350  MHz  to 
1.9  GHz. 


The  target  models  are  illustrated  in  Figure  18.  The  dimensions  shown 
are  the  approximate  full  scale  dimensions  of  the  aircraft.  The  dimensions  of 
the  models  are  summarized  in  Tables  3  and  4.  These  dimensions  appear  in  the 
headers  of  the  computer-generated  plotted  result^.  Size-scale  A  corresponds 
to  what  is  considered  to  be  the  lowest  practical  radar  frequency.  Even  at 
this  low  frequency,  the  targets  are  many  wavelengths  long.  Scales  B,  C,  D 
correspond  to  progressively  higher  frequencies. 

It  will  be  seen  that  the  general  characteristics  of  these  targets  are 
similar  at  al i  the  frequencies  considered;  the  main  difference  being  the  rate 
of  fluctuation  with  changi.ic  aspect.  No  new  results  are  expected  at 
higher  equivalent  frequencies  than  size-scale  D.  In  fact,  all  the  important 
conclusions  of  this  study  could  be  drawn  from  the  results  of  sizes  A  and  B. 
Computations  were  also  made  for  a  lower  ka  ^  2.4  (corresponding  to  160  MHz 
for  a  full  scale  model).  This  is  considered  near  the  lower  limit  of 
applicability  of  the  scattering  center  model  for  these  targets.  (No 
measurements  were  made  at  this  size  scale,  since  the  SNR  would  be  too  low). 


TABLE  3  TARGET  DIMENSIONS  IN  UNITS  OF  CYLINDER  RADIUS,  a 


r\ 


T 

B 


z 


Model 

H'a 

B/a 

DZ/a 

Zl/a 

tf;(degees) 

1  (F-5) 

10. 

4.71 

4.0 

-0.3 

30 

2  (M16-2J.) 

10. 

3.65 

7.4 

2.0 

60 

3  (MIG-25) 

10. 

5.85 

6.3 

-0.3 

30 

A. 

Jj 


H 


fct 


TABLE  4  TARGET  SIZE  SCALES 


Full  Scale 

Model  Dimension* 

Size  Scale 

Frequency 

Model 

k& 

£0  a  (cm) 

1 

5.24 

14.3 

A 

350  MHz 

2 

5.24 

14.3 

3 

8.14 

22.2 

1 

8.14 

22.2 

B 

550  MHz 

2 

8.14 

22.2 

3 

12.8 

34.9 

1 

12.8 

34.9 

C 

870  MHz 

2 

12.8 

34.9 

3 

18.6 

50.8 

1 

18.6 

50.8 

D 

1.27  GHz 

2 

18.6 

50.8 

3 

27.9 

76.2 

♦Model  dimension  reference  (cylinder  length)  for  measurement  at  35  GHz 

l 
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3.2  COMPUTED  AIRCRAFT  CHARACTERISTICS  AT  0*  ELEVATION 


We  will  first  examine  the  polarization  characteristics  of  Model  1  at 
the  low  frequency  ka  «  2.4  .  The  parameters  r.i,  2v,  Zy  are  plotted  versus 
aspect  In  figure  19.  The  target  Is  plane  symmetric  at  0*  elevation,  so  that 
the  scattering  matrix  al agonal  and  the  three  parameters  completely  describe 
the  target.  The  cross-poi  null  polarizations  (one  of  which  Is  the  co-pci  max 
polarization)  are  at  the  top  and  bottom  of  the  Poincare  sphere.  The  co-pol 
null  polarizations  (see  Section  2)  are  given  by  the  parameters  *2v  and  Zy  , 
which  are  respectively  the  longitude  and  co-latitude  of  the  point  on  the 
Poincare  sphere.  The  magnitude,  m.  Is  the  square  root  of  the  co-pol  max  power. 

From  nuw  on,  the  parameters  m,  2v,  2y  will  be  plotted,  as  In 
Figure  19,  as  functions  of  aspect  (0*  to  90*),  rather  than  plotted  as  a  locus 
on  the  Poincare  sphere  or  polarization  chart.  The  reason  is  that  the  locus 
fluctuation  Is  so  rapid  it  is  difficult  to  follow  Its  detail  on  the  Poincare 
sphere.  A  disadvantage  Is  that  the  plot  of  the  angle  2\>  shows 
discontinuities  as  It  passes  through  +  or  -180*.  It  should  be  remembered 
that  2\>  actually  rotates  continuously  and  smoothly  around  the  sphere. 

Model  1  consists  of  a  cylinder  (of  which  the  polarization 
characteristics  are  Illustrated  in  Figure  12)  with  the  addition  of  wings  at  a 
30*  rake  angle.  The  magnitude,  m,  has  additional  bumps  on  It  due  to  the 
wings,  compared  with  cylinder  alone.  The  angle  Zy  snows  considerable 
deviation  from  the  cylinder  alone  case,  particularly  at  those  aspects  where 
the  relative  contribution  of  the  wings  to  the  response  Is  strong  (<*  ^30*), 

The  wing  has  a  horizontal  return,  hence  the  polarization  maximum  is  at  H 
polarization  (top  of  sphere),  and  the  co-pol  nulls  are  In  the  lower  hemisphere 
( 2y  <  90*).  This  last  feature  thus  Identifies  a  horizontal  edge.  For 
o  >  40*  ,  the  target  looks  like  a  pure  cylinder.  It  should  be  noted  that  the 
rapid  little  spikes  In  Zy  occur  at  minima  in  the  target  response  (as  seon  In 
the  plot  of  m),  They  occur  when  the  cylinder  H  return  happens  to  hit  a 
deep  null  and  the  (small)  response  Is  temporarily  mostly  vertical.  The  effect 
of  the  wing  on  the  co-pol  null  polarization  is  complicated,  slrce  It  involves 


the  relative  phases  of  the  contributions  from  wing  and  cylinder.  When  the 
wing  contribution  is  not  much  larger  than  that  of  the  cylinder.  Its  effect  on 
the  null  polarizations  is  unclear.  Only  when  the  contribution  of  the  wing  is 
very  large  (as  seen  by  a  large  bump  in  m)  is  there  a  distinctive  (i.e., 
downward)  change  in  2y  . 

For  comparison,  the  characteristics  of  wings  alone  are  given  in 
Figure  20.  The  wing  response  at  e  »  0  is  purely  horizontal,  hence  the 
co-pol  null  is  at  V  (2?  =  0)  for  all  aspects.  (2v  is  undefined  and  here 
plotted  as  zero). 

As  was  already  noted  in  the  dumbell  investigation  of  Section  2,  the 
polarization  characteristics  of  the  compound  target  bear  little  if  any 
relationship  to  those  of  its  component  parts,  except  when  one  of  the 
scatterers  is  dominant.  The  characteristics  change  rapidly  with  aspect;  there 
are  no  invariants  with  aspect. 

The  response  of  the  wings  alone  for  Model  2  is  shown  in  Figure  21. 

The  leading  edge  of  these  wings  have  60*  rake,  causing  a  peak  in  the  response 
at  60*  aspect.  The  response  at  0*  is  due  to  the  trailing  edge  return. 

In  Figure  22  the  characteristics  of  Model  2  are  plotted  side  by  side 
with  those  of  Model  1.  We  see  that  the  characteristics  of  the  two  targets  are 
different  near  30*  and  60*,  where  the  respective  wings  have  an  appreciable 
effect.  One  could  utilize  this  difference  for  target  discrimination  if  the 
aspect  were  known.  Since  2v  varies  more  rapidly  with  aspect,  the 
characteristic  2*  appears  to  be  a  better  candidate  for  use  in  target 
discrimination.  We  note  that  even  in  these  regions  where  2y  show  a 
distinctive  departure  from  its  average  value,  it  Is  still  a  fluctuating 
quantity,  the  result  of  relative  phase  dependent  Interference  of  multiple 
scattering  centers. 
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FIG.  22  Polarization  characteristics  vs.  aspect  ot  model  1  and  model  2;  ka  *  2.4;  0*0 


The  frequency  scales  of  interest  in  this  study  were  higher  than 
considered  so  far,  namely  ka  ■  5.24  and  higher.  We  now  present  the  computed 
null  characteristics  for  all  the  aircraft  models  and  some  cylinders  at  all 
size  scales,  for  elevation  e  -  0  ,  in  Figures  23  through  36.  Let  us 
concentrate  first  on  Figures  23  through  26,  which  are  for  size  scale  A; 
respectively,  cylinder  1  or  2  only  (23),  Model  1  (24),  Model  2  (25),  and  Model 
3  (26).  The  latter  has  a  larger  cylinder,  which  is  actually  represented  by 
Figure  27,  by  coincidence  of  our  frequency  choices. 

The  cylinder  (Figure  23)  has  a  2y  characteristic  which  averages 
near  90*  (meaning  that  the  HH  and  VV  responses  are  about  equal  on  the 
average).  The  angle  2v  sometimes  makes  short  excursions,  sometimes  varies 
continuously  around  the  sphere.  The  co-pol  null  locus  plotted  on  the  Poincare 
sphere  for  aspects  between  25*  and  40*,  a  region  of  small  2\>  variation, 
would  look  like  Figure  16  (a  spiral);  the  locus  for  aspects  between  65*  and 
80*,  a  region  of  continuous  2v  variation,  would  look  like  Figure  13d 
(continuous  encirclement  of  the  sphere  with  up  and  down  oscillation).  The 
addition  of  the  Model  1  wings  (Figure  24)  changes  the  characteristics 
profoundly.  However,  aspect  has  to  be  well-known  to  distinguish  between  the 
cylinder  and  Model  1  on  the  basis  of  a  few  measurements:  Any  value  of  2v 
or  2 y  which  is  taken  on  by  the  cylinder  is  also  taken  on  by  Model  1  at  a 
nearby  aspect.  The  exception  is  near  30*  aspect,  where  the  leading  wing  edge 
is  dominant.  The  above  remarks  also  apply  to  Model  2  (Figure  25),  except  that 
the  distinctive  region  occurs  at  60*.  Model  3  (Figure  26)  is  larger  and 
consequently  the  null  characteristics  vary  more  rapidly  with  aspect.  Again,  a 
distinctive  region  for  2?  is  noted  near  o  *  30*,  corresponding  to  the  wing 
rake  angle.  The  characteristics  of  Model  1  and  Model  3  are  not  readily 
differentiable  except  by  the  rate  of  fluctuation. 

At  the  higher  frequencies  (size  scales  B  through  D  of  Figures  27 
through  36),  the  above  remarks  apply  virtually  without  change.  The  effect  of 
the  wings  stands  out  more  clearly  in  the  plots  of  m  and  2y  (the  bumps  at 
a  »  30*  in  Figure  36,  for  example),  which  is  due  to  the  fact  that  the  wing 
edge  response  Is  independent  of  frequency,  while  the  oblique  cylinder  response 
goes  down  by  the  square  root  of  the  frequency. 


-65- 


KA-  5.24  H.B.SZ.Zl*  lA.ee 


CHARACTERISTICS  H,  2HU.  2&A  US  ASPECT  FOR  cyUMOii* 


KA-  S .24  H.B.DZ.Zl*  16.66  3.86  7  46  3  M  BE1.8C2 


i 


« 


a 


X 


«-  r>  •  ID  CO  ft!  *4«OR*AI<*« 


FIG.  25  Polarization  characteristics  vs.  aspect  of  model  2,  Scale  A;  ka  =  5.24;  0  =  0‘ 
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It  is  apparent,  particularly  at  the  higher  frequencies,  that  at 
almost  all  aspects  the  polarization  characteristics,  as  they  are  presented 
here,  are  useless  for  discriminating  between  these  targets  (consider  Figures 
35  and  36,  for  example). 

3.3  COMPUTED  AIRCRAFT  CHARACTERISTICS  AT  10°  ELEVATION 

For  non-zero  elevation  angle,  an  aircraft  model  no  longer  appears 
plane  symmetric  to  the  radar.  Hence  the  scattering  matrix  is  not  diagonal 
(i.e.,  there  is  coupling  H-V).  The  polarization  maximum  no  longer  occurs  at 
one  of  the  poles  of  the  Poincare  sphere.  Five  rather  than  three  parameters 
are  needed  to  characterize  the  target. 

The  easiest  way  to  deal  with  this  added  complexity  is  to  find  first 
the  co-pol  maximum  polarization,  characterized  by  the  polarization 
parameters  4'  ,  x  (recall  that  4*  ,  x  are  the  orientation  and  ell iptici ty  of  the 
polarization  ellipse);  then  to  apply  the  transformation  to  the 

scattering  matrix,  which  diagonalizes  it.  This  transformation  corresponds  to 
the  rotations  24*  (about  s^)  and  2x  (about  s^)  of  the  Poincare  sphere, 
which  places  the  co-pol  maximum  at  the  top,  and  does  not  change  the  location 
of  the  co-pol  nulls  relative  to  the  co-pol  max.  We  then  obtain  m,  2v,  and 
2y  as  before  for  the  transformed  scattering  matrix.  This  theory  was 
described  in  Section  2.10.  [A  detail  is  the  convention  that  we  restrict  24< 
as  well  as  2x  to  between  ±90°.  The  co-pol  max  is  then  rotated  to  either 
top  or  bottom  of  the  sphere  and  2y  has  the  range  0°  to  180°.] 

As  has  been  described  in  Section  2,  the  rotation  by  241  on  the 
Poincare  sphere  is  the  same  as  rotation  by  4*  about  the  RL0S,  or  a  simple 
reorientation  of  the  radar  reference  plane.  Now,  rol 1- symmetric  targets  at 
any  orientation  in  3-space  have  a  plane  of  symmetry  which  contains  the  RL0S 
(as  one  can  convince  oneself  by  holding  up  a  cylinder  and  cocking  one’s 
head).  For  such  targets,  x  =  0  (i.e.,  no  x  -transformation  is  required  to 
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diagonalize  the  matrix).  Evidently,  2t  is  a  measure  of  the  asymetry  of  the 
target.  These  observations  on  t  were  made  by  Huynen  [14].  However,  it  will 
be  seen  that  2t  Is  also  not  an  Invariant  of  the  target  under  changes  of 
aspect  or  frequency. 

In  Figure  37  the  angles  24'  and  2t  are  plotted,  in  addition  to 
m  ,  for  Models  1  and  2,  at  the  low  frequency  ka  «  2.4,  at  elevation  angle 
o  *  10*.  The  corresponding  characteristics  m  ,  2v,  2y  are  plotted  in  the 
usual  manner  in  Figure  38.  The  fluctuation  of  24*  and  2t  is  seen  to  be 
even  more  erratic  than  that  of  2v  and  2y.  This  has  been  found  to  be 
typical  of  compound  targets.  It  can  be  understood  by  the  following: 

So  long  as  2y  is  unequal  to  90*,  the  co-pol  max  polarization  is  a  point  on 
the  Poincare  sphere,  which  wanders  about  slowly  for  (very)  small  changes  in 
aspect;  then  as  2y  becomes  equal  to  90",  the  co-pol  max  "flashes"  into  a 
great  circle  (*s  discussed  in  Section  2.5);  upon  coming  out  of  this 
degeneracy,  the  co-pol  max  has  switched  from  one  hemisphere  to  the  other.  All 
the  while  the  co-pol  nulls  are  well  defined  and  moving  smoothly.  [Note  tnat 
the  180*  jumps  in  24'  are  artifacts  of  our  definition  of  \v\  <  90*,  ¥ 
actually  changes  continuously  and  smoothly.  It  remains  true,  however,  that 
near  the  switch  of  from  one  hemisphere  to  the  other,  that  motion  is  very 

rapid.  This  phenomenon  occurs  frequently:  whenever  the  co-polarized  return 
switches  from  being  larger  to  smeller  than  the  orthogonal  co-polarized  return.] 
Similar  behavior  is  exhibited  by  2x  .  Because  of  this  erratic  behavior, 

24»  and  2t  (or  the  co-pol  max  locus)  are  not  useful  for  target 
discrimination  for  compound  targets. 

The  2v  and  2y  plots  for  the  e  -  10*  case  have  properties 
similar  to  the  o  =  0*  case. 

In  Figures  39  through  46,  these  characteristics  are  given  for  the 
cylinder  and  Models  1,  2,  and  3  at  e  ■  10*  for  size  scale  A.  For  the  sake 
of  much  needed  brevity,  the  plots  for  the  other  size  scales  are  omitted  from 
this  report,  as  they  show  no  new  information. 
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FIG.  48  Polarization  characteristics  vs.  a«»ct  of  modal  3,  Scaia  A;  ka  *  8,14;  8-10* 


3.4  AVERAGE  POLARIZATION  CHARACTERISTICS  OF  AIRCRAFT  MODELS 


The  rapid  fluctuations  In  the  polarization  characteristics  with 
variation  In  aspect  Is  due  to  the  changes  In  relative  phase  between  scattering 
centers.  Similar  fluctuations  at  any  one  aspect  occurs  with  frequency 
variations.  The  fluctuations  are  a  nuisance  for  target  ID  purposes.  We  must 
either  utilize  them  or  eliminate  them  by  averaging.  The  averaging  method  Is 
Illustrated  here  (refer  also  to  Section  2.8). 

It  Is  assumed  that  measurements  are  made  over  a  band  of  frequencies 
and/or  aspect  such  that  at  least  one  cycle  of  a  fluctuation  Is  observed  (the 
bandwidth  required  for  this  Is  just  the  same  as  required  to  resolve  the 
Individual  scattering  centers,  as  when  forming  an  Image  or  range  profile). 

The  powers  of  the  Individual  measurements  are  added  (Incoherent  sum).  That 
Is,  we  measure  the  average  Stokes  vectors  of  the  scattered  radiation.  The 
process  Is  similar  to  the  natural  observation  of  partially  coherent  light.  We 
may  form  the  expected  value  of  this  average.  This  removes  all  rapid 
fluctuations  due  to  relative  orientation  of  the  scattering  centers. 

The  results  for  the  cylinder,  Model  1,  and  Model  2,  for  size  scale  A, 
at  e  »  0  ,  are  shown  respectively  In  Figures  *'7,  48,  and  49  (m  Is  plotted  In 
dBsm  in  these  figures).  These  are  to  be  compared  with  Figures  23,  24,  and  25. 
We  see  that  the  2v  fluctuation  Is  eliminated.  In  fact  2v  >0.  The  rapid 
fluctuations  In  2?  are  eliminated.  What  Is  left  Is  the  distinctive  effect 
of  the  wing  contribution  to  2?  near  the  wing  specular.  The  fact  that 
2v  -  0  for  these  targets  Is  because  all  the  strong  scatterers  are 
single-bounce  and  there  are  no  strong  double-bounce  scattering  centers.  If  a 
strong  double-bounce  scatterer  were  present  over  some  range  of  aspect,  this 
would  show  up  on  the  2v  plot  as  «*n  excursion  toward  *90*. 

It  should  be  noted  {see  Section  2.8)  that  2v  and  2*  here  do  not 
represent  a  co-pol  null,  but  rather  a  co-pol  minimum.  The  effect  of  averaging 
Is  to  depolarize  the  received  power. 
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FIG.  40  Expected  value  of  iwny  polarisation  characteristics;  modal  2.  Scale  A;  ka  ■  5.24; 

0-0°. 


This  preliminary  investigation  shows  that  averaging  may  prove  to  be 
the  method  by  which  polarization  characteristics  become  viable  target 
discriminants.  It  would  require  a  wideband  radar,  hut  simple  incoherent 
processing.  The  adequate  treatment  of  this  approach  requires  a  statistical 
analysis  of  the  measured  polarization  characteristics.  This  is  beyond  the 
scope  of  this  report. 


SECTION  4 


MEASURED  POLARIZATION  PROPERTIES  OF  SIMPLE  AIRCRAFT  MODELS 
4.1  DESCRIPTION  OF  DUAL-POLARIZED  RADAR 

A  simple  yet  versatile  testbed  rader  system  has  evolved  at  Sperry 
Research  Center  to  meet  our  needs  for  target  and  clutter  data  gathering  and 
algorithm  testing.  Figures  50-52  show  block  diagrams  of  the  measurement 
system  (for  which  we  hold  two  patents,  with  others,  pending). 

The  operation  of  the  system  may  be  understood  by  referring  to  Figure 
50.  An  FMCW  (frequency  modulated,  continuous  wave)  radar  (Figure  50a) 
achieves  range  resolution  by  transmitting  a  chirp.  I.e..  a  signal  which  has  a 
frequency  versus  time  resembling  a  sawtooth  (Figure  50b),  mixing  the  received 
signal  with  the  transmitted  signal,  (resulting  In  an  l.f.  signal  (Figure  50c) 
for  which  frequency  Is  proportional  to  target  range),  narrowband  filtering  at 
l.f.,  and  detecting.  If  the  resulting  signal  is  Integrated  over  many  chirps. 
It  contains  broadband  target  Information  (the  transmit  bandwidth  Is  equal  to 
500  MHz  for  our  system);  however.  If  the  detected  video  Is  sampled  each  chirp 
at  the  same  point  In  the  chirp.  It  contains  target  Information  only  at  a 
single  transmit  frequency  (i.e.,  CW  Information).  Either  of  these  modes  is 
possible;  the  latter  was  used  for  measurements  under  this  effort. 

In  Figure  51,  the  FMCW  Sunn  oscillator  (the  transmitter),  together 
with  the  circulator,  antenna,  and  mixer  can  be  recognized  from  the  simple 
single-polarized  FMCW  system  of  Figure  50.  A  3  dB  power  splitter  has  been 
inserted  immediately  after  the  transmitter  (which  breaks  up  the  transmitter 
power  Into  two  branches,  vertical  and  horizontal);  also  added  are 
computer-controlled  attenuators  In  both  branches  and  a  manual  phase  shifter  In 
the  vertical  branch.  If  the  phase  shifter  Is  set  for  90*,  either  left-hand 
circular,  vertical,  or  horizontal  polarized  radiation  can  be  transmitted  under 
computer  control. 
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The  two  *.f.  signals  (the  outputs  of  the  two  mixers)  go  to  a 
dual -polarized  signal  processor  which  has  two  modes  of  operation.  In  the 
"sequential11  mode,  computer  controlled  attenuators  and  phase  shifter  In  the 
receive  signal  processor  are  used  to  change  (from  chirp  to  chirp)  the 
effective  polarization  of  the  receive  antenna.  Independent  of  the  polarization 
of  the  transmit  antenna.  Any  arbitrary  elliptical  receive  polarization  can  be 
synthesized  by  the  controls  shown.  It  Is  shown  in  Section  4.2  that  for  any 
transmit  polarization,  a  complete  set  of  Stokes  parameters,  and  therefore  one 
of  the  columns  of  the  relative  phase  scattering  matrix,  can  be  calculated 
using  detected  power  from  four  chirps  by  suitably  varying  the  receive 
polarization  (viz.,  vertical,  horizontal,  45*  linear,  and  left-hand 
circular).  In  fact,  six  receive  chirps  were  used  for  each  of  two  transmit 
polarizations  (the  above  four,  plus  135*  linear  and  right-hand  circular)  to 
provide  redundant  Information  for  reduction  of  noise  effects. 

Only  the  sequential  mode  was  used  for  the  measurements  In  this 

contract,  sines  the  simultaneous  mode  requires  bilinear  multipliers  (mixers) 

* 

and  precision  square  law  detectors  of  quality  exceeding  that  available  for 
easily  obtainable  l.f.  signal  processing  components. 

Figure  52  shows  the  overall  radar  system  Integrated  with  ancillary 
components  such  as  a  pedestal  with  controller,  a  boreslghted  TV  camera 
(together  with  video  monitor),  and  data  storage  means  (eight-inch  floppy 
diskette).  The  figure  Indicates  the  central  role  that  the  DEC  NINC  computer 
plays  In  system  operation  and  data  recording.  We  used  the  somewhat  slow  BASIC 
MINC  operating  system  for  hardware/ software  development  and  translated  to  the 
faster  FORTRAN  language  for  an  order-of-magnltude  Improvement  In  date-taking 
run  time. 

Table  5  gives  the  parameters  of  the  RF  subsystem. 
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TABLE  5 


Radar  System  Parameters 


Center  Frequency 
Swept  Bandwidth 
Antenna  Bean  Width 
Antenna  Gain 
Transmit  Power 

Receive  Noise  Figure 


35  GHZ 
500  MHz 
2  degrees 
35  dB 

25  MW  (linear  polarization) 
50  MW  (circular  polarization) 
4  dB  OSB 


Figure  53  shows  receive  antenna  output  power  In  dBm  as  a-  function  of 
range  to  the  target  In  meters.  Receive  noise,  transmitter  noise,  ground 
clutter  (assumed  to  be  In  a  -25  dB  sldelobe)  and  two  signal  levels  (radar 
cross  section  or  RCS  of  -40  to  -20  dBsm,  or  10“*  to.  10“*  n?)  are  shown. 

It  Is  clear  that,  especially  for  the  smaller  target  RCS,  a  range  of  20  m  or 
less  Is  desirable  to  ensure  an  adequate  S/C  level.  In  fact,  all  the  data  runs 
were  at  12  m  or  20  m  range.  The  larger  target  sizes  (above  0.1  m  )  should 
remain  adequately  above  both  clutter  and  noise  out  to  100  m  or  more. 


Figure  54  shows  the  measurement  test  geometry.  A  radar  fence  between 
radar  and  support  structure  Is  useful  In  blocking  backscatter  from  the 
pedestal  and  surrounding  ground. 


As  shown  In  Figure  54  for  the  short  range  (12  m)  0*  aspect  runs,  the 
transmitter  was  placed  at  position  Sj  ,  the  fence  at  Sg  ,  and  the 
pedestal  at  S3  .  For  10*  aspect  runs,  the  transmitter  was  raised  2.1  «•  to 
position  .  For  the  long  range  (20  m)  runs,  the  positions  Lj  ,  Lg, 
and  L3  ,  respectively,  were  used  for  transmitter,  fence,  and  pedestal.  In 
every  case,  the  fence  was  carefully  positioned  using  a  large  RCS  target  at  the 
top  of  the  support  cone  so  as  to  simultaneously  minimize  knife-edge- 
diffraction  distortion  of  the  signal  and  sldelobe  scatter  from  the  pedestal. 


The  target  support  pylon  (Figure  55)  was  a  comical  polyfoee  block 
shaped  so  that  Its  RCS  was  below  -50  dBsm.  The  polyfoam  cylinder  rested  on  a 
remote-controlled  pedestal  which  rotated  It  In  azimuth.  A  14-bit  digital 
shaft  encoder  (Figure  52)  reported  azimuth  angle  to  0.03*  accuracy* 
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FIG.  64  Twt  gaomatry  for  polarization  maaauramants. 


FIG.  55  Poiyfoam  support  cona  and  padastai. 


As  suggested  In  Figure  52,  a  central  computer  (a  Otgltal  Equipment 
Corporation  MINC  U)  was  Med  to  control  the  radar  configuration,  acquire  and 
condition  radar  scattering  data,  and  record  these  data  (along  with  calibration 
and  ancillary  data  in  header  blocks)  on  eight-inch  f loopy  diskettes.  The  data 
on  these  diskettes  was  In  the  form  of  square-law  detected  power  for  different 
transmlt-recjlve  polarization  pairs.  These  data  were  then  read  from  the 
floppy  Into  the  Unlvac  1108  where  power  data  were  mapped  into  co-pol  nulls  and 
cross-pol  nulls  for  each  aspect. 

Figures  86-59  show  photographs  of  the  RF  subsystem,  the  receiver 
signal  processor,  the  20  m  range  target  as  viewed  from  the  antenna,  and  two 
targets  mounted  on  polyVoam  cone  tips. 

Since  targets  which  are  symmetric  with  respect  to  earth  vertical  have 
zero  HV  resporse  (i.e.»  zero  off-diagonal  elements  In  the  earth-basis 
scattering  matrix),  horizontal  and  vertical  polarization  are  poor  choices  for 
transmit  polarization  for  relative  phase  measurements  between  the  HH  and  VV 
responses  (the  diagonal  elements  of  the  scattering  matrix).  Figure  56  shows  a 
simple  solution:  tilt  che  radar  45*  so  that  one  transmit  attenuator  really 
controls  linear  polarized  radiation  at  135*  with  respect  to  earth  horizontal, 
and  the  other  attenuator  really  controls  45*  linear  polarization.  These  two 
polarlz  .tioi.s  w“»re  referred  to  as  V  prime  and  H  prime,  respectively. 

Figure  87  Indicates  the  receive  signal  processor  layout,  with  60  MHz  l.f, 
components  (Avantea  voltage-controlled  amplifiers,  bandpass  filters,  3  dB 
power  splitters,  Olektron  phase  shifter,  and  RHS  log  detector)  on  the  lowest 
level  and  video  processing/ sampling  circuitry  ',n  the  uppermost  level. 

"i’jure  58  Indicates  the  relative  geometry  of  the  polyfoam  support, 
target,  radar  fence,  and  surrounding  ground,  as  viewed  by  the  radar. 

Figure  59  shows  two  targets,  scales  A  and  0  respectively,  sitting  on 
their  polyfoaut  cone  tips. 


-108- 


B2-M7 

FIG.  56  Radar  tranamittar  and  antanna. 
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FIG.  57  Radar  racaivar/signal  proeaaaor. 


82-8841 


FIG.  58  Support  pylon  and  target  viewed  from  radar. 
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FIG.  59  Two  target*  (Scale*  A  and  C)  and  cone  tip*. 


-107- 


4.2  MEASUREMENT  OF  STONES  VECTORS 


The  null  end  maximum  polarization*  and  the  Stokes  parameters  (SR)  are 
both  sets  of  polarization  descriptors  that  completely  describe  the  elements  of 
the  relative  phase  polarization  scattering  matrix  (RPSM)  of  an  object.  The 
polarization  nulls  and  maxima  have  the  advantage  that  they  are  Invariant  with 
respect  to  a  radar's  transmit  and  receive  polarization  for  a  given  target 
orientation.  The  SR  do  not  have  this  advantage.  However,  since  In  a 
practical  radar  system  the  transmit  and  receive  polarizations  are  usually 
known,  and  It  Is  the  elements  of  the  RPSM  that  are  measured;  and  since  the 
relation  between  these  measurements  and  the  SP  Is  more  apparent  than  that 
between  the  measurements  and  the  nulls  and  maxima,  the  focus  here  will  be  on 
the  SP. 


The  measurements  for  the  target  sets  required  to  satisfy  this 
contract  consist  of  two  groups:  one  at  a  0*  radar- to- tar get  depression  angle 
and  another  at  approximately  a  10*  degression  angle.  In  the  0*  depression 
angle  case,  the  aircraft  models  are  plane-symmetric  In  that  they  lie  In  a 
plane  of  symmetry  relative  to  the  radar  resulting  In  a  polarization  scattering 
matrix  which  has  only  diagonal  terms  that  are  In  general  non-zero.  Since  the 
RPSM  Is  desired.  It  Is  clear  that  three  real  numbers,  two  magnitudes  and  a 
relative  phase,  completely  describe  the  diagonal  matrix  for  this 
plane-symmetric  target  case: 


RPSM 


i*n  i 


l«22  l«J“ 


(40) 


where  a.,  and  a-9  are  the  complex  elements  of  the  scattering  matrix  and 

«U  \  J*i  2 


*11  *  I  *11  le  »  *22  ■  1  *22  I* 


and  ap  -  #22-  #u  . 


It  was  found  that  a  tingle  transmit  polarization  (135*  linear  was 
chosen)  and  a  minimum  of  four  receive  polarizations  are  sufficient  to 
determine  the  elements  In  the  matrix  above.  A  theoretical  minimum  of  four 
receive  polarizations  Is  necessary  Instead  of  three,  to  ensure  that  the 
relative  phase  Is  measured  unambiguously.  That  Is,  since  ad  Is  measured  as 
either  sin  ad  or  cos  ad  ,  both  must  be  measured  to  unambiguously  determine 
a#  .  Practically  It  was  decided  to  measure  the  power  In  six  receive 
polarizations:  V  and  H  linear,  135*  and  45*  linear,  and  LH  and  RH 
circular.  These  received  powers  are  denoted  by  P^,  P^,  P3,  P4,  Pg, 
and  Pg,  respectively.  The  SP  for  a  received  electric  field  may  then  easily 
be  constructed  as  follows: 

s0  ■  V(P2  -  P1)Z  *  (P4  -  P3>Z  ♦  <P6  -  PS>2 
sl's0  ’  (P2  *  Pi),s0 

(41) 

s2/s0  -  (P4  -  P3)/s0 

S3/SQ  ■  (Pg  ~  Pg)/$o  • 

This  approach  ensures  that  the  SP  s^,  and  s3  sum  square  to  sQ 
and  thus  that  the  degree  of  polarization  is  one  as  It  should  be  for 
monochromatic  electric  fields.  The  SP,  three  of  which  are  Independent, 
completely  describe  the  three  numbers  that  represent  the  symmetric  target  RPSM 
from  which  the  polarization  nulls  and  maxima  and  In  particular  m,  2y»  and 
2v  are  derived.  See  Appendix  8.10  for  a  description  of  the  relationship 
between  the  SP  and  the  elements  of  the  RPSM. 


-109- 


In  the  10*  depression  angle  case,  the  aircraft  models  do  not  lie  In  a 
plane  of  symmetry  relative  to  the  radar,  resulting  1h  a  RPSH  all  of  whose 
elements  are  In  general  non-zero.  That  Is*  for  such  an  asymmetric  target 


•u  ie 


where  once  again  the  complex  elements  of  the  scattering  matrix  are  a^^  «|a^j|e 


a12  "  la12'c 


a22  -  la22 I e  *  and  A*ll  “  *11  ‘  *12  and 


a*22  "  *22  ~  *12  *  F’ve  real  numbers  completely  describe  this  matrix. 

In  order  to  measure  these  five  nurr.bc  ;,  it  was  found  that  at  least 
two  transmit  polarizations  are  necessary.  The  two  chosen  here  are  the 
orthogonal  pair  45*  linear  and  135*  linear.  Once  again,  as  in  the  symmetric 
target  case,  it  was  decided  to  measure  the  powers  in  six  receive  polarizations 
V  and  H  linear,  135*  and  45*  linear,  and  LH  and  RH  circular,  for  each 
of  the  transmit  polarizations  resulting  in  a  tot  .1  of  twelve  power 
measurements  at  each  target  aspect.  Some  of  these  measurements  are  redundant. 
However,  by  combining  them  as  was  done  in  Equation  41,  it  can  be  ensured  that 
the  normalized  SR  for  each  transmit  sum  square  to  one  independently  of  each 
other.  Thus  all  twelve  measurements  are  used  to  determine  the  five  real 
numbers  that  completely  describe  the  RPSM  (see  Appendix  8.11). 


Two  different  minimal  sets  of  seven  measurements,  one  for  each;  of  the 
three  magnitudes  desired  |  axx  I  >1  ai2  !  v  and  I  a22  I  and  two  each 
avoid  any  sign  amgibulty)  for  each  of  the  relative  phases  desired,  were  also 
found.  It  was  not  clear,  however,  that  using  one  of  these  minimal  sets  of 
measurements  was  better  than  using  all  twelve  measurements,  which  has  the 
advantage  of  decreasing  the  effects  of  noise  and  sampling  jitter. 
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4.3  POLARIZATION  CALIBRATION 

The  calibration  of  the  SRC  dual-polarized  radar  Mas  separated  Into 
two  parts.  The  first  part  consisted  of  the  calibration  of  the  Individual 
components  In  the  system*  Including  the  Gunn  oscillator  power  source,  the 
log-detector,  the  transmit  and  receive  attenuators,  the  receive  amplifiers, 
and  the  transmit  and  receive  phase  shifters.  The  second  part  consisted  of  the 
system  polarization  calibration,  which  ensuied  that  Indeed  the  transmit  and 
receive  polarization  was  that  which  was  desired.  The  component  calibration  is 
described  further  In  Section  4.4.  The  focus  In  the  remainder  of  this 
subsection  will  be  on  the  system  polarization  calibration. 

A  one-way  calibration  procedure  was  performed  to  ensure  that  the 
transmit  polarization  was  as  desired.  Since  135*  linear  and  45*  linear 
polarizations  are  the  desired  transmit  polarizations.  It  was  decided  that 
these  two  polarizations  could  be  best  transmitted  by  rotating  the  radar  45* 
relative  to  earth  horizontal,  and  transmitting  vertical  or  horizontal 
polarization  (in  radar  coordinates).  If  one  thinks  In  terms  of  a  V,  H 
coordinate  system  for  the  target,  and  a  V',H'  coordinate  system  for  the 
radar,  where  V',H'  Is  rotated  by  +45*  relative  to  V,H  ;  then  V  transmit 
alone  Is  135*  linear  polarization  and  H'  transmit  alone  Is  +45*  linear 
polarization  In  the  target  coordinate  system.  This  procedure  allows  the  two 
desired  transmit  polarizations  to  be  realized  with  one  transmit  channel  ON 
while  the  other  Is  completely  OFF.  The  effects  of  cross-polarization  leakage, 
frequency  pulling  on  the  Gunn  oscillator,  and  uncertainties  In  the  path 
lengths  of  the  two  channels  which  affect  the  relative  phase  and  hence  the 
polarization  of  the  transmitted  signal,  are  thus  minimized.  Each  transmit 
attenuator  may  then  be  calibrated  Individually  to  ensure  that  the  power  is 
Identical  for  each  transmit  polarization,  45*  end  135*  linear. 
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The  receive  attenuators  and  amplifiers  were  Initially  calibrated 
Individually  with  a  laboratory  reference  60  MHz  Input  signal.  In  order  to 
calibrate  the  entire  system  receive  polarization;  that  Is,  the  receive 
amollflers  and  phase  shifter  which  determine  the  radar  receiver  polarization, 
a  two-way  automatic  calibration  procedure  was  developed.  Basically  this 
procedure  consists  of  Illuminating  a  known  cress  section,  1  square  meter, 
dihedral  corner  reflector  with  the  crease  oriented  at  22.5*  relative  to 
horizontal  with  either  a  45*  or  135*  linear  polarization  (-the  transmit 
polarizations  used  to  measure  the  elements  of  the  RPSM).  This  calibration 
target  returns  a  vertically  or  horizontally  polarized  signal,  respectively; 
either  of  which  has  equal  power  In  both  the  V  and  H'  channels  of  the 
radar  receiver.  Knowing  that  this  Is  the  case,  the  electronically  controlled 
digital  phase  shifter  may  then  be  stepped  through  all  possible  phases  under 
computer  control,  and  the  resulting  log-detected  power  stored  for  each  phase 
shift.  The  phase  shift  corresponding  to  maximum  received  power  must  then  be 
0*  relative  phase  between  the  two  channels,  while  the  minimum  received  power 
corresponds  to  a  180*  phase  shift  between  the  two  channels.  All  other  phase 
shifts,  and  In  particular  those  allowing  LH  and  RH  circular  polarization 
reception,  may  be  Inferred  from  these. 

The  power  difference  between  the  maximum  and  minimum  received  power 
gives  an  indication  of  the  balance  of  the  amplification  In  each  receive 
channel.  If  this  difference  (measured  In  decibels)  Is  small.  It  can  be 
Improved;  that  is,  made  larger  by  adjusting  one  of  the  receive  amplifiers  to 
balance  or  equalize  the  amplification  In  both  channels.  The  automatic  phase 
calibration  program  may  then  be  run  once  more  producing  the  phase  settings 
required  to  configure  the  receiver  to  the  receive  polarizations  desired. 

The  power  difference  between  the  maximum  and  minimum  received  power 
also  gives  an  Indication  of  the  polarization  purity  of  the  receiver.  It  was 
found  that  the  cross-polarized  power  was  from  25  dt>  to  30  dB  below  the 
co-pol ari zed  power. 
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Since  the  calibration  target,  the  dihedral  corner  reflector,  is  a  0 
dBsm  or  1  square  meter  reflector  at  35  GHz;  it  was  also  used  to  determine  a 
conversion  factor  relating  the  power  measured  by  the  log  detector  In  dBm  to 
the  actual  target  size  In  units  of  dBsm.  Comparison  of  predicted  and  measured 
responses  In  Section  4.6  shows  this  power  calibration  to  be  fairly  accurate. 

The  accuracy  of  this  two-way  calibration  procedure  strongly  depends 
on  the  accuracy  of  the  angles  involved  In  setting  up  the  calibration 
experiment.  The  radar  should  be  tilted  at  exactly  45*  relative  to  horizontal 
and  the  dihedral  should  be  tilted  at  exactly  22.5*  relative  to  horizontal.  A 
bubble  level  with  an  accuracy  of  *0.5*  was  use<1  to  set  these  angles  and  was 
found  to  be  acceptable. 

4.4  OTHER  CONCERNS 

Even  after  careful  polarization  calibration  there  are  several 
potential  sources  for  error  In  a  radar  system  measuring  polarization  effects 
which  sometimes  can  be  of  very  small  magnitude.  These  potential  error  sources 
can  be  groupedlnto  radar-system- Induced  errors,  errors  associated  with 
geometry,  errors  due  to  the  construction  of  the  target,  and  errors  from 
ambient  clutter. 

4,4.1  System 

System  errors  come  frcm  three  major  sources:  additive  noise,  system 
nonl Inearl  ties,  and  drift. 

Transmitter  noise,  which  dominates  receiver  noise  In  our  system,  (see 
Figure  53),  is  due  to  FM  sidebands  on  the  carrier  of  the  Gunn  diode 
transmitter,  which  leak  through  the  system  and  appear  at  the  l.f.  frequency  . 
(In  our  system,  the  major  leakage  path  is  a  bounce  off  the  subref lector  of  the 
Cassegrain  antenna).  By  using  radar  absorbing  material  on  the  Subreflector, 
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transmitter  noise  Is  kept  to  12  dB  above  receiver  noise;  averaging  ten  pulses 
further  reduces  both  noise  sources  by  5  dB,  resulting  In  a  noise  "baseline11 
which  corresponds  to  a  target  cross  section  of  -54  dBsm  (six  millionths  of  a 
square  meter)  at  a  range  of  12  m.  This  was  measured  as  part  of  run  P29R01, 
data  from  which  Is  shown  In  Figure  60.  The  total  backscattered  power  (in  dB 
relative  to  a  square  meter)  Is  plotted  as  the  pedestal  turns;  at  aspects  of 
25* »  45*,  50*,  52.5*,  55*,  75*,  and  90*,  the  "target*  was  qu1cM>  replaced. 

The  data  from  0*-52.5*  and  55*-90*  represents  several  varieties  of  polyfoam 
support  cone  tips,  while  the  data  from  52.5*-50*  and  90*-100*  represents 
essentially  system  noise;  note  that  It  has  an  average  level  r '  about  -54  dBsm. 
Actually,  this  "noise"  data  also  includes  ambient  clutter  f\  -:.  the  pedestal 
and  nearby  ground,  which  is  evidently  well  screened  by  the  radar' fence. 

The  other  system-induced  errors,  nonllnearltles  and  drift,  require 
careful  calibration  for  those  aspects  which  are  highly  nonlinear  but  stable 
(e.g.  log  detector  deviation  from  true  linear  volt.s-per-dB  characteristic, 
voltage  controlled  amplifier  gain-versus-volts  characteristic),  and  frequent 
calibration  for  those  aspects  which  indicate  drift  (e.g.  phase/gain 
characteristics  of  the  RF  mixers,  which  potentially  can  change  as  the  Gunn  VCO 
frequency  changes  with  time  or  temperature).  The  stable  nonlinear  elements  of 
the  radar  were  calibrated  both  with  a  CW  60  MHz  l.f. laboratory  source  and  with 
scattered  power  from  a  1  square  meter  reference  dihedral  target.  These 
calibration  runs  were  automated  with  the  MINC  computer  as  master  controller; 
typically  100  measurements  were  averaged  to  drive  measurement  noise  errors  In 
the  calibration  process  to  0.2  dB  or  less.  The  drift-prone  components  of  the 
system  (primarily  the  two  mixers)  were  calibrated  hourly  using  a  special 
dual-polarized  computer  controlled  algorithm  described  in  Section  4.3,  which 
simultaneously  closed-loop  corrected  both  relative  gain  and  relative  phase  in 
the  two  receiver  channels. 

4.4.2  Geometry 

i 

Potential  error  sources  associated  with  geometry  Include  target 
orientation,  antenna  orientation,  placement  of  the  radar  fence,  and  near  field 
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effects.  Target  aspect  was  controlled  via  the  following  chain:  the  MINC 
computer  interfaced  a  motor  controller*  which  determined  the  direction  and 
speed  of  a  dc  motor,  which  was  linked  by  a  pulley  to  a  synchro  driver.  The 
synchro  receiver  in  the  pedestal  then  turned  the  polyfoam  support  and  target. 
Any  backlash  or  other  angle  errors  would  be  measured  by  a  synchro-to-digital 
converter  which  was  accurate  to  0.03*.  Target  pitch  (which  was  nominal  zero 
degrees)  was  carefully  measured  each  run  with  bubble  levels  accurate  to  0.5*. 

Antenna  orientation  was  facilitated  with  the  use  of  a  co-mounted  TV 
camera  with  remote  monitor.  The  camera  was  carefully  bore-sighted  using  a 
physically  small  target  with  large  RCS  (the  4.1-c»-wide  dihedral  with  1  square 
meter  RCS)  at  the  range  of  interest  (either  12  m  or  20  m).  Later,  the 
radar-TV  combination  could  be  pointed  at  a  small  RCS  target  with  0.1* 
accuracy,  ensuring  that  the  peak  of  the  2*  beam  is  on  target.  The  45  tilt  of 
the  radar  (Figure  56)  was  ensured  using  a  bubble  level. 

The  radar  fence  (Figures  54  and  58)  was  carefully  placed  to  minimize 
the  sldelobe  response  of  the  metal  pedestal,  while  at  the  same  time  minimizing 
distortion  of  the  target  response  due  to  knife-edge  scatter  over  the  top  of 
the  radar  fence. 

Near  field  effects  (see  Section  8.9)  were  deemed  to  be  serious  enough 
at  12  m  range  for  the  largest  three  target  scales  (C,  D,  E)  that  runs  were 
made  at  20  m  range,  even  though  the  S/N  ratio  decreased  9dB  by  moving  out  to 
that  range. 

4.4.3  Target 

Potential  target- Induced  errors  Included  unwanted  scattering  from 
mounting  screws,  wing  slots,  and  Incorrect  assembly.  The  screw  and  slot 
errors  were  minimized  by  covering  them  with  conductive  tape,  which  was 
carefully  smoothed.  The  major  assembly  problem,  wing/body  alignment,  was 
identified  in  an  early  test  run  when  the  local  RCS  peak  due  to  the  wing 
leading  edge  appeared  at  an  azimuth  of  32*  Instead  of  30*;  for  all  runs  after 
this,  wing  alignment  accurate  to  0.5*  was  ensured  by  careful  assembly. 


4.4.4  Clutter 


Unwanted  responses  from  scettcrers  near  the  target  (1 .e.,  rider 
clutter)  dominated  the  Interference  at  the  12  a  range  (see  Figure  S3),  while 
at  20  m,  transmitter  noise  Is  dominant.  The  major  contributor  to  clwtter  Is 
the  polyfoam  support  cone;  of  this,  the  tip  Is  dominant.  Other  scatterers 
Include  the  pedestal  and  nearby  ground,  both  of  which  ere  reduced  to  below 
transmitter  noise  (even  at  12  m  range)  by  the  radar  fence.  Scatterers  at  long 
range  (for  example  the  tower  visible  In  Figure  58)  are  easily  range-gated  out, 
using  the  FMCM  3  MHz  range  gate  filters  (labelled  "NBF*  for  "narrowband 
filter"  in  Figure  51). 

The  upper  12-Inch  tip  of  the  60-Inch  polyfoam  support  was  detachable; 
In  fact,  several  tips  were  fabricated.  Two  types  are  shown  In  Figure  59.  The 
cone  tip  on  the  right  had  three  sections,  held  together  with  a  heavy 
application  of  low  RCS  glue,  giving  a  relatively  poor  RC$  of  -37  dBsm  to  -44 
dBsm,  as  shown  in  Figure  60  (0*-25*  azimuth).  Other  three-section  tips  had 
better  performance  (below  -45*  dBsm  typically,  as  shown  In  Figure  60, 

25*-50*).  The  best  low  RCS  cone  tip  consisted  of  only  one  section  with  no 
glue,  as  shown  to  the  left  In  Figure  59;  this  cone  tip  had  an  RCS  of  -50  dBsm. 

An  Indication  of  the  clutter-limited  sensitivity  of  the  system  at  12 
m  range  Is  given  In  Figure  61,  which  shows  polarization  data  for  a -half-inch 
ball  bearing  (with  an  RCS  of  -41  dBsm).  The  ball  bearing  was  removed  b*twt*n 
20*  and  40*  and  between  60*  and  80*,  leaving  only  the  cone  tip  (one  of  the  -45 
dBsm  low-glue  cones).  Note  that  the  presence  of  the  ball  bearing  Is  clearly 
evident  In  both  the  Sg  and  s2  plots,  as  It  should  be,  and  does  not 
appear  In  the  and  plots,  as  It  In  fact  should  not.  (Remember 
that  for  135*  linear  Illumination  and  a  specular  one-bounce  scatterer,  all  the 
scattered  radiation  is  at  135*;  s^  -  0,  Sj  -  0,  s^/Sg  •  -1.) 
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FIG.  61  RMponat  of  half-inch  iphar*  (-41  dftam) 


One  night  ask  just  how  Urge  should  the  slgnal-to-c latter  ratio  be  to 
cause  acceptably  low  distortion  to  a  measured  power  level.  For  signal  power 
S  ,  clutter  power  C  and  relative  signal /clutter  phase  o  the  S+C  power  Is 


p  -  |/S  ♦  /T  e*#  |2  (43) 

which  has  a  maximum  of  S  ♦  C  ♦  2  /?C  and  a  minimum  of  S  ♦  C  -  2 /SC  . 
These  bounds  are  shown  In  Figure  62.  Note  that  20  dB  S/C  ensures  error  less 
than  *1  dB,  while  10  dB  S/f  results  in  errors  between  -3.3  dB  and  +2.4  dB. 
Consulting  Figures  53  and  62,  we  see  that  -50  dBsm  cones  can  cause  *3  dB 
errors  In  -40  dBsm  target  measurements,  but  only  *0.3  dB  errors  In  -20  dBsm 
target  measurements. 


4.5  MEASUREMENT  RESULTS 

In  all,  over  fifty  target  runs  were  recorded.  Table  6  describes  the 
Information  recorded  In  the  disk  file  header  for  each  run,  and  also  describes 
the  file  name  convention  used.  The  first  letter  Is  always  P  ifor 
polarization).  The  next  two  numbers  describe  the  target  body  and 
configuration,  respectively.  (Target  body  scales  A-D  correspond  loosely  to 
body  types  1-5).  The  next  letter  Is  R  for  0*  single- transmit  (symmetric) 
targets;  A  for  10*  dual-transmlt  (asymmetric)  targets;  or  0  for 
dual-transnlt  runs  with  symmetric  targets.  The  last  two  digits  form  the  run 
number. 


Table  7  gives  a  brief  listing  of  the  files,  along  with  the  date 
recorded,  five  header  data  Items,  and  header  comments.  The  header  data  Items 
Include  attenuator  setting  (Identical  switched  attenuators  were  used  before 
both  horizontal  and  vertical  receive  amplifiers  to  preclude  Input  saturation), 
azimuth  step  size  of  the  pedestal  In  degrees,  target  range  In  meters,  and 
calibration  target  responses  In  vertical  and  horizontal  channels. 
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HEADER  INFORHATION**t************m********f ********* 

1  I  LINES  OF  DATA  ON  FILE  AFTER  HEADER 

2  H'  XHT  GAIN  (DB> 

3  V'  XHT  GAIN  (OB) 

4  XHT  PHAGE  SETTING 

5  FIXED  ATTENUATOR  <DB>  IN  H  RCV  PATH  (-100.0) 

6  FIXED  ATTENUATOR  <DB>  IN  V  RCV  PATH  (-100»0> 

7  HAX  H'  RCV  GAIN  (DB> .V'  XHT 

8  HAX  V'  RCV  GAIN  (DB).V'  XHT 

9  RCV  PHAGE  GHIFTER  SEETTING  FOR  IN-PHASE.V'  XHT 

10  TOTAL  AZIHUTH  SCAN  (DEG) 

11  AZIHUTH  STEP  SIZE  (DEG) 

12  AZIHUTH  START  ANGLE  (DEG) 

13  SAHPLE  TIHE  IN  NSEC 

14  NUHBER  OF  SAHPLES  AVERAGED 

15  RANGE  TO  TGT  (H) 

16  IF  PREAHP  OUTPUT  (DBH)  FOR  1  SQ  H  TGT»V'  XHT 

17  RCV  PHASE  SHIFTER  SEETTING  FOR  IN-PHASE.H'  XHT 

18  IF  PREAHP  OUTPUT  (DBH)  FOR  1  SC  M  TGT »H'  XHT 

19  HAX  W'  RCV  GAIN  (DB).H'  XHT 

20  HAX  V'  RCV  GAIN  (DB)'H'  XHT 


FILE  NAHE  CONVENTION********************************* 
PIJR01.DAT 
I  BODY  TYPE 


1  14  CH  AIRCRAFT  (SCALE  A) 

2  22  CH  AIRCRAFT  (SCALE  A-B 

3  35  CH  AIRCRAFT  (SCALE  B-C 

4  43  CH  AIRCRAFT  (SCALE  C-D 

5  64  CH  AIRCRAFT  (SCALE  D) 

6  1  SQ  H  DIHEDRAL 

7  1/2  IN  BB 

8  1  FT  SPHERE 

9  CONE  ONLY 


J  CONFIGURATION 

0  BODY  ONLY  (WITH  SLOTS  AND  TAPE) 

1  30  DEG  THIN  WING 

2  60  DEG  WING 

3  30  DEG  FAT  WING 

4  445  DEG  DIHEDRAL 

5  -22.5  DEG  DIHEDRAL 

6  VERT  DIHEDRAL 

7  HOR  DIHEDRAL 

8  SHOOTH  RCC  AIRCRAFT  BODY 

9  CONE  ONLY 


TABLE  6 

FILE  HEADER  AND  KANE  CONVENTIONS 
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TABLE  7 
TEST  RUN  LOG 


F I Lf * DATE  .ATTEN. AZ1TEF t  RANEE  » DDHV ’ . DDHH ' *  CONNENT 


DISK  4 


F32D02 

lf-JUL-t 

-20 

0.50 

20.000 

-12.700 

-14.200 

LONS  RAMOS 

F43D03 

lt-JUL-l 

-30 

0.50 

20.000 

-12.700 

-14.200 

I.ONO  RANOE  —  RADAR  FENCE 

F31D02 

10-JUL-S 

-20 

0.90 

20.000 

-12.700 

-14.200 

LONO  RANOE 

F42D02 

U-JUL-t 

-30 

0.50 

20.000 

-14.000 

-15.200 

LONS  RANOE  —  RADAR  FENCE 

P41D02 

14-JUL-t 

-30 

0.50 

20.000 

-14.000 

-15.200 

LONO  RANOE  —  RADAR  FENCE 

FS3D02 

U-JUL-t 

-30 

0.50 

20.000 

-14.000 

-15.200 

FIRST  DON  RUN-  RADAR  FENC 

P12D01 

U-JUL-t 

-20 

0.20 

12.000 

-4.400 

-4. tOO 

RADAR  FENCE  —  CALN 

F53D01 

U-JUL-t 

-40 

0.50 

12.000 

-4.400 

-4. *00 

RADAR  FENCE-CALN  DAY 

FILE .DATE, ATTEN. AZSTCP 

.RANOE 

.DDHV'. 

DiHH ' . CONNENT 

DISK  3 

F41D01 

U-JUL-t 

-40 

0.50 

12.000 

-4.400 

-4. tOO 

RADAR  FENCE—  CALN  AIR 

F31D01 

14-JUL-t 

-30 

0.5V 

12.000 

-5.400 

-5.300 

RADAR  FENCE 

P42D01 

14-JUL-t 

-40 

0.50 

12.000 

-5.400 

-5.300 

RADAR  FENCE 

F32D01 

14-JUL-t 

-30 

0.50 

12.000 

-5.400 

-5.300 

RADAR  FENCE 

P43D02 

14-JUL-t 

-40 

0.50 

12.000 

-5.400 

-5.300 

RADAR  FtNCE 

F43D01 

14-JUL-t 

-SO 

0.50 

12.000 

-5.400 

-5.300 

RADAR  FENCE 

PB0D01 

15-JUL-t 

-40 

1.00 

12.000 

-5.000 

-0.100 

HERDIE  IN  NAIN  DEAN 

PttAOl 

13-JUL-i 

0 

1.00 

12.000 

-4.000 

-4.100 

P43A01 

13-JUL-t 

-40 

0.50 

12.000 

-4.000 

-4.100 

STILL  VlNDr 

P42A01 

13-JUL-I 

-40 

0.50 

12.000 

-4.000 

-4.100 

VINDY 

P32A01 

13-JUL-I 

-40 

0.50 

12.000 

-4.000 

-4.100 
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P22A01 

13-JUL-t  -30 
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-5.400 
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A  LITTLE  VINDY 

P11A01 

13-JUL-t  -20 

0.20 

12.000 

-4.700 

-5. tOO 

NOT  VINDY 

P33A01 

13-JUL-t  -30 

0.20 
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-4.700 

-5. tOO 

FIRIT  ROOF  RUN-  HOT  SUN— 

P33R02 

Ot-JUL-t  -40 

0.20 
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P20R01 
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F2SR02 
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0.20 
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12.200 
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FlOROl 
T70R02 
FitROl 
F29R01 
*70R0l 
F2IR01 
F12R04 
FtOROl 
r 28R01 
P23JJ04 
* 33R01 
f 12R03 
r  23R03 
f  23R02 
P44R01 
f  23R01 
f 12R02 
f 43R01 
f 12R01 


30-JUN-I2  NEW  CON  -30  0.50  12.200  -11.200  0.000 

30-JUH-I2  NEU  CONE-  20  ON'  20  OFF  0  1.00  12.200  -11.200  0.000 

It- JON-02  ON  0000  CON  -30  0.90  12.200  -0.000  0.000 

It- JUN-92  CONE  ONL  0  0.20  12.200  -0.000  0.000 


-t.tOO 


•0.000 

-0.000 


HALF  IN  tFHERE  ON  OFF  EVERY  20  FOINTt  15-JUN-t  0  1.00  12.200  -t.tOO 

fCRCV  STICK*  UF  13-JUN-t  -30  0.20  12.200  -0.000  0.000 

NEW  TAFE-  VINOS  SQUARED  1S-JUN-0  -20  0.50  12.200  -0.000  0.000 

VINDY-  ON  TASLE  15-JUN-t  -30  0.10  12.200  -0.000  0.000 

RCC  RUN  ON  tTAND  15-JUN-t  -30  0.20  12.200  -0.000  0.000 

SQUARED  VINOS  -  .1  DEQ.  AZ  STEF  10-JUN-t  -30  0.10  12.200  -0.000  0.00 

FIRST  RUN  -  TAFED  -  3  NHZ  FILTER  03-JUN-t  -30  1.00  12.200  -0.000  0.0 

TAFtO  -  3  NHZ  FILTER  03-JUN-0  -30  1.00  12.200  -0.000  0.000 

10  DEO  TO  20  K6*  .1  DEO  STEFS  U3-JUN-0  -30  0.10  12.200  -0.000  0.000 

TAFE  ON  SCREVt  -  3  NHZ  FILTER  03-JUK-S  -30  1.00  12.200  -0.000  0.000 

REFRODUCE  MEANER  RESULTS  03-JVN-9  -30  1.00  12.200  -0.000  0.000 

FIRST  RUN  VITH  TOT  2-VINDY  AN»  NOT  27-NAY-S  -35  1.00  12.200  -0.000  0 

25  Dt  ATTEN  -  10  SANFLEO  A VO 8  -  20.0  -  22  27-NFY-Q  -20  1.00  12.200  -0. 

TEST  20-NAY-t  -30  1.00  0.000  0.000  0.000 

FIRST  TAROET  RUN  20-MAY-t  -30  1.00  12.200  -4.400  0.000 
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Figures  63  through  76  plot  measured  Stokes  parameter  data  for  targets 
at  0*  aspects  while  Figures  77  through  88  plot  these  data  for  targets  at  10* 
aspect.  As  described  In  Section  4.2,  one  transmit  polarisation  (135*  linear) 
Is  sufficient  for  the  first  set,  while  two  transmit  polarizations  (45*  and 
135*)  are  required  for  the  second  set.  The  model  and  scale  conventions  used 
can  be  Interpreted  by  consulting  Table  4;  for  example.  Figure  70,  which 
represents  Model  3,  scale  B,  In  turn  corresponds  to  -a  MIG-25  at  550  MHz  scaled 
frequency. 

Several  target  polarization  features  emerge  clearly  from  these  data. 
By  far,  the  most  striking  and  potentially  useful  Is  the  horizontal  dipole 
formed  by  the  leading  edge  of  the  target  wing.  The  characteristic  Stokes 
parameter  response  for  such  a  dipole  Is  s^/ sQ  >1,  Sj/Sq  -  0, 

Sj/Sq  -  0.  In  fact,  the  figures  with  Models  1  and  3  tend  to  show 
s^/Sq  near  +1  at  30*.  while  s2  and  s3  both  fluctuate  around  zero 

(e.g.  Figure  76);  the  figures  with  Model  2  show  this  behavior  at  a  -  60* 

(e.g.  Figure  72).  The  effect  Is  stonger  at  60*,  since  this  Is  further  away 
from  the  sldelobes  from  the  nose-on  disk. 

The  0*  aspect  runs  all  show  a  strong  specular  scatterer  near 
nose-on.  This  scatterer  has  a  Stokes  response  of  s^/Sq  -  0,  s^/ sQ  ■  t 

s^/Sq  -  0  If  the  Illumination  Is  135*  linear  (such  a  flat  plate  target 
leaves  linear  polarization  unchanged).  All  the  figures  from  63  to  76 
demonstrate  this  specular  response  near  nose-on. 

The  width  of  the  nose-on  response  as  a  function  of  aspect, 

particularly  measured  using  $2/Sq  »  40  Indication  of  the  diameter  of 

the  flat  plate  nose  of  the  aircraft  models.  As  the  scale  size  of  the  models 
Increases  from  A  to  D  ,  the  azimuth  extent  over  which  s2/Sq  remains 
near  -1  decreases  as  the  Inverse  of  the  scale  size. 


I  -  AIRCRAFT  S«. 1.2.3  VS  hSPECT  4SLIN  POLZII  DATA  P19R02  07-JUI-82 

FIG.  63  Measured  Stokes  parameters  vs.  aspect  of  cylinder.  Scale  A;  ka  ~  5.24; 
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•  PCA5  -  A/C-1000  SO.  1 . 2. 3  US  ASPECT  POLES  -1.1  DATA  P1SD01  16-JUL-8E 
FIG.  65  Measured  Stokes  parameters  vs.  aspect  of  model  2,  Scale  A;  ka  “  5.24; 
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a  PCA3  -  AIRCRAFT  $0. 1.2.3  US  ASPECT  4SLIM  POLZN  DATA  P23R05  09-JUL-82 
FIG.  M  Measured  Stokes  parameters  vs.  MP«et  of  model  3,  Scale  A;  ka  *  8.14,; 
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•  PCA3  -  AIRCRAFT  St. 1.2. 3  VS  ASPECT  4SLIH  P0L2H  DAT APS IRAS  M-JUL-S8 

FIG.  08  Measured  Stokes  parameters  vs.  aspect  of  model  1,  Scale  B;  ka  •  8.14; 
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FIG.  70  Mi— jrtd  Stokos  paramctarv  ya.  aapoct  of  modal  3,  Scat*  B;  lea  -  12.8; 
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FIG.  82  Measured  Stokes  parameters  vs,  aspect  of  model  3.  Scale  B;  ka  ■  12.8; 
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FIG.  84  Measured  Stokes  parameters  vs.  aspect  of  model  2,  Scale  C;  ka  *  12.8; 
9  -  10°. 


At  10*  viewing  aspect,  the  measured  results  found  in  Figures  77 
through  88  show  the  nose-on  total  power  response  sQ  ,  in  most  cases,  to  be 
decreased  in  magnitude  as  it  is  expected  to  be.  For  example,  consider  the 
results  corresponding  to  scale  size  A  and  aircraft  Model  1  found  in  Figure 
77.  Although  the  10*  aspect  results  show  the  nose-on  total  power  sQ  to  be 
approximately  4  to  5  dB  below  that  of  the  0*  aspect  run,  Figure  64,  the 
azimuth  extent  around  nose-on  over  which  Sg/Sg  is  approximately  -1  is 
roughly  the  same. 

A  strong  scatterer  is  also  seen  in  the  measurement  results  of  total 
power  Sq  near  broadside  for  all  target  scales  and  models.  This  scatterer 
is  identified  as  single-bounce  specular  most  convincingly  in  the  measurement 
results  corresponding  to  cylinder  only  (aircraft  without  wings)  for  the  scale 
sizes  measured.  It  is  clearly  seen  in  Figures  63  and  67,  corresponding  to 
scale  size  A  and  B  cylinder  only,  respectively,  that  s^/Sg  is,  on 
average,  near  -1,  while  both  s^ sQ  and  $3/sQ  are,  on  average,  0 
near  broadside,  indicating  a  single-bounce  specular  scatterer.  Somewhat 
surprising  is  the  fact  that  the  measurements  of  the  aircraft  models, 
regardless  of  scale  >  wing  type,  do  not  show  this  same  polarization  response 
as  clearly  as  the  cylinder  only  measurements.  The  wing  causes  a  significant 
amount  of  fluctuation  in  the  SP,  particularly  in  s^/ Sq  and  s^/Sg  , 
which  obscures  the  expected  single-bounce  specular  scatterer  response  it 
broadside. 

In  general,  at  target  azimuths  that  do  not  correspond  to  a  clearly 
identifiable  scatterer,  such  as  the  nose,  wing,  or  broadside,  the  polarization 
signature  as  represented  by  the  SP  shows  a  great  deal  of  fluctuation.  The 
Interaction  between  scatterers  is  extremely  complicated  and  seems  unlikely  to 
be  of  use  for  target  classification  purposes.  A  method  suggested  elsewhere  in 
this  report.  Section  6,  considers  using  an  average  response  for  target 
classification  which  eliminates  some  of  these  fluctuations.  Preliminary 
predicted  response  results  show  this  technique  to  be  worth  further 
consideration. 
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The  measurement  results  reported  in  Figure  89  represent  the 
scattering  response  of  a  0*  aspect  or  plane- symmetric  target  for  two  incid  it 
polarizations,  135*  and  45*  linear.  Since  the  relative  phase  polarization 
scattering  matrix  for  this  target  has  zero  off-diagonal  elements,  it  is 
expected,  and  as  indeed  Is  the  case,  that.  sQ  should  be  identical  for  each 
transmit  polarization,  135*  and  45*  linear.  Similarly,  Sj/ sQ  should  be 
identical,  and  Is,  for  each  transmit  polarization.  It  can  be  shown,  on  the 
other  hand,  that  s2/Sg  for  135*  linear  transmit  polarization  should  be 
the  negative  of  Sg/Sg  for  45*  linear  transmit  polarization.  The 
measured  data  clearly  shows  this  to  be  the  case  for  Sg/Sg  and  for 
s3^s0  ’  should  also  exhibit  this  same  phenomenon.  A  deviation  from 

this  relationship  between  the  sets  of  SP  could  be  an  indication  of  the  lack  of 
system  calibration  or  asymmetry  of  the  aircraft  target  model. 

The  near  field  effects  discussed  in  Section  8.9  are  clearly  visible, 
for  example,  in  Figure  76,  in  which  the  broadside  RCS  for  the  Scale  D  Model  3 
target  is  smaller  than  the  nose-on  RCS.  Similar  effects  are  seen  for  the 
other  large  scale  targets  when  measured  at  the  shorter  (12  m)  range  (see  Table 
7).  These  effects  are  analyzed  in  Section  8.9.  It  is  felt  that  even  with 
degradation  in  the  data  due  to  near  field  effects,  the  polarization  signatures 
of  even  the  large  models  are  qualitatively  correct. 


4.6  COMPARISON  OF  MEASURED  AND  CALCULATED  RESULTS 

Using  the  same  modeling  program  that  ger.srated  the  polarization 
characteristics  of  Section  3,  the  Stokes  vectors  were  computed,  so  that 
computed  and  measured  data  could  be  compared  directly.  In  Figures  90,  91,  92, 
and  93,  the  computed  Stokes  vectors  are  presented  for  the  cylinder,  and  Models 
1,  2,  3,  respectively  at  size  scale  A.  These  results  are  to  be  compared  with 
the  measured  Stokes  vectors  for  these  targets  given  in  Figures  63  through  66. 
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FIG.  91  Predicted  Stokes  parameters  vs.  aspect  of  model  1,  Scale  A;  ka  s  5.24 


Consider  first  the  response  of  the  cylinder.  The  maxima  In  Sq 
(representing  the  total  power)  at  0*  and  90*  aspect  agree  within  2  dB  with 
measurements.  The  lobe  peaks  In  the  region  40*  to  70*,  where  the  response  Is 
lowest,  are  predicted  5  dC.  lower  than  measured.  This  discrepancy  Is  thought 
to  be  due  to  noise  and  clutter  and  is  supported  by  the  unexpected  envelope 
fluctuation  In  the  lobe  structure  In  the  measured  data.  The  lobe  spacing 
agrees  exactly.  The  normalized  Sj,  s?»  *3  curves  exhibit  similar 
rapid  fluctuations  for  the  measured  and  predicted  cases.  The  s^  curves 
(which  measure  the  cosine  of  the  relative  phase  between  the  horizontal  and 
vertical  components  of  the  scattering  matrix)  exhibit  similar  average 
structure  for  the  measured  and  predicted  cases.  The  details  of  the 
fluctuations  (relative  phase  with  respect  to  aspect)  of  the  predicted  s^, 

Sg,  S3  curves  differ  from  the  measured  cases.  This  Is  considered  to  be 
explained  by  the  fact  that  the  fluctuations  are  highly  sensitive  to  relative 
scatter  position  at  the  high  frequencies.  We  thus  consider  that  agreement  of 
the  overall  properties  of  the  responses  constitutes  good  agreement. 

Similar  observations  hold  for  the  comparisons  between  measured  and 
computed  Stokes  vectors  for  Models  1,  2,  and  3.  Note  that  the  effect  of  the 
wings  (at  30*  or  60*)  shows  up  particularly  in  s\  . 

The  polarization  characteristics  m,  2v,  2y  were  derived  from  the 
measured  Stokes  parameters.  These  results  were  also  compared  with 
predictions.  Figures  94,  95,  96,  and  97  compare  the  computed  and  measured 
characteristics  respectively  for  the  cylinder.  Models  1,  2,  3,  at  size  scale  A, 
at  e  -  0*.  The  left  columns  of  these  figures  are  the  predicted  m,  2v,  and 
2y;  these  duplicate  the  results  already  given  in  Section  3.  We  see  that  the 
plots  of  m  compare  well,  subject  to  the  same  comments  made  above  for  Sq. 

The  plots  of  2v  agree  In  their  overall  properties  (rate  of  fluctuation),  but 
disagree  completely  In  detail.  This  parameter  Is  most  sensitive  to 
uncertainties  In  the  exact  prediction  and  measurement  of  rapidly  fluctuating 
targets.  The  2y  curves  show  fair  agreement,  although  they  again  disagree  In 
the  exact  detail  of  the  fluctuations.  In  particular,  note  that  all  2y 
curves  are  fairly  constant  near  nose-on;  and  note  the  effect  of  the  wings  at 
30*  In  Figure  95  and  97,  and  at  60*  in  Figure  96. 

These  curves  also  Illustrate  the  need  for  averaging  out  the  rapid 
fluctuations  In  the  characteristics,  as  discussed  In  Section  3.4. 
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FIG.  94  Predicted  snd  measured  polarization  characteristics  of  cylinder.  Scale  A; 
ka  -  5.24;  9  -  0°. 
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FIG.  95  Predicted  and  measured  polarization  characteristics  of  modal  1,  Scale  A; 
ka  -  5.24;  9-0°. 
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SECTION  5 


EFFECTS  OF  NOISE,  CHAFF,  AND  JAMMING  ON  POLARIZATION  MJLIS 


In  this  section  the  effects  of  noise,  chaff,  and  jamming  on 
polarization  null  estimation  Is  Investigated  by  obtaining  a  bound  on  the  area 
of  the  concentration  ellipse  of  the  estimate  probability  distribution.  This 
bound  Is  used  to  evaluate  the  deterioration  In  accuracy  of  estimation  of 
co-polarizatlon  null  coordinates  of  a  simple  symmetric  target.  A  technique  to 
minimize  the  effect  of  a  jammer  Is  then  postulated.  In  which  receive 
polarizations  are  chosen  parallel  and  perpendicular  to  the  jammer,  and 
transmit  polarizations  are  chosen  in  orthogonal  pairs  at  six  antipodal  points 
on  the  Poincare  sphere.  The  gain  afforded  by  this  technique  is  evaluated  for 
six  different  jammer  polarizations,  and  is  shown  to  be  significant  for  some 
jammers,  but  a  loss  for  others.  The  loss  Is  attributed  to  Improper  choice  of 
transmit  polarizations. 


5.1  INTERFERENCE  EFFECTS 


The  return  signal  is  first  modeled  when  interference  consists  of 
receiver  noise,  chaff,  and  jamming.  This  leads  to  a  probability  distribution 
for  the  observed  detector  outputs  of  a  polarization-sensitive  radar.  From 
this,  a  bound  on  the  area  of  the  concentration  ellipse  of  any  unbiased 
estimate  of  the  co-polarization  null  coordinates  of  a  symmetric  target  Is 
obtained  by  evaluating  the  Fisher  information  matrix  for  this  distribution. 

The  bound  is  then  evaluated  for  backgrounds  of  receiver  noise  plus  chaft  and 
receiver  noise  plus  jamming. 

5.1.1  ^Return  Model 

The  complex  voltage  V  measured  by  a  coherent  polarization-sensitive 
radar  when  chaff,  jamming,  and  receiver  noise  are  present  can  be  represented 
as: 

V  -  rT[(S  ♦  C)t  +  x£]  ♦  n  ,  (44) 
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where 


£  and  r  are 'the  transmit  and  receive  polarization  vectors 
of  the  radar; 

^  is  the  jammer  polarization  vectov; 

S  Is  the  target  scattering  matrix 

C  is  the  coherent  sum  of  the  scattering  matrices  of  all  chaff 
scatterers; 

x  is  a  complex  jammer  waveform  sample; 
n  Is  a  complex  receiver  noise  sample  . 


For  simplicity,  we  will  restrict  attention  In  this  section  to 
symmetric  targets,  for  which  the  target  scattering  matrix  S  can  be  written, 
in  an  H,V  basis,  as 


s  * 


0  )  ** 
e~^2vtan2Y  / 


(45) 


where  v  and  y  are  the  coordinates  of  the  co-polarization  nulls  of  the 
o 

target,  and  m  is  proportional  to  the  signal  return  power.  The  phase 
angle  ot  is  determined  by  the  two-way  propagation  distance  to  the  target, 
which  is  a  very  large  number  of  wavelengths.  Model o  2*,  «t  can  therefore  be 
considered  to  be  a  random  phase,  uniformly  distributed  on  a  2w  Interval. 


The  chaff  voltage  Is  the  coherent  sum  of  contributions  from  all  chaff 
scatterers,  and  so  Is  the  overall  scattering  matrix: 


c  ■ 


(46) 


If  a  chaff  scatterer  Is  modeled  as  a  dipole  at  orientation 
horizontal,  then  the  chaff  scattering  matrix  from  the 
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where  p.  Is  proportional  to  the  Man  return  power  from  such  a  dipole  for 
H-H  transmit/receive,  averaged  over  the  distribution  of  *  and  where 
P^  Is  analogous  to  ,  and  Is  uniformly  distributed  on  (0,  2«).  In  a 
cloud  of  many  chaff  scatterers,  the  sum  of  the  returns  may  be  well 
approximated  as  a  complex  Gaussian  random  variable. 

The  janroer  Is  modeled  as  having  a  fixed  polarization  vector  ,  with 
a  common  random  modulation,  x  .  The  modulation  Is  modeled  as  a  complex 
Gaussian  random  process.  The  receiver  noise,  n  ,  Is  also  modeled  as  complex 
Gaussian. 


The  observable  quantity  Is  the  magnitude  or  envelope  of  V,  for  each 
of  several  pairs  of  transmit  and  receive  vectors: 

z  -  [  V  j  .  (43) 

Since,  by  the  above  models,  V  is  the  sum  of  a  fixed  but  unknown  signal  and  an 
additive  complex  Gaussian  random  variable,  z  has  a  Rician  distribution,  with 
parameters  A  and  o: 


P(z) 


exp[  “~y  (2^ 


♦  A2)]  I 


(49) 


where 


A  and  a 


are  given  by 


z  ,  signal  only, 

(50) 

-y 

z  ,  no  signal. 

(51) 

For  large  signal -to- interference  ratios  (large  A/o  ),  this  distribution  is 
well -approximated  as  a  Gaussian,  with  mean  value  A  and  variance  o2  . 


When  the  target  Is  symmetric,  as  in  Equation  (45),  A  is  given  by 

A  »  m  j  r^tj  +  r2t2  tan2?  e"^v|  (52) 
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where  the  transmit  and  receive  polarizations  are  written  (H-V  basis): 


(?) 

(?) 


The  variance,  o'  ,  Is  the  sum  of  three  components.  The  component 
due  to  chaff,  when  the  chaff  dipole  orientation  is  presumed  to  be 
uniform  on  (0,  2w),  can  be  shown  to  be 

?  1  (K1)  2|rltll?*  2|r2l2|2*  ,rl‘z  *  rZtll2i 

•c-  - 7  j< 

where  1/ 2 ( Pc )  is  the  variance  due  to  chaff  when  H-H  transmit/receive 
polarization  is  used  (i.s.,  r^  »  t^  ■  1,  r2  ”  ^  *  0).  comP°nen^ 
due  to  jamming  is 

°0  *  "7"  PJ  'rljl  +  r2j2  I2  (56 


where 


and  the  jammer  polarization  vector  is 


0 


The  component  due  to  receiver  noise  is  simply 

2  l  p 
°N  w  T  PN 


where 
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When  z  is  observed  sequentially  for  N  different  pairs  of  transmit 
and  receive  vectors,  and  when  sequential  observations  are  sufficiently 
separate  in  time  to  assure  independent  chaff,  jammer,  and  receiver  noise 
samples,  the  observations  are  independent  Gaussians,  whose  joint  distribution 
can  be  written 

1  1  0 

p<i>  -  n  277  e*p[  - 1?.  <zi  -  v2  ^  •  leu 

i=l  1  1 

where  and  correspond  to  the  i  transmit  and  receive 
polarization  vectors.  Target  polarization  null  classification  then  requires 
estimation  of  the  target  parameters  v  and  y  ,  given  these  N  observations. 
For  example,  one  such  estimation  technique  would  be  to  calculate  the  maximum 
likelihood  estimates,  which  are  the  values  of  v  and  y  which  maximize 
Equation  (61). 


5.1.2  Fisher  Information  Matrix  for  Estimating  v  and  y 

A  convenient  mechanism  for  investigating  the  accuracy  with  which  v 
and  y  may  be  estimated  is  the  Fisher  information  matrix,  J  ,  obtained  from 
the  probability  distribution  of  the  voltage  envelopes  z^  ,  given  in  Equation 
(61). 


The  Fisher  information  matrix  , 

,  hi  ji2 

J  a 

J21  hz 


is  defined  as 


(62) 


where 


(63 


(64) 


This  matrix  is  useful  in  the  following  way:  Any  estimates  v  and  y  will  be 
random  variables  themselves,  since  they  are  functions  of  the  random 
observations  z  .  If  the  mean  value  of  v  and  y  are  the  true  values  of  v 
and  y  ,  then  the  estimates  are  said  to  be  unbiased.  From  trial  to  trial, 
the  estimate  values  will  be  randomly  scattered  about  their  mean  values 
according  to  some  probability  distribution.  The  scatter  can  be  described  by 
the  concentration  ellipse,  which  is  defined  by 


T 

e 
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(66) 


where 


(67) 


is  the  estimate  error  vector,  A 

c 

Ac  =  E[  £T  £  ] 


is  the  estimate  error  covariance  matrix; 

,  (68) 


2 

and  c  is  a  constant.  If  the  error  distribution  is  joint  Gaussian,  the 

probability  P  that  the  error  vector  is  within  the  concentration  ellipse  Is  a 
2 

function  of  c  : 

.2 

P  -  1  -  exp(  -  )  .  (69) 

The  Fisher  information  matrix  is  useful  because,  for  an£  unbiased 
estimate  of  v  and  y  ,  the  concentration  ellipse  defined  by  Equation  (66) 
lies  either  outside  of  or  on  the  bound  ellipse  defined  by 


c T  J  c  -  c2  (70) 

where  J  is  the  Fisher  information  matrix  (see,  for  example,  page  81  of 
reference  15).  In  particular.  If  we  use  as  a  measure  of  estimation  accuracy 
the  area  of  the  concentration  ellipse  of  the  estimates,  this  area  can  be  no 
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smaller  than  the  area  of  -the  ellipse  defined  by  Equation  (70).  This  Is  the 
approach  taken  here.  Equation  (70)  will  be  evaluated  for  c2  «  2  (P  -  1  -  e"1) 
and  the  area  of  the  concentration  ellipse  Investigated  and  compared  for  various 
types  of  Interference. 


Application  of  Equations  (63),  (64),  and  (65)  to  Equation  (61),  when  thi 
target  Is  symmetric,  results  In  a  Fisher  Information  matrix  of  the  form: 
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In  Equations  (73)  and  (74),  the  1-dependence  of  r^,  r2,  t^,  and  t2 
has  been  suppressed  for  notatlonal  simplicity. 


The  area  of  the  concentration  ellipse  of  Equation  (70),  for  c^  =  2, 
can  be  shown  to  be  given  by 


ii  11J09  **  j 
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Thus,  small  concentration  ellipse  areas  are  produced  by  large  diagonal 
elements  In  J  (j^  and  j22)  and  small  off-diagonal  elements  (jj2). 

Large  diagonal  elements  are  produced  by  large  slgnal-to-lnterference  ratios 
(m1 la])  ,  and  by  large  values  of  the  partial  derivatives  of  Equations 
(73)  and  (74).  Small  off-diagonal  terms  are  obtained  by  selecting,  where 
possible,  transmit  and  receive  polarization  vectors  In  complementary  pairs,  so 
that  the  off-diagonal  terms  have  equal  magnitude  but  opposite  sign,  and  thus 
cancel  when  the  J^'s  are  summed  over  all  polarizations  i  .  This 
technique  will  be  Illustrated  in  Section  5.2. 


5.1.3  Effects  of  Chaff  and  Jamming  Using  Nominal  Polarization  Sets 


In  the  presence  of  receiver  noise,  a  nominal  set  of  N  «  12  pairs  of 
transmit  and  receive  polarizations  has  been  selected.  They  are  sufficient  to 
allow  estimation  of  target  polarization  nulls  for  an  arbitrary  target.  In 
this  section,  the  effect  of  chaff  and  jamming  will  be  evaluated  on  the  area  of 
the  concentration  ellipse  of  estimates  of  y  and  v  ,  when  the  target  is 
symmetric  and  the  nominal  set  of  twelve  pairs  of  polarizations  are  used.  The 
equations  developed  in  the  previous  section  will  be  used.  . 


The  nominal  set  of  polarizations  are  the  twelve  pairwise  combinations 
of  two  transmit  polarizations  (45*  and  -45*  linear) 


i'C)  an‘  (-!) 


and  (  I  (76) 

and  six  receive  polarizations  (45*  and  -45*  linear,  H  and  V,  and  RC  and  LC): 


The  target  is  a  synmetrlc  target,  with 

y  -  45*  , 


(77) 


v  ■ 


0 


(78) 


p 

i.e.,  a  sphere  or  flat  plate.  The  signal-to-receiver  noise  ratio  (m  /PN) 
is  +20  dB.  Figure  98  is  a  plot  of  the  concentration  ellipse  area  of  Equation 
(70)  as  a  function  of  chaff-to-recelver-noise  ratio,  when  the  Interference 
background  is  receiver  noise  plus  chaff.  The  effect  of  chaff  on  the  ellipse 
area  is  essentially  the  same  as  an  equal  amount  of  receiver  noise.  For 
example,  at  a  chaff-to-receiver-noise  ratio  of  0  dB,  the  total  background 
power  is  twice  that  of  receiver  noise  only,  and  the  ellipse  area  is  also  twice 
that  for  receiver  noise.  This  property  is  not  surprising  in  view  of  the  fact 
that  the  chaff  model  used  contains  no  preferred  dipole  orientation. 

When  the  interference  is  a  jamming,  however,  the  result  is  a  function 
of  the  jammer  polarization  vector.  Figure  99  shows  seven  jammer  polarizations 
chosen  for  calculations.  Figure  100  shows  the  variation  of  concentration 
ellipse  area  with  jammer-to-noise  ratio,  when  the  interference  background  is 
receiver  noise  plus  a  jammer.  The  target  is  again  a  sphere  or  flat  plate,  and 
the  signal-to-receiver-noise  ratio  is  +20  dB.  In  this  case,  each  jammer 
polarization  has  a  somewhat  different  effect  on  the  estimate  accuracies. 

5.2  ECCM  TECHNIQUES 
• 

In  this  section,  the  properties  of  the  Fisher  information  matrix 
which  lead  to  small  concentration  ellipse  areas,  which  were  identified  in 
Section  5.1.2,  will  be  exploited  to  obtain  an  ECCM  technique  for  countering 
the  effects  of  jamming  when  estimating  y  and  v  of  a  symmetric  target. 
First,  a  specific  example  will  be  discussed  to  motivate  the  postulation  of  a 
counter-jamming  technique  for  a  more  general  situation.  The  technique  will 
then  be  evaluated  in  the  more  general  situation. 

5.2.1  Example  Case 

When  the  target  polarization  parameters  are  estimated  from  a  sequence 
of  N  independent  voltage  envelope  measurements,  i  ,  the  Fisher  information 
matrix,  0  ,  as  we  have  seen,  is  the  sum  of  N  terms,  one  for  each 
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FIG.  100  Polarization  null  estimation  accuracy  in  jamming. 
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measurement .  Each  term  Is  dependent  on  the  slgnal-to-lnterference  ratios  and 
the  partial  derivatives  of  the  mean  return,  for  the  associated  transmit  and 
receive  polarizations.  The  concentration  ellipse  area  can  be  made  small  If 
the  diagonal  terms  of  the  sum  J  are  large,  and  the  off-diagonal  terms  are 
small.  If  at  least  one  of  the  N  terms  in  the  sum  has  a  large  diagonal  term 
for  each  parameter,  and  If  the  terms  can  be  grouped  In  pairs  such  that  the 
off-diagonal  terms  are  of  equal  magnitude  but  opposite  sign,  the  desired  net 
effect  will  be  achieved.  Therefore,  the  key  to  good  performance  in  an  ECCM 
environment  is  to  examine  the  J^'s  for  various  transmit  and  receive 
polarizations,  and  select  polarizations  to  produce  the  desired  effect  on  the 
sum. 


As  an  illustration,  consider  the  symmetric  target  again,  with 
y  *  45* ,  v  »  0*,  in  a  background  of  receiver  noise  plus  jamming.  Table  8 
is  a  tabulation  of  the  terms  of  the  B.  matrix  (Equation  (72))  for  each  of 
36  pairings  of  transmit  and  receive  polarizations  chosen  from  the  set: 

45*  linear,  -45*  linear,  H,  V,  I.C,  and  RC. 


Suppose  we  were  to  choose  from  these  polarizations  the  set 
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If  o?  .  O  ,  the  desired  cancellation  of  off-diagonal  terms  occurs. 

1  c  2  2 

However,  since  r ^  ^  r^  ,  we  may  have  a\  t  °2  *  ^  we  a<^  *° 

the  set,  tg  ■  -45*  linear,  rj  ■  LC  ,  t^  «  -45*  linear,  r^  ■  RC  , 

then  two  more  terms  are  added: 
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is  diagonal,  and  hence  so  is 
'4  0 


(83) 


An  alternate  choice  of  r  and  t.  pairs  would  be  to  interchange  the 
transmit  and  receive  polarizations  in  the  above  selection,  so  that  transmit 
polarizations  are  LC  and  RC  ,  and  receive  polarizations  are  *45*  linear. 

are  unchanged  by  the  interchange  of  t^  and 
Thus,  we  get 

4  0 
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The  partial  derivatives  of 
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£  ,  but  the  values  are. 
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The  better  technique  depends  on  the  jammer  polarization, 
jammer  polarization  is  45* 
linear. 


If  the 
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then 
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For  the  ±45*  linear  polarization  reception  case,  we  get 
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For  all  0/N  ratios,  we  have 
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so  that  the  latter  choice  of  polarization  sets  Is  better. 
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When  the  jammer  Is  LC  or  RC,  the  opposite  conclusion  Is  true.  This 
suggests  the  following  counter  jamming  scheme: 

1.  Choose  receive  polarization  vectors  parallel  and  orthogonal  to  the 
jammer  polarization  vector. 

2.  For  each  of  these  two  receive  polarizations,  transmit  on  pairs  of 
polarization  vectors  chosen  so  that  the  off-diagonal  terms  of  the  pair 
of  corresponding  6^  matrices  are  of  equal  magnitude  but  opposite 
sign,  and  thus  cancel  when  the  terms  are  summed. 


Step  1  requires  knowledge  of  the  jammer  polarization,  which  should  be 
feasible  as  long  as  the  jammer  polarization  does  not  change  rapidly.  For  Step 
2  It  Is  not  clear  how  to  choose  the  appropriate  transmit  polarizations.  The 
above  example  suggests  that  If  transmit  pairs  are  orthogonal,  the  dnclred 
property  will  be  achieved  for  the  particular  target  and  jammer.  The  correct 
transmit  polarization  set  choice  is  not  fully  understood  for  a  general  target 
and  jammer  combination. 

5.2.2  Calculation  of  ECCM  Effectiveness  for  Six  Jammer  Polarizations 

To  test  the  above  conjecture,  the  concentration  ellipse  area  was 
calculated  when  measurements  are  made  with  twelve  pairs  of  transmit  and  receive 
polarizations,  for  six  janrner  types.  In  all  cases,  the  transmit/ receive  pairs 
are  obtained  by  combining  six  different  transmits  with  two  different  receive 
polarizations.  The  six  transmit  polarization  vectors  are  the  six  antipodal 
polarization  vectors  used  on  receive  in  the  nominal  set:  ±45*  linear,  H,  V, 
LC,  and  RC.  The  two  receive  polarizations  are  determined  by  the  jammer 
polarization:  one  is  parallel  to  the  jammer,  and  the  other  is  orthogonal.  The 
target  is  a  symmetrical  target,  y  -  45* ,  \>  =  0*.  The 

signal-to-receiver-noise  ratio  In  all  cases  is  +20  oB.  Figure  101  plots  the 
ellipse  area  versus  jammer-to-noise 

ratio  for  all  six  cases.  Comparison  to  Figure  99  shows  the  gain  (or  loss) 
achieved  by  the  technique.  Curves  1,  2,  4,  and  6  show  a  significant  gain  at 
high  jammer-to-noise  ratio.  These  curves  correspond  to  a  LC  jammer  and  three 
different  linearly-polarized  jammers.  At  +20  dB  0/N,  the  gain  ranges  from  3  dB 
for  LC  to  13  dB  for  30*  linear.  For  the  two  cases  of  elliptical  jammer 
polarizations,  however,  there  is  a  net  loss  of  about  4  dB.  This  loss  is 
believed  to  have  occurred  because  of  improper  selection  of  transmit 
polarizations  for  the  jammer-driven  receive  polarizations,  and  for  the 
particular  target.  The  desired  cancellation  of  the  off-diagonal  terms  probably 
did  not  occur.  More  work  needs  to  be  done  to  determine  how  the  transmit  set 
should  be  chosen.  However,  the  correct  choice  may  be  target-dependent,  thereby 
defeating  the  utility  of  the  technique  for  classification  of  an  unknown  target. 
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SECTION  6 

SUMMARY  AND  CONCLUSIONS 


We  have  given  the  theoretical  background  of  polarization  null  and 
maximum  properties  of  simple  and  compound  targets  as  functions  of  aspect  and 
frequency.  A  measurement  system  was  described  which  obtains  Stokes  vectors  of 
tne  scattered  radiation,  from  which  the  polarization  null  and  maximum 
properties  are  derived.  Computed  and  measured  results  are  given  for  three 
simplified  aircraft  models  at  several  size  scales.  The  effect  of  noise, 
clutter,  and  jamming  on  polarization  characteristics,  and  the  use  of 
polarization  agility  to  reduce  the  effect  of  jamming  was  Investigated. 

The  object  of  this  study  w^s  to  assess  the  efficacy  of  polarization 
null  and  maximum  concepts  for  aircraft  identification  at  practical,  single 
radar  frequencies.  It  was  found  that  simple  targets  have  highly  distinctive 
polarization  characteristics,  but  that  the  dynamic  behavior  of  the  polarization 
nulls  and  maxima  of  compound  targets  Is  very  complex.  The  conclusion  of  this 
study  Is  that,  in  general,  single-frequency,  instantaneous  polarization 
characteristics  are  not  useful  for  discrimination  of  extended  targets  (see 
Section  3.2).  An  exception  could  be  certain  distinctive  targets  over  certain 
restricted  aspect  bins.  The  reason  for  the  failure  is  purely  due  to  the 
complexity  of  the  dynamic  null  and  max  loci.  It  remains  true  that  the 
polarization  characteristics  of  targets  are  highly  distinctive  and  contain  more 
information  about  a  target  than  does  data  taken  at  a  single  polarization. 

It  is  necessary  to  find  another  approach  to  utilize  the  polarization 
information.  Two  methods  are  suggested  here.  One  way  is  to  overcome  the  rapid 
fluctuations  in  polarization  data  by  averaging  over  frequency  and/or  aspect, 
and  is  discussed  in  a  preliminary  fashion  in  Section  3.4  The  result  is  that 
different  aircraft  models  are  readily  distinguishable  over  certain  aspect 
bins.  A  more  complex  but  promising  approach  is  to  subresolve  the  target,  that 
is  utilizing  frequency  and/or  aspect  diversity.  (It  should  be  noted  that  the 
nulls  of  simple  scattering  centers  are  highly  characteristic  and  do  not  change 
rapidly  with  frequency  or  aspect.)  Full  analysis  of  these  two  methods  are 
beyond  the  scope  of  this  report  and  are  suggested  for  future  Investigation. 


The  restriction  of  this  study  to  practical  frequencies  and  aircraft 
targets  resulted  in  compound  targets  which  are  many  wavelengths  long.  We  may 
expect  quite  different  behavior  of  polarization  nulls  and  maxima  at  much  lower 
frequencies  or  for  much  smaller  targets,  that  is  in  the  resonance  region  where 
the  largest  target  dimension  is  of  the  order  of  a  wavelength.  Such  targets 
are  expected  to  show  slow  variation  of  null  characteristics  with  changes  in 
frequency  or  aspect.  (This  could  have  application  to  discrimination  against 
small  decoys.)  To  analyze  complex  targets  in  the  resonance  region  requires 
(computationally  intensive)  target  modeling  by  means  of  integral  equations, 
since  scattering  center  concepts  do  not  apply.  Such  methods  are  available, 
e.g.  [6,13]  and  could  be  pursued. 

In  the  presence  of  jamming,  the  accuracy  of  the  measurement  of 
polarization  characteristics  will  deteriorate.  By  proper  chDice  of 
measurement  polarizations,  the  effect  of  a  jammer  may  be  overcome  in  many 
cases.  It  was  found  that  for  some  jammer  polarizations  a  simple  prescription 
for  measurement  polarizations  gives  good  results;  while  at  other 
polarizations,  the  search  for  optimum  measurement  polarizations  is  complex. 
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SECTION  8 


APPENDICES 


8.1  UNITARY  TRANSFORMATION  PROPERTIES 

Let  b  ■  Qa  ,  where  Q  Is  2  x  2.  (93) 

#  T  it 

Then  to  preserve  the  Inner  product,  (aj.ag)  -  aj8g  -  aja^ 
we  must  have  bjbg  -  ajQ'Q  a2  -  ajag  .  This 
Is  satisfied  If 

QTQ*  -  I.  or  Q'1  -  Q*T;  (94) 

that  Is,  Q  Is  unitary.  In  particular,  Q  will  preserve  an  orthonormal 
basis.  Writing  out  Equation  (94)  In  full.  It  Is  easily  found  that  the  rows  of 
q  are  orthonormal  and  that  the  magnitude  of  det(Q)  -  1.  We  will  further 
restrict  the  phase  of  the  determinant  to  be  zero. 


If  b  and  a  are  any  unit  polarization  vectors,  then  Equation  (93)  can  be 

*  * 

solved  for  Q  In  the  above  form,  namely  q^  -  b^  ♦  bga2  and 
q2  »  -a^b2  +  bj«2.  Hence,  there  Is  a  Q  for  every  polarization 
basis  transformation.  We  have  that 

1  2  2 

l<*ll  +  1*2  I  “  1 

Det(Q)  .  l 

and  that  the  product  of  two  such  matrices  Is  a  matrix  of  the  same  form. 


It  is  not,  in  general,  true  that  a  Q  of  form  (Equation  (95))  can  be 
found  which  transforms  any  orthogonal  pair  of  basis  vectors  into  any  other 
orthogonal  pair.  For  that  we  would  have  to  maintain  a  more  general  form  with 
det(Q)  ■  e^.  The  reason  is  that  a  pair  of  basis  vectors  is  still 
orthogonal  if  one  of  them  is  multiplied  by  a  constant  phase.  Suppose 
are  orthogonal  with 

f  *  V 


and  A- 


*2  e 


* 

-a,  e 


V1 


and 


that  Bj,B2  are 


orthogonal  with 


1  ■! 

V1 


\ 


v 


* 

b2  e 


and  B2  - 


„*  ej‘b 

-bi e 


for  the  more  general  Qq  ,  obtaining 

r 

% a 


) 


y-q2e 


*  J(V^a}  *  jMaJ 

qle 


Then  we  can  solve 


(96) 


with 


det 


J(‘b-  ♦,) 


In  our  case,  we  will  always  deal  with  relative  phase  polarizations,  so  that  we 
can  use  the  more  restricted  form  (Equation  (95)). 
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8.2  THE  ELEMENTARY  TRANSFORMATIONS  OF  POLARIZATION  AND  STOKES  VECTORS 


/ bi^  ( 

Let  5  I-  Ql  I,  with  respective  Stokes  vectors 

W  Vi/ 


7bii2+ib2r 

k i2  ♦  i  «2  r 

i  bi  i2- 1  b2 12 

12  -  1*2  I2 

2Re(bjb2*) 

and  s.  - 

2Re(aja2*) 

-2Im(bjb2*) 

-2Im(aja2*) 

We  want  to  find  V  such  that  sb  •  Vsa  ,  in  terms  of  qj,  q2  of  Equation  (95). 
This  can  be  done  directly  by  writing  out  the  matrix  multiplication  and 
Identifying  terms.  The  result  Is 

fl  0  0  0  1 


0  lql  |2”  lq2 12  2Re(Vz*> 

0  -ZRefqjqg)  Re(qj2-  q22) 

0  2Im(qjq2)  -Im^2-  q22) 


Zlmfq^*) 
Imtqj2^  q22) 

Re(q]L2+  q22) 


(97) 


We  now  note  that  rotations  about  the  axes  Sj,  s2,  s3,  respectively,  (refer  to 
Figure  2  In  Section  2.1),  have  the  matrices 
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These  correspond  respectively  (as  can  be  seen  from  Equation  (97)  to 


The  decomposition  V  >  \^  Vv  corresponds  to 

The  elementary  V  transformations  applied  successively  are  Euler  rotations 
(about  the  transformed  axes).  The  Indicated  angle  restrictions  ensure 
uniqueness,  for  a  fixed  order  f,  T,  v  . 
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A  wore  elegant  method  of  relating  the  transformations  of  polarltatlon 
and  Stokes  vectors  and  of  relating  the  scattering  matrix  to  the  Hueller 
matrix.  Is  by  means  of  the  Kronecker  product.  This  may  be  found  In  [11].  It 
has  the  advantage  of  providing  a  formalism  for  transforming  from  one 
representation  to  another,  from  which  properties  can  be  rigorously  derived. 


•ft)" 


One  form  of  the  Kronecker  product  of  two  complex  vectors  a 
is  defined  as 


(:;)• 


In  particular,  for 


*2bl 

*2b2 


(*)■ 


E®E*. 


(100) 


L^yf  J 

We  can  then  define  the  Stokes  vector  corresponding  to  E  (In  accordance  with 

Equation  (4))  as  _  . 

s  -  T  E  ©  E  ,  (1Q1) 


where 


1  0 
1  0 
0  1 
0  j 
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for  «  right-hand  coordinate  system.  (For  the  left-hand  coordinate  system, 
which  also  appears  In  our  scattering  formulation,  we  must  use  -  T  . 
discussion  In  Appendix  8.3. 

We  also  need  the  Identity  (A  Q  B)(a  Q  o)  -  Aa  Q  Bb 
Kronecker  product  of  two  matrices  Is  defined  as 

l 

AiiCB]  J  A„[B] 


See 


where  the 


I  12l 


A  B 

\  "?ltDJ 

I  J 

Let  the  scattered  polarization  vector  be  given  by 

E$  -  AEt  . 

Then  the  scattered  Stokes  vector  (defined  In  the  left-hand  system)  Is 


♦  it 

ss  -  T  Es  ®  Es 

-  T*(A  Et)  ®  (A  Et)* 

-  T*(A  ®  A*)(Et  ®  E*) 

S$.  Hst 

where 

st  "  T  Et  ©  Et 

and 

M  «  T  A  ®  A  T  1 

(102) 

M  Is  the  Mueller  matrix  corresponding  to  A  .  It  Is  readily  verified  that 


T"1  -  1/2(T#T) 


For  a  coordinate  transformation  (no  change  of  handedness) 


we  derive  similarly  the  corresponding  Stokes  vector  transformation 


S  -  *  *.  •  (103) 

where  V  -  T  Q  ®  Q  T  *  . 

The  latter  equation  can  be  worked  out  to  give  Equation  (97).  Equation  (103) 
can  be  used  to  derive  MQ  of  Equation  (20)  In  terms  of  the  components  of 
SQ  of  Equation  (16). 

8.3  SYMMETRY  RELATIONS  FOR  SCATTERING  MATRIX  AND  TRANSMIT-RECEIVE  CONVENTIONS 
The  linear  scattering  matrix  is  written 


N  i  as  ‘"'i 

Sign  conventions  for  the  ••monostatic  convention"  are  that  V  -  V  ;  H  -  H 
for  receiver  and  transmitter  at  the  same  location.  Identifying  H,  V  with 
x,y  and  writing  £  for  the  direction  of  propagation,  we  see  (Figure  102a) 
that  v\  k’  form  a  left-handed  triplet  if  Hs,  Vs,  ks  are 
chosen  to  be  right-handed. 

Since  much  of  the  theory  of  polarization  phenomenology  treats  the 
symmetric  case  aHV  .  «yH  ,  It  Is  important  to  Identify  when  this  relation 
holds.  This  is  done  carefully  in  [5].  First  consider  when  aHy  ■  »  0. 

define  a  "bistatic  plane"  which  contains  fc1  and  £s  .  Define  the 
"polarization  reference"  planes  for  and  §s  Independently  as  the 

plane  containing  £  and  H  .  Then  a^y  ■  *y^  ■  0  when 

(1)  the  bistatic  plane  Is  a  symmetry  plane  of  the  target,  and 

(2)  the  bistatic  plane  Is  the  polarization  reference  plane  for  both 

and  Es  . 


SCATTER 


TAANMIT 


M  MONOtTATIC  CONVENTION: 


A,  AAA  A  A 

V-V'.hW,  PORk^-k* 


M  MULTIPLE  ECATTERtNO  IN  A  PLANE  bm 

FIG.  102  Propagation  and  acattsring  sign  convtntion*. 


The  aLove  are  shown  to  be  ‘rue  readily  by  reciprocity.  For  the  monostatic 
case,  given  a  plane  symmetric  target.  It  Is  always  possible  to  align  all  the 
refererce  planes  such  that  aHV  -  aVH  -  0  .  Based  on  the  above,  one  can 
readily  see  the  general  condition  for  a^y  -  a^  4  0  ,  namely  perform  the 
spatial  rotation  Qy  defined  In  Section  2,  which  rotates  both  polarization 
reference  planes  by  ¥  .  The  conditions  for  which  a^y  -  aVH  are 

(1)  the  bistatic  plane  is  a  symmetry  plane  of  the  target,  and 

(2)  transmitter  and  receive  polarization  reference  planes  are 
rotated  by  the  same  angle  with  respect  to  the  symmetry  plane. 

In  general  it  Is  true  that  a^y  -  aVH  If  and  only  If  there  exist  a 
polarization  for  which  the  transformed  A  Is  diagonal.  This  observation 
agrees  with  the  results  of  Section  2.3.  For  the  monostatic  case,  reciprocity 
directly  results  in  a^y  «  aVH  for  any  target. 

The  above  holds  for  linear  Isotropic  media  and  targets,  and  breaks 
down  when  reciprocity  does  not  hold. 

Because  of  the  reversal  of  the  basis  triplet,  one  has  to  be  careful 
when  multiplying  scattering  matrices  in  order  to  represent  multiple 
scattering.  The  simple  case  of  all  propagation  vectors  in  the  same  plane  is 
Illustrated  in  Figure  102b.  A  left-hand  triplet  is  Incident  upon  Aj  and 
scattered,  by  convention,  as  a  right-hand  triplet..  Before  applying  the  next 
scattering  matrix  A?  ,  it  Is  necessary  to  reverse  the  basis  back  to  a 
left-hand  one,  using  the  matrix  ^ j .  The  total  transformation  Is  then 

A  .  A2(;  >  (105) 
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For  example,  for  double-bounce  scattering  from  two  specular  surfaces  with 
A,  -  A,  -  I  ,  we  obtain  A  «  / 1  0  \  . 

V#-v 

The  use  of  the  monostatic  convention  Introduces  a  confusion  In  the 
definition  of  s.  In  the  Stokes  vector  (Equation  (4)  In  Section  2.1).  For 
the  right-handed  system  (H,  V,  k)  we  defined  s3  •  -2Im(sHey).  But 
then  In  a  left-handed  system  (H,  V,  -k),  we  must  take  s,  -  +2Im(eue!!). 

\  3  n  V 

j  |  and  s3  •  1  In  a 

right-handed  system  results  In  a  specularly  back-scattered  field  given  by  ( j 
In  a  left-handed  system  with  s3  >  -1  (RC).  On  the  other  hand,  45*  linear' 
polarization  t  -  f  j  J  with  s2  -  1  ,  upon  back-scattering  results  In  (  j  J 
also  with  s2  «  1. 

8.4  DIAGONALIZATION  OF  SCATTERING  MATRIX  AND  EIGENVALUES 


The  object  Is  to  diagonalize  the  scattering  matrix  with  a  unitary 
transformation.  This  will  be  achieved  by  solving  the  eigenvalue  problem 

Ae  -  xe*  (106) 

First  we  note  that  If  x,e  Is  an  eigenvector,  eigenvalue  pair,  then  so  are 
\  •  xe^,  e  -  e  e^.  See  this  by  multiplying  through  by  e^  : 

Ae  e’**  -  (xe^2^)(e*e-^); 

★ 

hence  Ae  -  x  e  . 

We  can  therefore  require  the  eigenvalues  Xp  x2  to  be  positive  real. 

To  solve  Equation  (106),  premultiply  by  A*: 

He  -  (A*A)e  .  xA*e*  -  |x|2e  (107) 
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(113) 


Then  Q  is  unitary  (Q*T  »  Q"1  with  |det(Q)  |  -  1).  Equation  (106)  can 
then  be  written 

QTAQ  -  A 

The  phases  of  x,  are  arbitrary  so  far.  Let  us  choose  some 
phase  such  that  x^x2  -  det(A)  -  det(A)  -  e“^B  |det(A)  |  .  Then  we  must 
have  det(Q)  »  1,  which  makes  it  a  unitary  transformation  In  our  restricted 
form.  Further,  choose  m  - 1  x^  |  >  |x2  (  ■  s  and  take  the  phase  factor  out: 

s0-("  “)-«NT*0.  («•> 

There  are  several  Invariants  under  a  unitary  transformation  that  are 
immediately  apparent  from  this  derivation:  If  matrix  B  Is  related  to  A  by 
a  unitary  transformation: 


B  -  QTAQ  , 

then  It  has  the  same  eigenvalues  as  A  ,  which  satisfy  the  same  characteristic 
equation,  namely,  the  relations  (108)  and  (109).  That  Is,  the  norm  and 
determinant  are  invariant. 

There  are  two  degenerate  cases.  The  first  Is  det(A)  -  0.  Then 

2 

m  «  norm(A)  and  s  «  0  .  This  happens  when  l^a^ !  ■  I*i2l 
and  Ui?*h#22)  1$  real.  The  other  case  Is 

2 

norm  (A) 

Jdet(A)}- _ . 

2 

In  that  case,  m2  «  s2  •  Jdet(A)j.  This  happens  when  Ja^J  -  |a22  j 
and  a^2  *  0  or  when  •  *22  “  Th®  degenerate  cases  are 
discussed  In  Section  2.4. 
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8.5  EQUAL  AREA  PROJECTION  FOR  POLARIZATION  CHART 


The  Huynen  polarization  chart  [4,7]  Is  an  orthogonal  projection  of 
the  Poincare  sphere  with  circular  polarization  at  the  center.  T*?o  charts: 
one  for  left-sense  polarization  and  another  for  right-sense  represent  the 
entire  sphere.  This  chart  is  Ideal  for  the  purpose  of  plotting  null  loci, 
except  for  the  fact  that  with  an  orthogonal  projection  the  area  nearest  the 
outer  perimeter  is  compressed.  Since  in  the  present  study  we  want  to  examine 
the  separability  of  null  loci,  it  is  desirable  to  avoid  that  distortion.  Me 
therefore  use  a  polar  equal  area  projection,  so  that  the  separation  between 
two  loci  near  the  outer  perimeter  appears  visually  the  same  as  for  loci  near 
the  center. 

The  equal  area  polarization  chart  is  showr,  in  Figure  (103).  At  its 
center  is  the  axis  s3  and  LC  polarization.  The  perimeter  is  linear 
polarization  with  H  (axis  Sj)  at  the  top  and  V  at  the  bottom.  The 
rotation  corresponds  to  the  polarization  tilt  angle  y.  The  radius  from 
the  center  Is  related  to  the  ell ipti city  by 

r  -  yi  +  cos  2  t  -  ^1  -  cos  Zt  (115) 

(By  contrast,  in  the  Huynen  polarization  chart,  r  -  cos  2t  .)  Except  for  the 
distortion  in  the  radius  and  our  usage  here  of  placing  H  on  top,  the  chart 
is  the  same  as  the  Huynen  chart. 


8.6  MODELED  RESPONSE  OF  RIGHT-CIRCULAR  CYLINDER 

The  cylinders  used  In  this  report  all  have  length,  L,  much  greater 
than  the  radius,  a,  with  the  radius  In  the  upper  end  of  the  resonance  region 
and  much  higher.  Approximate  analytic  expressions  exist  for  this  case  which 
exhibit  the  proper  polarization  dependence.  The  response  curves  as  a  function 
aspect,  a,  are  divided  Into  three  sections:  near  nose-on  {«  »  0),  near 
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broadside  («  s  90),  and  -oblique.  We  will  give  the  responses  as  coefficients 
6,  normalized  such  that  the  RGS  Is 


o  ■  v 


*  I  ®l! 


(116) 


Near  nose-on,  we  have  simply  the  response  of  a  flat  disc.  Creeping 
waves  can  be  neglected  due  to  the  great  length  of  the  cylinder.  A  physical 
optics  approximation  can  be  made,  because  we  will  use  this  model  only  In  the 
main  lobe  of  the  return.  Then 


(jka)  cos  a 


2J, (x) 


(117) 


where 


-  2ka  sin  a 


cos  a 


(The  factor  j  ,  often  not  written  In  published  results.  Is  Included  to  give 
the  proper  90*  phase  shift.)  The  phase,  0,  Is  with  respect  to  the  cylinder 
center.)  This  expression  Is  used  over  0  <  x  <  3.83,  where  3.83  Is  the  first 
zero  of  J^x). 

At  oblique  angles  the  GTD  expressions  are  useful,  as  given  In  [9]. 

For  the  RCC,  for  the  monostatic  case,  the  expression  for  scattering  from 
cylinder  edge,  k.  Is 

Gk  ■  -f5,n »  Vrrc-b^  v  *  sra 

where  ^  <118> 

v  -  120*, 

the  top  sign  Is  for  horizontal  (TM)  polarization  (H  i  to  symmetry 
plane). 

For  aspect  a  measured  with  respect  to  the  cylinder  axis,  uk  In  the  above 
has  the  values 

■  a  ,  <>2  ■  -a  ,  -  180*  , 


the  edge  In  the  third  quadrant  Is  shadowed  and  does  not  contribute.  The  total 
cylinder  response  Is  then 


#2  -  2ka(-s1n  «  ♦  -jj^cos  a), 

#4  -  2ka(s1n  a  -  -j^cos  a). 

Note  that  (118)  only  has  a  (ka)~^z  dependence  on  frequency,  and  will  fall 
at  low  frequency.  Also,  typically  for  GTD,  the  expressions  fall  near  grazing 
angles. 


Near  broadside,  we  use  an  approximation  derived  from  the  exact 
Infinite  cylinder  result  [10].  This  Is  a  good  approximation,  because  the 
length,  l,  Is  much  larger  than  the  radius,  a  .  The  expression  Includes 
polarization  dependence  In  the  specular  part  of  the  response  and  a 
polarization  dependent  creeping  wave  contribution.  In  the  following,  three 
terms  must  be  retained  of  the  series  expansion,  which  Is  then  applicable  for 
ka  >  2.  For  horizontal  (TM)  polarization: 


whare 

(cretp)H  - 

#  - 


cr 


F  . 


,-(5.05*j2.M)(k«)l/3  ,J*cr> 

2  ka  sin  a  , 

k«(2  ♦  ,)  (  V0*8"  }  • 

s1n{ka-j-cos  •) 

- C--J  — 


(120) 


ka~  cos  « 


The  factor,  F,  represents  the  effect  of  looking  slightly  off-broadside;  the 
creep  phase,  i  ,  contains  an  adjustment  for  locking  off-broadside.  The 

cr 

creeping  wave  magnitude  Is  not  quite  correct  for  a  4  90*  ,  but  its 
contribution  is  small. 

For  vertical  (TE)  polarization: 

«*  -  F  -  la|(”t)?)-  I5(3k,»'1'2 . Merely*. 

where  '  / 

(creep),  .  J.ttUk.-1-'8  J*cr  . 

These  expressions  are  used  In  the  main  lobe  of  the  broadside 
specular,  namely,  0  <  (ka)(L/a)cos  ;  <  w. 

The  broadside  specular  expressions  go  smoothly  into  the  6TD 
expressions  [9],  but  not  If  polarization  dependence  Is  Included  In  the 
broadside  return.  Near  nose-on  the  6TD  expressions  do  not  at  all  converge  to 
the  specular  expressions.  To  obtain  a  smooth  response  versus  aspect,  the 
cylinder  model  computer  program  Interpolates  between  the  two  expressions  over 
a  small  Interval  about  the  transition  point. 

For  reference,  the  responses  of  a  I/O  ■  10  cylinder  are  given  In 
Figures  104  and  105  at  ka  -  5,24  and  8.15,  respectively. 
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reneral  features  ot  the  RCc  response  that  ere  pertinent  to 
polarisation  null  behavior  are:  for  W  the  leading  edge  scatterer  Is 
dominant;  the  other  two  edges  are  weak  becoming  larger  In  their  respective 
specular  directions,  for  HH  the  two  edges  at  opposite  ends  of  the  cylinder 
have  strengths  on  the  same  order  of  magnitude.  Consequently,  the  HH  nulls  are 
deeper  than  the  W  nulls,  because  the  Interference  Is  stronger.  At  broadside, 
the  polarization  dependence  of  the  specular  and  the  creeping  wave  contribution 
are  strongest  for  W  polarization.  Consequently,  the  VV  broadside  response 
oscllates  slightly  with  Increasing  ka,  and  Is  usually  somewhat  smaller  than 
the  HH  response.  At  higher  frequencies  (ka  >  G),  the  HH  and  W  responses  are 
almost  the  same. 

8.7  MODELED  RESPONSE  OF  FINS 

The  aircraft  wings  were  modeled  as  flat  plates  of  general 
quadrilateral  contour.  There  are  no  closed-form  expressions  which  give  the 
response  at  all  aspects.  Because  of  the  large  number  of  response  calculations 
that  are  to  be  made.  It  was  necessary  to  use  simple  expressions.  On  the  other 
hand.  It  was  desired  that  these  expressions  give  a  good  approximation  to  the 
polarization  dependent  response  to  assure  meaningful  Interpretation  of  the 
computed  polarization  characteristics.  The  expressions  used  were  basically 
those  of  Physical  Optics  (PO)  and  the  geometrical  theory  of  diffraction  (STD), 
but  modified  In  accordance  with  experimental  results  and  results  from  exact 
calculations  In  the  regions  were  the  PO  and  GTO  expressions  fall. 

In  summary,  6TD  was  used  for  Incidence  normal  to  an  edge.  This  gives 
a  response  dependent  on  polarization  and  elevation  with  respect  to  the  plane 
of  the  wing.  For  aspects  in  the  plane  of  the  wing,  a  (sin  x)/x  dependence 
was  imposed  (as  predicted  from  PO  Interpretation  of  edge  currents).  At  large 
elevations,  the  expressions  for  leading  and  trailing  wing  edges  combine  to 
form  the  simple  PO  expression  for  a  flat  plate.  Near  grazing  angle,  the 
simple  STD  results  fall  drastically  [12,  13]  for  the  trailing  edge  of  a  plate 
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and  are  incorrect.  for  the  trailing  edge  even  at  larger  elevations.  The  effect 
may  be  Interpreted  as  due  to  the  existence  of  traveling  eaves  on  the  plate. 
Using  Space  Time  Integral  Equation  (STIE)  results  [13].  a  single  expression 
was  formed  giving  the  elevation  and  polarization  dependent  correction  to  the 
GTD  expressions.  The  expression  was  based  on  educated  guesses  as  to 
functional  behavior  and  forced  to  fit  the  data  found  In  [12.  13]. 

A  summary  of  the  plate  model  follows,  given  In  terms  of  the 
coefficient  G  ,  where  the  RCS  Is  o  »  w|G|*  (refer  to  Figure  106).  In  the 
main  lobe  of  the  flat  plate  specular  return  (e  close  to  90*) 


G  - 

j  *  ( s,"x  x  > 

where 

A  - 

the  area  of  the  plate 

x  - 

kL  cos  e 

L  - 

the  dimension  projected  In  the 

observation  plane. 

(This  has  no  polarization  dependence.  It  applies  only  in  the  main  lobe,  since 
It  does  not  include  any  dependence  on  the  plate  contour.)  Out  of  the  main 
specular  lobe,  GTD  expressions  are  used  for  the  individual  edges.  For  the 
leading  edge: 

6  -  jjp(  *  IK  •  <1?2> 

* 

where  the  top  sign  Is  for  E  In  the  plane  of  the  plate, 
the  bottom  sign  for  the  perpendicular  polarization 
0  «  edge  length 
x  -  kD  s1n(«  -  an), 

(a  -  an)  Is  the  aspect  In  the  plane  of  the  plate  with 
respect  to  the  edge  normal. 
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FIG.  106  Flat  Wing  Modal,  Showing  fading  and  trailing  edge  responses 
traveling  wave  interaction. 


This  expression  has  been  shown  by  exact  (SHE)  calculations  to  be  very 
accurate  for  all  o  .  Not  so  accurate  is  the  expression  for  the  trailing  edge: 


6 


1 

cos  o 


*  1)( 


sin  x 
x 


(123) 


For  e  -  0  this  gives  a  large  value  for  perpendicular  polarization,  whereas 
actually  a  sharp  null  is  observed;  and  gives  zero  for  parallel  polarization, 
whereas  actually  a  finite  response  is  observed.  To  this  is  added  a  "traveling 
wave  contribution": 


G 


f±  2  cos2o 


-jkL(l+COSO)±j^-  1n 
e  (— x“> 


(124) 

where  L  and  ±  have  the  same  meaning  as  above,  the  factor  f  is  chosen  as 
2 

f+  *  e-kDsin  °^2  to  obtain  agreement  with  exact  calculation  for 
parallel  polarization  and  f  -  1  (perpendicular  polarization).  The  above 
gives  good  agreement  also  with  perpendicular  polarization  data,  but  still 
fails  at  q  «  0  ,  for  which  a  null  is  desired.  Hence  the  entire  (sum)  result 
for  perpendicular  polarization  is  multiplied  by 

9  .  (1  -  e-k0s,"2<"'2)  .  (125) 


The  factors  f±  and  g  are  not  based  on  theory,  but  appear  to  result  in 
plate  responses  which  agree  with  experiment  and  exact  computations  using 
STIE.  Figure  107  shows  the  computed  response  of  a  plate  versus  aspect  at  0* 
and  10*  elevation.  The  perpendicularly  polarized  return  Is  zero  at  0*  and 
small  at  10*. 
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8.8  MULTIPLE  SCATTERING  -  AIRCRAFT  MODELS 


For  wings  attached  to  a  cylinder,  the  wing  Is  reflected  In  the 
cylinder,  resulting  In  an  even-bounce  return.  This  may  be  Important  for 
polarization  dependent  modeling,  since  the  signatures  of  even-  and  odd-bounce 
scattering  are  different.  For  small  elevation  angles,  the  following 
approximate  expressions  were  used  (applied  only  to  leading  wing  edges): 


Deff  20m*  m1n(s1n  a,  cos  a)  , 
Dm  -  m1n(D,  Lc  cos  o  tan  a)  , 


(126) 


where  L  is  the  fuselage  length  available  for  mirroring; 

0  Is  the  edge  length. 

Then  the  response  Is  given  by  the  specular  GTD  leading  edge  expression  for  an 
effective  length,  °eff  ,  and  with  opposite  sign  for  perpendicular 
polarization.  If  the  wing  edge  Is  slanted,  the  main  lobe  of  the  even-bounce 
return  Is  bistatic  and  the  response  is  multiplied  by  a  (sin  x)/x  factor. 

There  are  also  triple  bounce,  partial  mirroring,  and  shadowing 
effects:  for  the  thin  fuselages  considered  here,  these  effects  nearly 
cancel.  Hence  the  total  odd-bounce  return  from  two  wings  was  modeled  as  a 
single-bounce  return  from  only  one  wing  and  its  full  image. 

There  are  also  double-bounce  returns  for  e  >  0  due  to  the  cylinder 
reflection  in  the  plate.  These  were  also  modeled,  but  turned  out  to  be  small. 

8.9  FRESNEL  REGION  RESPONSES 

It  was  unfortunately  not  possible  to  satisfy  the  far  field  conditions 
for  the  largest  of  the  targets.  We  attempt  to  estimate  the  severity  of  the 


Fresnel  region  effects  In  this  appendix.  The  exact  transmitted  field  equation 
for  a  uniform  source  J$(F')  -  -  x  over  a  square  region  2a  x  2a  In  the 
x',y'  plane,  for  transmission  In  direction  z  Is  given  by 


H(x,  y,  z) 


^)(x  x  G)  e‘JkR  dx'dy'  ,  (127) 


where 

In  the  Fresnel  region: 
and 

this  can  be  written 


R  .  I  r  -  r1 


(x  -  x')2+  (y  -  y')2  «  z2 


X  «  z  , 


H(x,y)  «  y  Hfar 


1 

Ta 


l 


e  “  to' 


.V 


zi'ix-x:.] 


m2 


xz 


|dz'  ! 

! 


where 


a-Jkz 


(128) 


Hfar  -  ^wZ (2a)2  Is  the  Ideal  farfleld. 


Writing  one  of  these  In  standard  form  for  y  -  0  : 

/ 

I 


>0 


j  e"^'  dy'  -y  xz/2 


X0  -j-y-t2  _ 

e  c  dt  »  Z^TzJT 


c(t0)  -  jS(tQ) 


where 


t  -  y'/^Tz72”  ,  tQ  «  a/y  xz/2 

C  and  S  are  the  Fresnel  Integrals 

„tf 


C(t( 


))  m  f  cos(-J-t2)  dt.  , 
0 


s(tn) 


•/ 


°s1n(-y-  t2)  dt  . 
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I  /2a  is  the  on-axis  Fresnel-zone  correction.  We  can  then  rewrite 
y0 

Equation  (128)  at  a  point  x  with  y  -  0  : 


wnere 


H(x) 

1- 


I. 


r0 


’far 


yn 


7T 


Ao 

7 T 


F(x) 


and 


F(x) 


(1-x/a) 


-t0(l+x/a) 


dt  * 


7xz7T" 


» 


using  t  »  (x*  -  x)/  /  Vz/2 

F(x)  is  the  off-axis  correction  factor,  normalized  such  that  F(Q)  « 
Writing  this  out 


(129) 


F(x) 


Am 


'{tctyi-  f)>  *  c(t0(i.  |))]  -  j[s(t0(i-  +  s(t0(i+  f))]  | 

(130)' 


This  factor  applies  to  both  receive  and  transmit  patterns,  for  scattering  from 
the  point  x  .  The  magnitude  of  the  two-way  amplitude  pattern,  |F|  ,  Is 

plotted  in  Figure  108  as  a  function  of  the  displacement  x  in  cm  off-axis  for 
the  antenna  (a  »  15  cm)  used  in  these  measurements,  at  the  distance 
z  *  12  m.  Indicated  on  the  plot  are  the  maximum  extents  (broadside  view)  of 
the  targets  at  sizes  A,  B,  and  C  .  it  is  seen  that  target  size  B 
extends  significantly  Into  the  edge  of  the  radar  beam.  Target  C  at 
broadside  is  poor. 


The  effect  of  the  Fresnel  pattern  Is  principally  the  apparent 
shortening  of  the  cross-range  dimensions  of  the  target.  The  result  of  this  Is 
to  decrease  the  response  and  to  Increase  the  spacing  of  the  lobes  near 
broadside.  The  effect  can  be  estimated  and  was  observed  for  target  size  B  . 
At  aspects  closer  to  nose-on,  the  target  cross-range  extent  is  smaller  and 
these  effects  do  not  matter  much. 
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FIG.  106  Two-way  Framd  pattern  factor  vs.  off-axis  position  at  Z  ■  12  m. 
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Th«  larger  targets  were  measured  at  a  range  of  21  m,  In  an  effort  to 
reduce  the  Fresnel  effects  (the  limitation  Mas  reduced  SNR).  The  factor 
|F(x)  I2  Is  shown  for  this  case  In  Figure  109.  It  Is  seen  that  target  size 
0  Is  still  poor,  and  that  E  (signifying  target  model  3  at  size  scale  0) 
extends  hopelessly  out  of  the  beam.  The  measured  responses  of  the  larger 
targets  arc  therefore  only  valid  for  aspects  of  less  than  about  45°. 

Although  some  details  of  the  measured  responses  are  Inaccurate,  such 
as  amplitude  and  lobe-spacing,  the  general  nature  of  the  polarization 
characteristics  will  still  be  exhibited.  This  justifies  the  use  In  this  study 
of  Fresnel-reglon  conditions. 

8.10  RELATION  BETWEEN  SYMMETRIC  PSM  AND  SIX  POWER  MEASUREMENTS 

For  a  plane-symmetric  target,  the  relative  phase  polarization 
scattering  matrix  (RPSM)  may  be  written  as 


I  “ill* 


I  *22 1 


(131) 


where  a^  -  (a^je  11  ,  a^  •  ja^e  22  are  the  diagonal  elements  of  the 
scattering  matrix  and  ay  •  4^  -  i^z  •  Nbw*  The  polarization  of  the  signal 
incident  on  such  a  target  is  135*  linear,  the  scattered  electric  field  1$: 


a,,|ejA 


l*22l 


(132) 


The  Stokes  parameters  describing  the  polarization  of  this  scattered  field  are: 

*0  "  l*lll2*  l*22|2 
*1  "  l*lll2-  I  *22! 2 
h  "  -*l  *iil  1*22 160* 

s3  "  -2I  all!  1*22  ls1n  Af 
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(133) 


and  It  It  not  difficult  to  then  show  that 


l«lll“  <*o+*l),2 
"  <*0  ‘  %\)fi 

AT  -  tan"l(t 3/*2)  .  (134) 

It  was  docldod  to  Measure  the  power  at  six  different  receiver 
polarizations,  V  and  H1  linear,  135*  and  45*  linear,  and  LH  and  RH  circular  and 
call  these  powers  P^,  Pg,  Pj,  P4,  P$,  and  Pfi,  respectively.  The  SP 
May  then  easily  be  related  to  these  powers  as  fellows: 


*o  -  V«y  pi>z  *  oy  p3>z  * 


V‘o  * 

<V  P1)/S( 

V*o  - 

oy  p  3>/s, 

Vso  ' 

<V  P5)/v 

»») 


By  coMblnlng  the  relationship  between  the  elements  of  the  RPSH  and  the  SP  and 
the  expressions  above,  the  elements  of  the  RPSM  May  be  expressed  In  terns  of 
the  six  measured  powers  for  each  target  aspect. 


8.11  RELATION  BETWEEN  PSK  ANO  TWELVE  POWER  MEASUREMENTS 


The  purpose  of  this  appendix  Is  to  describe  the  procedure  used  to 
determine  the  elements  of  the  relative  phase  polarization  scattering  matrix 
for  an  asymmetric  target  In  the  target  coordinate  system.  Remember  that  the 
radar  coordinate  system  Is  rotated  +45*  relative  to  the  target  coordinate 
system.  That  Is,  given  a  polarization  scattering  matrix  (PSM)  In  the  target 
coordinate  system  H,V 
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(136) 


*h,v  " 


the  PSM  In  the  rader  coordinate  system  Is 


'll 

*12 

12 

*22 

N'.V"  R  ^H,VR 


(137) 


where  R  1$  the  rotation  matrix 

COS  a  -Sin  a 

(138) 

sin  a  COS  a 

»  _ 

and  where  a  Is  specified  to  be  +45".  Performing  the  matrix  multiplication 


V.V  "  ~7~ 


♦  2.,;,  *  lK 

‘*11  *  *22 

P12 

■*n  * 

*1  r  2*12+  *22 

m 

J*12 

*>22_ 

(139) 


where  the  pj  j's  are  complex  and  defined  as  above. 


The  procedure  used  to  measure  the  elements  of  the  relative  phase 

scattering  matrix  (RPSM)  e  v  where  a^  -  ja^ie  » 

1 J  -  1,2  consists  of  four  parts  or  steps. 

First,  Illuminate  the  target  with  two  transmit  polarizations,  135* 
and  +45*  linear  (In  the  target  coordinate  system)  and  receive  the  power  In  six 
polarizations,  V'  and  H'  linear,  135*  and  +45*  linear  and  LH  and  RH  circular 
(In  the  radar  coordinate  system)  for  each  transmit  polarization  for  a  total  of 
twelve  power  measurements  at  each  target  aspect.  Call  these  powers  P^.,  P2V«, 

P3V'*  P4V*  ^sv'*  i,ld  P6V'»  the  135*  linear  or  V 

transmit  polarization  and  similarly  P1H,,  P2H.»  Pjj^t 

P6H'*  resPect<vely»  f°r  th<  +^5*  Hn«4r  or  H*  transmit  polarization. 
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Second ,  given  these  twelve  power  measurements  form  two  sets  of  Stokes 
parameters 


so,  •  pi/  *  <%„-  ”3/  *  (pu-  V* 

*U,S Ox  "  *P2x~  Plx*/s0x 

(140) 

*2x/s0x  *  *P4x“  P3x*/s0x 

s3x/s0x  “  *P6x“  P5x^s0x 

where  x  denotes  either  V  or  H\  Notice  that  all  twelve  measurements  are 

used  and  that  the  SP  In  each  set  s^x,  s2x,  and  s3x  are  guaranteed  to 
sum  square  to  one  as  tney  should  for  monochromatic  radiation. 

The  third  step  Is  to  relate  these  two  sets  of  SP  to  the  elements  of 
A^,  y 1  ,  up  to  a  relative  phase.  Assuming  that  p12  f  0  ,  then  It  Is 
possible  to  show  that 


1 P11 1  •  1 

^  *S0H'*  S1H'"‘ 

I  P22 1  “  } 

^  (s0v«-  siV'^ 

1  P12l 

f  (sov,+  *1V^ 

tan  af22  • 

_S3V,/S2 H1 

tan  afjj  ■ 

S3H'/S2H' 

(141) 


where  p^j 
1  -  1,2. 


|  Pf  j  |e  ,  l,j  ■  1,2  and  af^  ■  -  p32  , 


Once  these  p^'s  ,  the  elements  of  the  RPSM  In  the  radar 
coordinate  system,  have  been  found,  then  the  a^‘s  ,  the  elements  of  the 
RPSM  In  the  target  coordinate  system,  may  be  determined.  Specifically  It. can 
be  shown  that 
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and  then  that 
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Letting 


where 


*22  "  q22  * 

*q12  -  tan”1 
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(144) 
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The  two  relative  pheses  a#^  and  a#22  of  the  diagonal  elements  of  the 
scattering  matrix  In  the  target  coordinate  system  and  the  three  magnitudes 
|a^j  |,  la^l,  and  |a^2|  of  that  matrix  are  written  In  terns  of  the 
twelve  powers  actually  measured  by  the  radar.  Using  these  relative  phase 
scattering  matrix  elements,  the  polarization  nulls  and  maxima  may  be 
determined  as  described  In  Sections  2.10  and  3.3. 
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HISTORICAL  DEVELOPMENT  OF  RADAR  POLARIMETRY 


The  history  of  polarlmetrlc  radar  can  be  traced  from  the  phenomenon 
of  polar  electromagnetic  energy  dating  back  to  1669,  when  the  first  work  on 
the  subject  was  published  by  Erasmus  Barthollnus.  It  contained  his 
observations  that  objects  viewed  through  a  crystal  doubled  and  that  an 
Incident  light  ray  split  Into  ordinary  and  extraordinary  rays  (see  Kftnnen, 
1982). 


Barthollnus  was  followed  by  Christian  Huygens,  who  contributed  most 
significantly  to  the  field  of  optics  by  proposing  the  wave  nature  of  light  and 
by  discovering  polarized  light  (1677).  E.  Louis  Mallus  proved  Newton's 
suggestion  that  polarization  Is  an  Intrinsic  property  of  light  and  not 
something  added  by  a  crystal  (1808). 

The  next  significant  contribution  to  this  field  was  added  by 
Augustine  Fresnel  (1788-1827),  who  proposed  that  light  could  be  considered  as 
a  transverse  wave.  His  reflection  formulas  are  still  In  use  today  and  have 
been  rederived  using  electromagnetic  theory  by  Drude  (twentieth  century). 

The  last  of  the  early  pioneers  was  Sir  David  Brewster  who,  by 
extending  the  work  of  Mallus,  discovered  the  relationship  between  the 
polarizing  angle  and  the  relative  refractive  power  of  dielectric  materials 
(1816). 


The  work  of  all  of  these  scientists  Is  available  In  a  translated  form 
In  the  literature  by  William  Swindell  (1975). 

The  modern  application  of  radio  waves  began  with  the  postulation  of 
electromagnetic  theory  by  Michael  Faraday  In  1832,  and  James  Clerk  Maxwell's 
mathematical  formulation  of  the  behavior  of  these  fields  In  1873.  Soon  to 
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follow,  In  1886,  was  Heinrich  Herts's  demonstration  of  the  application  of  the 
theory  as  it  applies  to  radio  waves,  but  the  usage  was  mainly  restricted  to 
ship  range  detecting  in  darkness  and  fog  throughout  the  first  three  decades  of- 
the  twentieth  century.  A  good  paper  on  the  first  applications  of  radio  waves 
was  done  by  Marconi  (1922). 

The  use  of  RF  waves  for  aircraft  detection  and  the  design  of  the 
first  radar  was  accomplished  by  the  1930's  and  further  advanced  during  World 
War  II  (1939). 

The  first  extensive  polarimetric  wave  propagation  analyses  were 
carried  out  for  the  description  of  optical  transmission  devices  by  R.  C.  Jones 
(1941-1944)  under  the  guidance  of  Professor  Hans  Mueller  (MIT);  and  extensive 
studies  on  radar  polarization  were  initiated  in  1946  by  G.  Sinclair  (1948)  at 
the  Ohio  State  University,  ElectroScience  Laboratory. 

These  studies  were  further  pursued  by  Rumsey  (1949-1951)  and 
particularly  by  Kennaugh  (1948-1954). 

Basic  papers  appear  in  a  series  of  papers  in  the  Proc.  IRE,  May  1951, 
in  which  the  paper  by  Deschamps  (1951)  is  of  particular  use.  Utilizing  matrix 
algebra,  Kennaugh  Introduced  a  new  approach  to  radar  and  developed  the 
"Optimal  Target  Polarization"  concept  for  the  monostatic  relative  phase  case 
(Kennaugh,  1950-1954)  which  became  of  particular  Interest  to  meteorological 
radar  studies  (circular  polarization  clutter  cancellation).  There  were  a 
number  of  isolated  other  studies  going  on  as,  for  example,  GIT-Project  A235 
(July  1955)  for  the  purpose  of  using  polarization  to  distinguish  between 
targets  and  clutter,  which  is  well  reviewed  and  summarized  in  Root  (1980). 

The  decade  of  the  fifties  closed  without  any  real  recognition  for  the  need  of 
decoy  discrimination,  and  polarimetric  radar  still  remained  highly 
underdeveloped. 
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The  ballistic  missile  defense  program  gave  birth  to  the  first  real 
need  for  decoy  discrimination  and  target  classification,  and  in  May  1960,  J. 

R.  Huynen,  of  Lockheed  Aircraft  Corporation,  supported  by  AFCRL  Hanscom  Field, 
reported  the  first  attempts  of  decoy  discrimination  utilizing  Kennaugh's 
optimal  target  polarization  concept.  A  very  extensive  amount  of  measurements 
on  the  relative  phase  scattering  matrix  were  made,  and  the  most  important 
finding  of  Huynen  is  the  phenomenon  of  aspect  invariance  of  a  radar  target's 
characteristic  (co-polarization)  angle,  i.e.,  that  the  spherical  angle  between 
the  two  co-polarization  null  locations  on  the  Poincare  sphere  is  almost  aspect 
independent  for  simple  targets.  Unfortunately,  his  studies  were  terminated, 
and  we  refer  to  his  excellent  dissertation  which  summarizes  the  deep  insight 
Huynen  gained  into  radar  target  polarization  properties.  We  note  here  that 
Huynen  (1970,  1978)  developed  in  his  clutter  decomposition  theorem  a  very 
powerful  tool  for  cl  utter/ target  analysis  which  requires  extensive  further 
studies.  Again,  we  note  that  in  the  Russian  literature  (Kanareykhin,  et  al. 
1965,  1968;  Varshavchuk  and  Kobak,  1971;  Zhivotovski,  1973-1978)  the  potential 
applicability  to  target/ clutter  analysis  was  recognized  long  ago.  We  also 
note  here  that  Copeland  (1960),  under  the  guidance  of  Kennaugh,  developed  a 
practical  scheme  for  "classification  and  possible  identification  of  targets", 
based  cn  radar  polarlmetric  concepts  for  purely  symmetrical,  reciprocal 
targets . 


Instead  of  vigorously  pursuing  polarimetric  radar  technology,  it  was 
set  aside  in  favor  of  electromagnetically  incomplete  methods  of  wake  analysis, 
Doppler  studies,  and  the  like;  all  concerned  with  the  recognition  of  re-entry 
bodies.  The  main  bulk  of  polarimetric  studies  still  existing  during  the 
sixties  and  seventies  stem  from  programs  initiated  within  the  Army  Missile 
Research  and  Development  Labs  for  specific  missions  of  discriminating  between 
airborne  targets  and  decoys  (Root,  1980).  However,  we  must  note  that  several 
Important  papers  on  the  subject  of  polarimetric  radar  were  published  in  the 
Russian  literature,  and  we  refer  to  the  two  excellent  text  books  by 
Kanareykhin,  et  al.,  (1966,  1968),  which  should  be  published  in  English; 
immediately! 
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Systematic  polarlmetrlc  analyses  of  radar  backscatter  from 
hydrometeors  In  Isolation  and/or  distribution  supplemented  t  -  rigorous 
experimental  and  theoretical  studies  were  most  likely  first  initiated  by 
Glendon  C.  McCormick,  Laverne  E.  Allan,  and  Archibald  Hendry  at  NRC/ Ottawa  In 
conjunction  with  the  radar  meteorological  field  studies  at  Penhold,  Alberta, 
of  Brian  L.  Barge  and  Robert  G.  Humphries  of  the  Alberta  Research  Council. 
Based  upon  extensive  background  research  studies  on  available  polarization 
diversity  radar  techniques,  McCormick /All an/ Hendry  systematically  developed 
the  most  economical  polarization  radar  systems  applicable  to 
hydrometeorological  backscatter  analysis  using  circular  polarization  base 
vectors  [Allan  and  McCormick,  1978].  Particular  emphasis  Is  being  placed  on 
experimental,  theoretical,  as  well  as  computer-verlflcatlonal  radar 
backscatter  studies  of  oblate,  prolate  dielectric  spheroids  In  Isolation  and 
distribution.  For  these  specific  classes  of  targets  the  relative  phase 
backscatter ing  matrix  of  a  "reciprocal  body  of  revolution  (S^g  -  w  «** 
written  in  terms  of  a  system  of  coordinates  and  selected  parameters  useful  In 
describing  this  class  of  lossy  dielectric  spheroidal  scatterers.  In  addition. 
In  utilizing  properties  of  reduced  polarization  phase  extraction  of  a  circular 
polarization  base  vector  measurement  system,  most  likely  the  simplest 
polarlmetrlc  radar  measurement  system  for  hydrometeorological  scatter  analysis 
was  designed  [Allan  and  McCormick,  1978,  1980;  McCormick  and  Hendry,  1979, 
1981;  McCormick,  1981],  which  yet  represents  an  Incomplete  polarlmetrlc 
measurement  system. 

We  have  had  many  discussions  with  this  very  talented  group,  and  we 
were  disappointed  to  learn  that  they  did  not  consider  the  optimal  polarization 
null  theory  useful,  although  they  had  been  aware  of  the  studies  of  Kennaugh 
[1948-1952,  and  later],  Huynen  [1965,  but  net  his  Doctoral  thdtls  or  later 
publications],  and  Kanareykhln,  et  al.,  [1966,  1968].  McCormick  flatly 
opposed  until  recently,  the  fact  that  the  optimal  polarization  null  theory  Is 
of  any  usefulness  to  radar  hydrometeorological  Analyses.  However,  this 
personal  view  on  his  part  is  slowly  subsiding,  and  he  has  rewritten  his 
results  In  terms  of  the  optimal  polarization  null  concept,  which  dearly 
establishes  the  usefulness  of  the  concept  to  radar  meteorology  (private 
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communications,  1982):  In  truth,  our  elaborate  meticulously  deep  studies  show 
that  they  “disguised"  Kennaugh’s  optimal  polarization  null  theory  Into  their 
approach 1 


We  also  note  that  for  hostile  radar  target  discrimination, 
classification  and  Identification  in  meteorological  clutter,  the  polarimetric 
r^ar  measurement  system  proposed  and  designed  at  NRC-EMD  is  not  suitable; 
because,  for  man-made  targets,  the  co-pol  nulls  may  occur  close  to  linear  an 
well  as  circular  polarization;  whereas,  the  assumption  of  their  measurement 
system  squarely  lies  on  the  assumption  of  near  circular  co-pol  null  location 
on  the  polarization  sphere. 

In  our  opinion,  it  is  this  very  fact  of  not  recognizing  the  great 
potential  of  the  optimal  polarization  null  theory  that  device  technology 
related  to  radar  polarimetry  is  so  retarded.  Although  the  research  of 
McCormick,  et  al.,  has  found  full  recognition  among  meteorologists,  it  still 
needs  to  be  translated  into  the  formulation  competitive  with  utilizing  the 
"dynamic  polarization  fork  concept"  (Huynen,  1970;  Poelman,  1980;  Root,  1980; 
Bnerner,  1982).  We  emphasize  that  It  is  the  specific,  unique  property  of  the 
optimal  polarization  nulls  below  the  decorrelation  threshold  that  will  become 
so  decisive  in  target-versus-cl utter  discrimination  in  the  problem  considered 
here:  (i)  "The  electromagnetic  wave  propagation  assessment  in  the  ocean 
environment  for  the  marine  boundary  layer  characterization";  (11)  “The 
electromagnetic  wave  interrogation  assessment  in  ground-based  battlefield 
environment,  where  dust,  smog,  smoke,  muzzle  blast,  and 

precipitation/hydrometeorlc  clutter  screen  the  radar  vision",  and  last  but  not 
least,  ( 1 i i >  "The  electromagnetic  wave  Interrogation  with  general-motional  or 
stationary  ground-based  targets  observed  from  an  aerospace  platform 
incorporating  the  naval,  as  well  as  ground-base  clutter  aspects".  We 
emphasize  here  that  po lari met? 1c  radar  technology  has  advanced  to  a 
state-of-art  that  broadband  polarimetric  radar  systems  can  be  designed  which 
recover  complete  scattering  matrix  Information  within  time  frames  below  the 
decorrelation  (reshuffling)  threshold  times  of  the  vector  scattering  centers 
of  c 1 utter/ chaff » 
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During  his  dissertation  studies  at  the  Technical  University  Delft 
under  Professor  J.  P.  Schouten,  Huynen  strongly  Influenced  (1970  to  1975)  a 
young  technician  at  SHAPE-TC  and  engineering  student  at  THD,  Mr.  Andre  J. 
Poleman,  Mho  then  realized  the  dreams  of  Huynen  by  creating  a  complete 
polarlmetrlc  X-band  radar  (9.95  GHz)  with  real-time  optimal  polarization 
post-processing  facility  of  his  own.  This  Polarlmetrlc  Radar  System  of 
SHAPE-TC  Is  certainly  one  of  the  most  advanced  measurement  facilities  now 
available  for  analyzing  radar  target  acquisition,  tracking,  discrimination, 
classification  and  Imaging  problems  (Poelman,  1981)  In  most  adverse  clutter 
conditions.  Unfortunately,  that  system  Is  stationary,  and  we  would  require  a 
mobile  system  to  analyze  various  target  shapes  as  well  as  various  kinds  of 
precipitation,  land  and  sea  clutter,  and  such  systems  need  to  be  made 
available  Immediately  within  all  relevant  dm-to-mm  wave  propagation  windows  of 
Interest. 


In  recent  years,  there  has  been  a  renaissance  of  research  in  radar 
polarlmetry  from  both  a  theoretical  and  experimental  point  of  view  directed 
toward  the  determination  of  the  characteristic  properties  of  radar  targets 
utilizing  a  scatterer's  optimal  polarization  properties.  It  has  been 
established  experimentally  (Daley,  1878-1981,  Welsbrod  and  Morgan,  1979-1981), 
that  the  null  polarizations  can  be  used  in  order  to  discriminate  targets 
against  scattering  sea  clutter  which  Is  of  particular  Importance  to  naval 
operations.  This  meets  greatly  the  existing  tremendous  need  for  improved 
clutter  rejection  methods  In  order  to  detect  accurately  small  targets,  which 
Is  also  of  great  significance  to  avionics,  geophysical,  and  environmental 
remote  sensing.  It  has  been  shown  (Welsbrod,  et  al.,  1979),  In  the  case  of 
sea  clutter,  that  its  non-random  behavior  manifests  Itself  as  a  characteristic 
clustering  of  co-polsrlzed  and  cross-polarized  nulls  as  plotted  on  the 
Poincare  sphere.  This  clustering  was  noticeably  disturbed  with  the  presence 
of  a  target.  This  phenomenon  could  lead  to  low  false  alarm  rate  discriminants 
with  the  use  of  theoretical  models  extending  existing  clutter  statistics 
sensitive  to  the  changes  In  the  clustering  of  co-pol  and  cross-pol  nulls. 
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Again,  we  emphasize  that  Poelman  (1981)  clearly  verified  these  aspects  and  has 
without  any  doubt,  demonstrated  with  his  polar Imetrlc  radar  at  SHAPE-TC,  that 
target  versus  clutter  discrimination  Is  possible  under  the  most  adverse 
clutter  conditions  as  long  as  the  target  and  clutter  co-polarization  nulls  are 
separated  on  the  Poincare  sphere. 

Other  recent  studies  Include  polarization  effects  for  millimeter  wave 
propagation  In  rain  done  by  W.  S.  Vogel  at  the  University  of  Texas  In  1975. 

The  work  of  S.  I.  Metcalf  (1977),  formerly  of  GIT,  provides  a  survey  of 
research  relevant  to  the  measurement  of  atmospheric  parameters  with 
polarization  diversity  radar. 

A  study  of  volumetric  effects  in  depolarization  of  electromagnetic 
waves  scattered  from  a  rough  surface  was  done  by  Andrew  J.  Blanchard  (1977), 

In  which  he  verified  the  thesis  that  depolarization  of  backscatter  from  a 
finitely  conducting  media  Is  predominately  a  multiple  scatter  volumetric 
effect.  He  developed  a  mathematical  model  capable  of  describing  the 
volumetric  effect  in  the  depolarization  of  EM  waves  scattered  from  rough 
surfaces,  collected  data  In  the  microwave  and  optical  frequency  regions,  and 
compared  model  performance  with  depolarization  results  accumulated  In  the 
measurement  program.  Another  more  advanced  study  of  the  same  problem  was 
carried  out  by  Fung  and  Eom  (1980,  1981)  and  Is  being  further  extended. 
Incorporating  the  optimal  polarization  null  concept  (Fung  and  Eom,  1982), 
which  proved  to  be  very  successful.  We  note  here  that  Blanchard  and  Thelss 
have  been  able  to  substantiate  this  finding  in  more  detail  (Blanchard  and 
Thelss,  1982). 

Recently,  another  area  to  which  polarization  has  been  applied  Is  its 
Incorporation  Into  the  equations  of  transfer,  or  to  be  more  specific, 
equations  of  radiation  hydrodynamics.  The  work  on  the  subject  Is  summarized, 
for  example.  In  Ishlmaru  (1978),  Bom  and  Wolf  (1964),  Wolf  (1976),  Zubalry 
and  Wolf  (1977),  as  well  as  In  G.  C.  Pomraning  (1973). 
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As  far  as  the  analytical  theory  of  polarization  Is  concerned,  one  of 
the  latest  summaries  was  done  by  R.  M.  Azzan  and  N.  M.  Bashara  of  the 
University  of  Nebraska  In  1977.  This  work  offers  a  beginner  In  polarization 
studies  a  concrete  foundation  Into  the  basic  theory  of  polarization.  We  also 
refer  to  Gerard  and  Burch  (1975)  In  which  the  transformations  from  the  Mueller 
matrix  to  the  Jones  matrix  representation  are  provided  In  the  appendix,  and  In 
which  the  case  of  partially  polarized  wave  Interaction  Is  treated.  A  succinct 
but  good  Introduction  for  the  radar  case  Is  given  In  R.  E.  Collin  and  F.  J. 
Zucker,  Antenna  Theory,  McGraw  Hill  Co.,  1969,  providing  a  discussion  on  how 
the  Stokes  parameters  are  useful  to  polarlmetrlc  antenna  design. 

Ellipsometrlc  measurements  for  laser  light  scattering  from  metallic 
and  dielectric  rough  surfaces  of  Identical  surface  profiles  were  carried  out 
and  reported  by  Gough  and  Boerner  (1978),  showing  that  the  difference  In 
measured  values  for  the  various  Mueller  matrix  elements  as  function  of 
incidence  angle  are  primarily  dependent  on  the  volumetric  scatter  from  the 
underburden  rather  than  due  to  surface  roughness.  We  note  that  Blanchard  et 
al.  (1981),  and  Fung  and  Eom  (1982),  came  to  the  same  conclusions. 

Or.  M-M.  Boerner  and  his  research  associates  of  the  Communications 
Laboratory  at  UIC  (January  1981),  investigated  the  basic  theory  for 
polarization  utilization  in  radar  target  reconstruction,  which  Is  presented 
and  verified  by  computer  computation.  The  study  clearly  demonstrated  that  the 
optimal  polarization  concept  Is  very  useful  In  radar  target  analysis. 

This  work  was  followed  by  a  second  report  (Boerner,  et  al.,  September 
1981),  In  which  the  objectives  of  the  analysis  carried  out  are  directed  toward 
extracting  useful  target  classification  algorithms  utilizing  experimental  and 
amplitude-plus-phase  matrix  data  measured  by  Welsbrod  and  Morgan  of 
Teledyne-MIcronetlcs  and  comparing  those  with  amplitude-only  data  obtained  In 
the  late  fifties  by  Huynen  of  Lockheed  (1960).  For  this  purpose,  the 
rudimentary  properties  required  with  its  transformation  invariants,  and  the 
optimal  polarization  descriptors,  are  summarled  In  the  previous  report 
(Boerner,  et  al.,  June  15,  1981). 
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The  measurements  by  Welsbrod  and  Morgan  of  Taledyne-MIcronatlcs  were 
carried  out  for  18  buoy  targets  located  above  a  sea-bed  surface*  Two 
specific,  rather  different  buoy  shapes  are  used,  where  six  of  the  targets  are 
"4-corner"  dihedral  reflectors  with  a  yaw  angle  of  90*,  and  the  rest  of  the 
targets  are  horizontal  open  pipe  sections  with  varying  sizes  and  displacement 
height  above  sea-bed  level.  We  note  here  that  Huynen  In  his  May  I960  report 
presented  most  excellent  data  on  24  different  Isolated  missile  decoy  shapes 
measured  at  9.735  GHz. 

The  scattering  matrix  data  of  Welsbrod  and  Morgan  were  measured  at 
only  one  frequency  of  3.150  GHz  as  functions  of  the  aspect  angle  and  presented 
to  us  in  printed  and/or  tape  format.  Their  results  and  own  analyses  are 
summarized  In  (Morgan  and  Welsbrod,  "RCS  Matrix  Studies  of  Sea  Clutter",  Final 
Report,  March  1982),  which  are  reproduced  here  with  their  consent:  (1)  In 
general,  RCS  matrix  response  of  a  target  mounted  over  water  will  be 
significantly  different  from  those  of  free  space.  This  effect  Is  attributed 
to  Illumination  taper  caused  by  Interference  between  the  direct  and  the 
reflected  components.  (2)  Of  various  targ  ts  which  were  put  together,  only 
the  dihedral  was  readily  Identifiable  under  most  conditions.  Other  targets, 
which  were  primarily  cylinders  mounted  on  a  vertical  buoy,  were  Identifiable 
part  of  the  time,  but  the  buoy  support  structure  tended  to  dominate  and  mask 
inhere  signatures  of  other  targets.  (3)  Single  parameter  results  as 
l°llJ  I°rr  1  **  I°lrI  turned  out  quite  effective  in  separating  a 
dihedral  from  any  other  simple  target.  A  flat  plate  type  of  target  could  also 
be  successfully  separated  from  other  targets  using  a  single  parameter 
classifier.  Utilizing  single  parameter  classifiers  and  sifting,  it  was  found 
to  be  possible  to  classify  targets  as  dihedrals,  flat  plates,,  or  neither  with 
about  70  to  80  percent  probability  of  being  correct.  (4)  When  used  alone,  the 
single  *^equ?  RCS  matrix  Is  primarily  useful  In  classifying  objects  whose 
seafaring  oxers  are  separated  by  a  small  fraction  of  a  wavelength.  Larger, 
more  complex  objects  should,  however,  be  classifiable  using  multifrequency  RCS 
matrix  data.  This  can  be  done  by  breaking  the  complex  target  Into  simpler 
elements  throui  ansformatlon  to  the  time  domain  and  Identifying  the 
components.  amatively,  the  signature  In  the  frequency  domain  could  be 
matched  to  a  library  of  reference  signatures  with  the  different  elements  of 
the  RCS  matrix  being  used  to  multiply  the  dimensions. 
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In  our  NAV-A1R  effort  (Boerner,  et  *1.,  September  15,  1981),  we 
developed  computer-assisted  program  to  determine  the  opt  1m  1  pol  aril  at  Ion 
characteristics,  ana  we  compered  the  Measured  RCS  with  that  obtained  fro* 
approx 1 eate  scattering  formulae  developed  for  these  shapes.  All  of  our 
results  are  compared  with  those  obtained  by  0.  Richard  Huynen  In  the  late 
fifties  and  as  reported  In  his  Nay  1960  report,  and  with  results  reported  In 
previous  reports  by  Morgan  and  Uelsbrod  (1960,  1981,  1982).  Our  specific 
novel  contributions  to  the  Interpretations  of  the  results  of  Uelsbrod  and 
Morgan  (1982)  are  contained  In  the  N.Sc.  thesis; 

B-Y.  Foo,  "A  High  Frequency  Inverse  Scattering  Model 
to  Recover  the  Spectral  Point  Curvature  froai  Polarlmetrlc 
Scatter  Data",  CL-EM ID-82-02,  M.Sc.  thesis.  College 
of  Graduate  Studies,  UICC,  May  1982. 

In  another  very  thoughtful  and  excellent  M.Sc.  thesis, 

C-Y.  Chan,  "Studies  on  the  Power  Scattering  Matrix  of  Radar 
Targets,  Cl-EMID-81-02,  M.Sc.  thesis,  College  of 
Graduate  Studies,  UICC,  May  1982. 

various  complete  poWimetric  measurement  techniques  as  compared  to  the 
Incomplete  techniques  proposed  In  Huynen  (1960,  1970,  1978), 
McCormick/Allan/Hendry  (1975-1981),  Root  (1980)  are  identified.  These  Initial 
studies  require  Immediate  advancement. 

TARGET  VECTOR  SCATTERING  CENTER  ANALYSIS 
Scalar  Case 


The  scalar  scattering  center  discrimination  technique  was  derived 
from  geometrical  optics  considerations  to  a  large  extent  at  the  Cornell 
Aeronautical  Laboratories  (CAL:  Kell,  1965;  Bechtel  and  Ross,  1966;  Hammer, 
1967),  and  It  draws  heavily  from  the  monostatlc-bl static  equivalence  theorem 
derived  by  Kell  (1965)  and  extended  to  the  general  bistatic  case  In  (Blckel, 
1965).  Starting  from  the  Stratton-Chu  vector  diffraction  Integral  (Stratton, 
1941),  Kell  (1965)  derives  the  "blstatlc  scattering  Integral* 

O  -  ^T(z)  exp[J2k0z  cos(p/2)]  dll2 

where  T(z),  not  precisely  known.  Is  a  composite  expression  of  surface  geometry, 
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of  illuminating  and  observing  ray  geometry,  and  of  surface/ creeping  wave 
propagation  effects.  A  cylindrical  coordinate  system  (a,  e,  z)  was  assumed  to 
relate  nonostatic  and  bistatic  properties  such  that  Its  z-axls  Is  coincident 
with  the  bisector  of  the  bistatic  angle  i  formed  by  the  transmitter,  target 
centroid,  and  receiver.  It  Is  found  that  the  above  Integral  can  be  subdivided 
Into  a  sum  of  Integrals,  each  of  which  Is  taken  over  the  nnge  of  z  within 
which  Its  Integrand  Is  continuous  and  contributions  only  come  at  end  points  of 
each  subregion  Integral  which  define  the  scattering  centers.  The  contribution 
of  each  scattering  center  to  the  total  bistatic  Integral  corresponds  to  the 
neighborhood  of  the  scattering  center  over  which  the  chase  of  the  net  phase  of 
the  Integrand  In  the  above  equation  remains  with  «/2\  of  Its  value  at  the 
scattering  center  and  which  applies  to  both  monostat1c\and  bistatic  scattering. 
Phase  changes  are  Indicated  by  the  factor  kQ  cos(e/2),  and  It  follows  that 
small  bistatic  angles  (less  than  10*)  have  little  effect  on  wavelength,  whereas, 
at  larger  bistatic  angles  the  reradiation  characteristics  of  the  scattering 
centers  are  most  Important. 

Evaluation  of  the  above  equation  over  the  various  separated  scattering 
centers  or  "cophased  areas"  yields  the  total  RCS 

H 

•  -  iX)  J*-i2 

m-1 

with 

•a-  2k„  ».  c« 

where  oR,  #R,  z|R,  and  ^  denote  the  RCS.  Individual  phase, 
distance  to  first  phase  center,  and  residual  phase  contribution  of  the  mth 
center  Including  creeping  wave  path-length  phase  corrections,  where  /a^t 
zR,  and  Cn  are  Insensitive  to  the  bistatic  angle  •  over  the  range  of 
e  considered,  for  those  centers  which  are  significant  In  «  ,  It  then 
follows  that  the  bistatic  cross-section  of  aspect  angle  o  and  bistatic 
angle  a  Is  equal  to  the* monostatic  RCS  measured  on  the  bisector  at  a 
frequency  lower  by  a  factor  cos  >/2  . 
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Using  a  narrowband  reconstruction  approach.  It  was  shown  molnty  by 
CAL  (Homer,  1967)  how  the  scattering  centers  can  be  determined  and  allocated 
from  Lack-scattered  data.  Using  principles  of  synthetic  aperture  radar 
(Harger,  1970),  It  Is  shown  In  (Graf,  1972)  and  In  (Gnlss,  Magura,  i976)  how 
Doppler  phase  Information  can  be  used  to  Identify  scattering  centers  of 
rotating  bodies.  It  should  be  noted  that  narrowband  reconstruction  requires 
considerably  Increased  efforts  In  data  processing  as  compared  to  wideband 
short-pulse  techniques  In  which  rase  the  complexity  can  be  built  Into  the 
radar  Itself  (Ross,  1978),  though  both  narrowband  (MITRE,  CAL,  RRE-UK-Malvern, 
RCA-Hoorestown,  NYU)  and  wideband  (Westlnghouse-Baltlmore,  SURC/Rome  APB, 
MIT-LIncoln  Labs/Kwajeleln,  and  Pattern  Recogn.  Inc. /Rome)  systems  are 
required  simultaneously  In  efficient  pattern  recognition  schemes  (Gogglns, 
1978;  Kslenski,  et  al.»  1973;  Repair,  1970).  It  should  be  noted  that 
narrowband  Interferometric  methods  (Tomljanovlch,  et  al.,  1968)  have  been 
developed  also  In  other  fields  such  as  aperture  synthesis  (Harger,  1970), 
speckle  Interferometry  (Dainty,  1975)  and  In  synthetic  interferometer  radar 
(Porcello  and  Allan,  1974) . 

VECTOR  CASE:  HUYNEN'S  N-TARGI.T  DECOMPOSITION  THEORY 

Of  particular  relevance  to  this  electromagnetic  target  scattering 
problem  is  the  Interaction  of  polarization/depolarization  sensitive  scattering 
centers  on  a  single  closed  target  of  Irregular  shape,  which  was  first  attacked 
rigorously  by  Huynen  (i960)  In  his  dissertation.  In  this  masterpiece,  he 
developeo  his  little  understood  "N- tar get  decomposition  theorem",  utilizing 
canonical  properties  of  Kennaugh's  optimal  target  polarization  null  theory, 
which  specifically  applies  to  clutter  analysis  and  multiple  vector  scattering 
center  interaction  of  single  targets.  This  theory  Is  of  paramount  Importance 
to  further  advancement  In  radar  polarlmetry.  Although  it  still  requires 
extensive  extension.  It  clearly  paves  the  way  to  a  single  unique  method  of 
complete  polarlmetrlc  description  of  an  extended  target  In  clutter,  an 
Immensely  complicated  electromagnetic  Inverse  problem. 


In  collaboration  with  Morgan  and  Welsbrod  (Teledyne-MIcronetlcs),  In 
our  studies  of  polarlmetrlc  CW  radar  target  characteristics*  we  are  following 
the  direction  of  examining  a  set  of  simple  canonical  target  shapes  such  as  the 
sphere,  the  linear  wire  target,  the  n-bounce  corner  reflector,  the  left/right 
winding  helices,  the  cone-tlp/oglval  and/or  spherical  capped  cylinders  with 
and  without  fins,  bumps,  protrusions,  etc.  These  are  treated  In  (Boemer,  et 
al.„  January  15/September  15,  1981).  Second,  In  consultation  with  Kleras  and 
Bennett,  Sperry  Research  Center,  (see  present  RADC  contract  report),  a  CW 
vector  dumbell  scattering  center  (matrix)  Interaction  was  chosen.  Both 
methods  have  proven  to  provide  useful  results  and  can  be  used  for 
Interpretation  of  the  motion  of  the  Huynen  polarization  fork  as  a  function  of 
frequency,  relative  aspect  angle  (with  respect  to  the  line  joining  the  vector 
scattering  centers)  and  the  electric  separation  of  component  scattering 
centers. 

DYNAMICS  Of  POLARIZATION  FORK 


Specifically,  for  the  dumbell  target,  we  observe  that  for  llrtear 
(H,  V)  polarization  base  pair  anchoring,  the  cross-polarization  nulls  move 
only  whenever  the  principal  target  symmetry  axis  Is  rotated  about  the  line  of 
sight  orthogonal  to  (H,  V);  and  that  the  co-pol arlzatlon  null  locations  move 
on  quasi-circular  spiral  non-closing  paths  as  functions  of  differential  change 
In  aspect  angle:  for  small  differences  In  component  scattering  center 
strengths,  the  circles  remain  within  Isolated  patches,  whereas  for  large 
differences,  the  nulls  move  on  large  circles  encircling  the  total  polarization 
sphere.  Furthermore,  the  electric  separation  of  the  scattering  centers 
dictates  the  relative  differential  speed  with  which  these  loci  are  traversed 
(slow  versus  rapid)  as  functions  of  differential  aspect  angle.  Me  also  note 
that  for  a  large  ensemble  of  closely  packed  vector  scattering  centers,  the 
loci  of  the  co-pol arlzatlon  nulls  remain  within  rather  small  Isolated  patches 
on  the  polarization  sphere  which  Is  Indicative  of  clutter-type.  Furthermore, 
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the  specific  quasi-circular  paths  drawn  are  Indicative  of  clutter  notion.  Me 
note  that  the  analytical  result  was  verified  experimentally  by  Poelman 
(1980-62)  as  explained  in  detail  In  (Boerner,  STAS  1914,  September  30,  1961), 
and  this  specific  phenomenon  of  the  dynamic  fork  motion  of  time  frames  of 
below  the  vector  scattering  center  reshuffling  time  requires  further  extensive 
analytical  and  experimental  studies.  In  general,  the  cross-polarization  null 
location  for  linear  symmetric  targets  Is  of  slow  precession  type,  and  the 
rapid  quasi-circular  path  motion  of  the  co-pol arlzatlon  null  location  Is 
nutatlve  gyroscopic  In  nature.  We  note  that  this  specific  dynamic 
polarization  fork  behavior  Is  well  described  by  Huynen's  "single  target" 
decomposition  Into  five  target  characteristic  parameters  (pm,  t*  v, 

However,  the  electromagnetic  Inverse  problem  of  decomposing  a  single 
radar  target  into  Its  characteristic  polarlmetrtc  target  vector  scattering 
centers  Is  very  complicated  and  still  not  resolved.  In  a  joint  research 
effort  with  ESL-OSU,  EML-UIU,  SRC,  BAC,  we  at  the  EMID-CL-UIC  are  planning  to 
attack  this  problem  In  the  near  future. 
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